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Kurzfassung 


Weltweit werden ca. 1% der gesamten elektrischen Energie für die Chlor-Alkali Elektrolyse 
aufgewendet. Die für die Elektrolyse benótigte elektrische Energie kann um ca. ein Drittel 
reduziert werden, indem auf Silberkatalysator basiertende Gasdiffusionselektroden — sog. Sauer- 
stoffverzehrkathoden (englisch: oxygen depolarized cathodes, ODC) - in den Prozess integriert 
werden. Ungeachtet ihres großen Energieeinsparpotentials sind die Prozesse und Limitierungen 
innerhalb der ODC weder verstanden noch quantifiziert. Die Einflüsse der Elektrodenstruktur auf 
die Leistungsfähigkeit und die Dynamik der ODC sind kaum ergründet. Zur Beantwortung der 
genannten Forschungsfragen wird das erste Modell einer ODC vorgestellt, das das Phasengle- 
ichgewicht der Gas-Flüssigkeitsgrenzfläche und strukturbedingte Inhomogenitäten in der Elek- 
trolytverteilung berücksichtigt. 

Es wird gezeigt, dass die Stromdichte durch die Verfügbarkeit an gelóstem Sauerstoff im flüssi- 
gen Elektrolyten limitiert wird. Bereits im stromlosen Zustand ist nur wenig Sauerstoff in dem 
Elektrolyten lósbar, der Grund liegt in der hohen Ionenkonzentration und der damit verbundenen 
niedrigen Wasseraktivität. Durch die Sauerstoffreduktion wird mit zunehmender Stromstärke 
Wasser verbraucht und Hydroxidionen produziert, die sich nahe der Gas-Flüssigkeitsgrenze an- 
reichern. Dies führt zu einer weiteren Verringerung der Wasseraktivität und somit zu einer weit- 
eren Verschlechterung der Sauerstofflöslichkeit und schließlich zu einer Sauerstoffverarmung. 
Auf Grundlage dieser Erkenntnisse wird im zweiten Teil der Einfluss des Wassermassentrans- 
ports und der Wasseraktivität detaillierter untersucht. Dazu werden zwei vergleichbare ODCs 
in Messaufbauten mit unterschiedlichen Stofftransportbedingungen in der Flüssigphase mithilfe 
des dynamischen Dreiphasenmodells analysiert. Es wird gezeigt, dass ein konvektiver Elek- 
trolytstrom nicht nur zu einer hóheren Elektrodenperformance, sondern auch zu einer weitaus 
schnelleren Dynamik führt. Beide Effekte sind auf kürzere Diffusionsschichten in der Flüssig- 
phase und somit auf einen schnelleren Wasser- und Ionenmassentransport zurückzuführen. Die 
Wasseraktivität beeinflusst die Sauerstoffreduktion nur geringfügig direkt, aber signifikant indi- 
rekt über die Sauerstofflóslichkeit. 

Im dritten Teil wird der Einfluss der ODC-Struktur auf die Performance und Dynamik un- 
tersucht. Ein modelgestützter Vergleich von ODCs mit unterschiedlichen Binder-Katalysator- 
Verhältnissen mittels elektrochemischer Impedanzspektroskopie führt zu der Erkenntnis, dass 
Unterschiede in der Performance im Wesentlichen auf die unterschiedlich ausgeprägte Gas- 
Flüssigkeitsgrenzschicht zurück zu führen sind, die die Menge an gelóstem Sauerstoff signifikant 
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beeinflusst. Die Impedanzspektren aller untersuchten Elektroden legen den Schluss nahe, dass 
sich die Dreiphasengrenzfläche inhomogen über die gesamte Elektrodentiefe erstreckt, jedoch 
scheinen die stärker gefluteten Elektrodenteile elektrochemisch nahezu inaktiv zu sein. 
Zusammenfassend wird mithilfe der makroskopischen dynamischen Modellierung aufgezeigt, 
dass zwei Faktoren besonders entscheidend für eine hohe ODC-Performance sind: Eine große 
Gas-Flüssigkeitsgrenzschicht, die besonders durch die ODC-Struktur beeinflusst werden kann, 
und ein schneller Wasser- und Ionenmassentransport in der Flüssigphase, der sowohl durch 
eine konvektive Strómung des Elektrolyten als auch durch die ODC-Struktur begünstigt werden 
kann. Auf der methodischen Ebene wird demonstriert, dass die dynamische Modellierung ein 
geeignetes Tool ist, um komplexe Dreiphasensysteme wie das der ODC zu analysieren. Dabei ist 
besonders die dynamische Ermittlung von Zeitkonstanten zur Separation der sich überlagernden 
Prozesse unabdingbar. 


Abstract 


Worldwide approximately 1 % of the total electrical energy is spent on chlor-alkali electrolysis. 
The demand of electrical energy for the electrolysis can be reduced by about one third by in- 
tegrating silver catalyst-based gas diffusion electrodes — so-called oxygen depolarized cathodes 
(ODC) - into the process. Regardless of its large energy saving potential, processes and limita- 
tions within ODCs are neither understood nor quantified. Also, the influences of the electrode 
structure on performance and dynamics are poorly investigated. To address the mentioned re- 
search questions, the first model of an ODC is presented that considers the phase equilibrium of 
the gas-liquid interface and structure-related inhomogeneities in electrolyte distribution. 

Itis shown that the current density is limited by the availability of dissolved oxygen in the liquid 
electrolyte. Even at open-cell state, oxygen is poorly soluble in the liquid electrolyte due to the 
high ion concentration and low water activity. With increasing current density, the oxygen reduc- 
tion reaction consumes water and produces hydroxide ions which accumulate near the gas-liquid 
interface. This leads to a further reduction in water activity and thus to a deterioration in oxygen 
solubility and finally to oxygen depletion. 

Next, the influences of water mass transport and water activity are investigated in detail. For this 
purpose, ODC measurements of two experimental setups differing in mass transport conditions 
in the liquid phase are analyzed with the dynamic three-phase model. It is demonstrated that a 
convective electrolyte flow does not only lead to a higher electrode performance, but also to far 
faster dynamics. Both effects are attributed to shorter mass transport length in the liquid phase 
and thus to faster water and ion mass transport. 

In the third part of this dissertation, the influence of ODC structure on performance and dynamics 
is examined. By means of electrochemical impedance spectroscopy it is demonstrated that the 
ODC performance is mainly governed by sizes of the gas-liquid interface, which significantly 
affects the amount of dissolved oxygen. The model-based analysis of all electrodes suggest that 
the three-phase interface extends inhomogeneously over the entire electrode depth, but the three- 
phase area more close to the gas side of the ODC is electrochemically almost inactive. 

In conclusion, macroscopic dynamic modeling reveals that two factors are particularly critical 
for a high ODC performance: A large gas-liquid interface, which can be mainly influenced by 
the ODC structure, and a fast water and ion mass transport in the liquid phase, which can be fa- 
vored by a convective flow of the electrolyte and also by the ODC structure. Methodologically it 
is demonstrated that dynamic modeling is a suitable tool to analyze complex three-phase systems 
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such as ODCs. Especially the dynamic determination of time constants for the separation of the 
overlapping processes is indispensable. 
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1 Introduction 


1.1 Motivation and purpose 


One of the greatest challenges for current and future generations is to mitigate climate change 
by significantly reducing greenhouse gas emissions. Two key levers for facing this challenge are 
the climate-friendly provision of energy and the reduction of energy consumption. This thesis 
focuses on the latter aspect: A significant share of global electricity consumption is accounted 
for chlorine production. As a basic chemical for 60 % of all chemical products [3], the world 
wide demand of chlorine has been increased from 42 million tons in 2001 [4] to 52 million tons 
in 2003 [5] and 89 million tons in 2017 [3]. Using the traditional NaCl-membrane electrolysis, 
which is the most energy-efficient of all established chlorine production processes [6], the chlo- 
rine production in 2017 corresponded to a consumption of 195800 GWh electric energy — which 
equals about | % of the global consumption of electrical energy [7] — and the emission of 82 mil- 
lion tons of CO, [3]. However, further energy optimizations of the traditional NaCl-membrane 
electrolysis are not attainable, since the savings potential of this process has been practically 
exhausted [4]. 


An thermodynamically driven increase of energy efficiency is enabled by a fundamental re- 
design of the membrane-process: The chlor-alkali electrolysis using oxygen depolarized cath- 
odes (ODC). Here, in contrast to the standard membrane-process, no hydrogen is produced as 
co-product, but oxygen is consumed at the cathode [8]. The conversion of the process reduces 
the required electrical energy by about 30 % [6]. The ODC, as the centerpiece of this technique, 
is a silver catalyst-based gas diffusion electrode (GDE). In ODCs, oxygen is provided in the gas 
phase, which dissolves into the 30 wt—% NaOH electrolyte, where it reacts at the silver surface 
with water under the acceptance of electrons to hydroxide ions [5]: 


021504: 48" ^5 40H (1.1) 


1 Parts of this chapter have been published in Róhe et al., ChemSusChem, 2019 [1] and Rohe et al., ChemElectroChem, 
2019 [2]. 
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To enable both, the liquid and the gas phase within the ODC, it consists in addition to hydrophilic 
silver particles also out of hydrophobic polytetrafluoroethylene (PTFE). This structure leads to 
a complex three-phase system, in which many processes interact (cf. fig. 1.1): Mass transport 
in the gas phase, mass transport in the liquid phase, processes at the gas-liquid interface such as 
water evaporation of the oxygen solution and the actual oxygen reduction reaction (ORR). 


o Hydrophilic 
Ag-Catalyst 
Hydrophobic 
PTFE 


Gas phase 


Electrolyte: 
10 M NaOH 


e Mass transport gas phase (2) Mass transport liquid phase 


© Processes at gas-liquid interface o Electrochemical reaction 


Figure 1.1: Schematic view on processes in oxygen depolarized cathodes. 


Irrespective of the great importance of ODCs, the processes and limitations are not well under- 
stood, especially not quantitatively. Essential properties like the local distribution of the gaseous 
and liquid phase within the electrode are unknown [9]. This missing insights harm a knowledge- 
driven improvement of ODCs. 

One approach to analyze the processes separately as well as their mutual interaction and their 
contribution to the electrode performance is through dynamic modeling. 

This approach has already been used successfully in many studies: Dynamic analysis was used 
for discriminating between reaction kinetic models e.g. for oxygen reduction [10] or methanol 
oxidation [11], for quantifying the impact of reactant transport, e.g. of oxygen through liquid 
electrolyte on ORR performance [12], for state estimation [13] or concentration sensing [14]. 


The objective of this thesis is derived from the open research questions of ODCs and the oppor- 
tunities of dynamic modeling: The goal of this dissertation is to gain deeper insights into the 
processes in ODCs for advanced chlor-alkali electrolysis in order to identify limitations and to 
provide the fundamentals for knowledge-based electrode improvements. 'To reach the goal, the 
suitable approach of dynamic modeling is used to identify and separate overlapping processes 
and limitations in ODCs. 


1.2 Scope and structure 


1.2 Scope and structure 


This dissertation documents the development and analysis of dynamic three-phase models of 
ODCs. Initially, an overview of ODCs and their processes is given based on the current liter- 
ature. Subsequently, ODC models are presented, parameterized, analyzed and extended on the 
basis of intermediate results. This dissertation is structured as follows: 

The target of the second chapter is to reveal and describe the processes as well as the process 
conditions of ODCs. After a brief overview of conventional chlor-alkali electrolysis processes, 
the chlor-alkali electrolysis with ODC and the ODC itself are discussed in detail. At the end of 
the chapter, the characteristic processes and properties of ODCs and therefore the essentials for 
ODC model development are summarized. 

In the third chapter, the processes and system properties previously classified as relevant are 
considered more in detail. The aim of this chapter is to lay the base for the actual model design. 
In detail, approaches to model the gas-liquid interface, the processes at the gas-liquid interface, 
the reaction kinetics and the mass transport are presented. In addition, suitable dynamic meth- 
ods for the analysis of ODCs are discussed. The findings of these chapters are summarized in a 
concretization of the research question and in a list of concrete model requirements. 

The fourth chapter includes the actually model design of the dynamic one-dimensional three- 
phase model, the validation as well as a model analysis. The model is parameterized using 
stationary literature data. A subsequent dynamic analysis shows that the dynamic behavior of 
the ODC is dominated by the water and ion transport in the liquid phase. This leads to the 
conclusion that in the common understanding of the ODC the influence of the water activity is 
underrepresented. The results are supported by a sensitivity analysis, which shows that especially 
parameters influencing the water and ion mass transport in the liquid phase have a high influence 
on the electrode performance. 

In the fifth chapter the role of water mass transport and water activity is analyzed in-depth. For 
this purpose the model is expanded to more accurately represent the influences of mass transport 
in the liquid phase. Two different measurement setups with different mass transport conditions 


in the liquid pase were investigated with the dynamic three-phase model by steady-state and 
dynamic data. It will be shown that a poor mass transport in the liquid phase leads not only to con- 
siderably slower dynamics, but also to a much lower performance of the ODC. It is demonstrated 
that water mass transport and water activity do not significantly influence the ORR directly, but 
indirectly via the water activity-dependent solubility of oxygen. 

Chapter six focuses on the structure of the ODC. This also includes further extensions of the 


? Measurements performed by Institute of Chemical and Electrochemical Process Engineering, Technische Universität 
Clausthal, Germany and Institute of Analytical Chemistry — Center for Electrochemical Sciences, Ruhr University 
Bochum, Germany, published in [2]. 
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model regarding inhomogeneities in electrolyte distribution. The influence of different PTFE- 
silver ratios to electrode performance and dynamics is evaluated by means of electrochemical 
impedance nd It is demonstrated that the three-phase interface extends inhomoge- 
neously over the entire electrode depth. The simulation results indicate that the more heavily 
flooded electrode parts are almost electrochemically inactive. A comparison of ODCs with dif- 
ferent structures shows that not only the size but also the position of the three-phase boundary is 
crucial for high electrode performance. 

The dissertation is completed with a conclusion in the seventh chapter, which assesses the im- 
pact of the findings and provides an outlook for future work. 


5 Measurements performed by Institute of Chemical and Electrochemical Process Engineering, Technische Universitit 


Clausthal, Germany, published in [15]. 


2 Advanced Chlor-Alkali Electrolysis 
and Oxygen Depolarized Cathodes' 


In this chapter, the processes within ODCs as well as the entire chlorine-alkali electrolysis are 
described based on the present literature. The aim of this chapter is to gather information about 
the system, the operating conditions and the processes taking place. This is essential to develop 
the model and to understand the interaction of processes, and the influence of process conditions. 
In the following, an overview of different chlor-alkali electrolysis processes — including the ad- 
vanced electrolysis with ODCs — is given. Subsequently, the processes in ODCs are considered 
in detail. The essential findings of this review for the model design are summarized at the end of 
this chapter. 


2.1 Processes of chlor-alkali electrolysis 


In all conventional chlor-alkali electrolysis processes water and sodium chloride (NaCl) are used 
as reactants and chlorine (Cl5) and caustic soda (NaOH) are obtained as products. Additionally, 
hydrogen is produced as a byproduct: 


5 
Where Eeen is the open circuit voltage (OCV) of the full cell under operating conditions. In the 
novel energy-saving membrane process with ODCS, this reaction is changed. While oxygen is 
consumed as an additional reactant, hydrogen is not produced: 


2 NaCl -+H,0 +0.50, > 2NaOH + Cl, f Ecen = 1.01 V (2.2) 


^ Parts of this chapter have been published in Róhe et al., ChemSusChem, 2019 [1] and Róhe et al., ChemElectroChem, 
2019 [2]. 

Voltages and potentials given in this chapter are referred to the operating conditions of the membrane process with 
and without ODC. Depending on the source, the values differ slightly due to varying process conditions. The data 
are taken from [6]. 
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Due to the reduced cell voltage of this innovative process, the need for electrical energy for 
chlor-alkali electrolysis is reduced by about 30 % — in comparison to the most energy-efficient 
conventional process [6]. As discussed later, the energy savings due to more favorable thermo- 
dynamics are partly counteracted by the kinetic losses of the ODC. 

In the following sub-sections, the conventional methods are introduced shortly, in-depth descrip- 
tions are given by Schmittinger et al [16]. Afterwards, the membrane process using ODCs is 
explained more in detail. 


2.1.1 Conventional methods 


In all three conventional processes — namely, the diaphragm process, the amalgam process and 
the membrane process — sodium chloride in an aqueous solution is used as the anolyte, where the 
NaCl is dissociated into its ions: 


NaCl = Na* + Cr” (2.3) 
As the actual anode reaction, chloride ions are oxidized to gaseous chloride: 
2Cr CbTt-c2e Eo = 1.33 V vs. SHE (2.4) 


Eo is the open circuit potential (OCP) of the half cell reaction. 
At the cathode, the kinetically preferred reaction is the hydrogen evolution reaction (HER), 
wherein water is reduced to hydroxide ions and hydrogen [4]: 


2H,0+2e >H,f+20H Eo = —0.90V vs. SHE (2.5) 


In case that produced hydroxide ions get to the anode side, the formation of hypochlorite (CIO , 
eq. (2.6)) and chlorate (CIO;, eq. (2.7)) take place as undesired side reactions: 


Cl 20H -OO C4 HO Q.6) 


6CIO +3H,0 ^ 2CIO; +4Cl +6H*4+1.50)+6e (2.7) 


These side reactions would result in contamination of the product stream and in a reduced elec- 
tricity yield [4]. The differences between the common chlor-alkali electrolysis techniques and 
the respective cathode reactions are mainly characterized by the way to ensure the absence of 
hydroxide ions in the anode chamber and thus to prevent the side reactions (2.6) and (2.7) [6]. 


2.1 Processes of chlor-alkali electrolysis 


In the diaphragm process, both half cells are separated by a porous diaphragm, which tradition- 
ally consists of asbestos. In newer applications, also PTFE and other inorganic substances are 
added to the asbestos [5]. At the cathode side, HER (eq. (2.5)) takes place in a caustic soda 
catholyte. To avoid hydroxide ions being transported to the anode, a convective flow of the 
anolyte is driven through the diaphragm to the cathode side. An advantage of this process is a 
relative low electrical energy demand of 2300 to 2900 kWh fci [6]. The main disadvantage is 
the sodium chloride contamination of the caustic soda solution [4]. 


In the amalgam process — also called mercury process — the alkaline HER and thus the formation 
of hydroxide ions is suppressed by the addition of mercury to the electrolyte. No separation 
between the anode and the cathode chamber is necessary. However, this process requires two 
separated electrochemical cells, the actual electrolysis cell and an amalgam decomposer. In the 
electrolysis cell the known reaction (2.4) takes place at the anode side. At the cathode, the Na* 
ions react with mercury to amalgam, no OH is formed [5]: 


Na*+e +Hg, — NaHg, (2.8) 


In the downstream decomposer, the mercury is recovered by electrochemical decomposition of 
the amalgam: 


2 NaHg, +2H,O — 2 NaOH +H; f + 2x Hg (2.9) 


This process requires with 3100 to 3400 KWhta, an increased amount of electrical energy [6], 
but a high purity of the caustic soda is obtained [4]. 


The membrane process combines the advantage of the low electrical energy demand (2100 to 
2600 kWhtg [6]) of the diaphragm process and the high caustic soda purity of the amalgam 
process. In addition, no environmentally harmful substances such as asbestos or mercury are 
required [4]. 
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Figure 2.1: Schematic illustration of the chlor-alkali electrolysis membrane processes, a) without and b) with ODC, 
based on [3, 6]. 


As shown in fig. 2.1 a), the anode and the cathode side are separated by a membrane that is 
permeable to positive charged Na* but not to negative charged CI’ and OH”. Since the mem- 
brane must meet the requirements of long-term stability at high electrolyte concentrations and 
high temperatures, this process was a theoretical concept until about 1975; the realization was 
only made possible by the development of the Nafion™ membrane. To the anode chamber, a 
concentrated NaCl solution is continuously added. While the CI” reacts to chlorine gas (cf. eq. 
(2.4)), the Na* ions diffuse through the membrane into the cathode chamber. To the cathode side, 
a diluted sodium hydroxide solution is fed. At the cathode itself, the alkaline HER takes place 
(cf. eq. (2.5). While the produced hydrogen is outgassing, the hydroxide ions are removed 
together with the sodium ions as concentrated caustic soda solution [6]. 


2.1.2 Membrane process with oxygen depolarized cathodes 


The demand of electrical energy could be further reduced by revolutionizing the membrane pro- 
cess by utilization of ODCs [6]. Even though the advanced chlor-alkali electrolysis using ODCs 
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was already theoretically presented by Butler in 1950 [17], the first plant in pilot scale was just 
set up in the 1990s [18]. Currently, it is already used industrially, inter alia in Germany and China 
[3]. As illustrated in fig. 2.1, compared to the conventional process, the cathode is replaced by an 
ODC, which is fed with the liquid electrolyte from the one side and with gaseous oxygen from 
the other side. The cathode reaction is changed from HER to ORR in alkaline milieu, where 
water reacts with oxygen to hydroxide ions: 


0,+2H,0+4e — 40H Ey = 0.32 V vs. SHE (2.10) 


The advantage of this process lies in the cathodic cell potential: While in the conventional process 
the HER (eq. (2.5)) has a OCP of about —0.90 V vs. Standard hydrogen electrode (SHE), the 
reversible cell potential is considerably higher for the oxygen reduction (eq. (2.10)) with 0.32 V 
vs. SHE. Although the overpotential of ORR is > 0.7 V [18] and higher than for HER, under 
industrial operation conditions the total cell voltage can be reduced from about 3 V to about 2 V 
and therefore by about 30 %. Accordingly, the demand of electrical energy is reduced to about 
1470 to 1820kWh tor: However, when considering the total energy consumption, it must be 
taken into account that oxygen with high purity is additionally required and the production of 
hydrogen is omitted [4]. 

A sustainability analysis of this process is presented by Jung et al. [19], deeper insights into the 
construction of industrial electrolyzers are given by Kintrup et al. [3]. An extensive overview of 
chlor-alkali electrolysis with ODC literature until 2008 is presented by Moussallem et al. [6]. 


2.2 Processes in oxygen depolarized cathodes 


As illustrated in fig. 2.1 b), the ODC in an industrial chlor-alkali electrolyzer is fed with gaseous 
oxygen from the one side and with liquid electrolyte from the other side. At the gas-liquid 
interface within the ODC, oxygen dissolves into the liquid electrolyte and is further transported 
to the catalyst surface where the ORR (cf. eq. (2.10)) takes place. To enable the industrially 
required current densities of 4 to 6kA m? [20], various demands must be ensured [18]: 


* A catalytic active and electronically contacted electrode surface 
* A high oxygen partial pressure in the gas phase 
* A sufficient amount of oxygen dissolves into the electrolyte 


* A satisfactory transport of water, oxygen and OH in the liquid phase 
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As a catalyst, typically silver is used. Silver is electronically highly conductive [18] and has a 
high catalytic activity for the ORR that is comparable to platinum when operating under indus- 
trial conditions [21] of NaOH-electrolyte concentrations exceeding 10 M and temperatures of 80 
to 90°C [6]. 

However, under these conditions the solubility of oxygen into the liquid electrolyte is poor [22]. 
To fulfill the criterion of a high dissolution rate of oxygen, the poor solubility must be compen- 
sated by a large size of the gas-liquid interface [18]. For a large interface between both phases, 
the proportions of pores filled with gas and pores filled with liquid electrolyte must be tared out 
well. To accomplish this, hydrophobic PTFE is used as second component in ODCs in addition 
to the hydrophilic silver. How different PTFE and silver loading influence the performance of an 
ODC is discussed by Franzen et al. [23] and Moussallem et al. [9]. 

The dissolution rate of oxygen into liquid electrolyte depends also on the partial pressure of 
oxygen in the gas phase at the gas-liquid interface. A high oxygen partial pressure is achieved 
by a high mass transport rate of oxygen in the gas phase. This implies that the evaporated water 
has to be removed through the gas-filled pores of the ODC [18]. 

In the liquid phase, the dissolved oxygen has to be transported from the gas-liquid interface to the 
catalyst surface. The mass transport resistance of this process is minimized, when the electrolyte 
forms a thin film of maximum a few micrometers on the catalyst surface. However, the nature of 
the liquid phase must also enable a sufficient transport of water as a reactant and hydroxide ions 
as the reaction product between the reactive zone and the liquid bulk phase [18]. 


At locations within the ODC where the requirements of a wetted, electrically connected catalyst, 
a nearby gas-liquid interface and sufficient mass transport in the gas and liquid phase are fulfilled, 
the ORR can take place. Although the reaction kinetics of ORR in the alkaline environment with 
a silver catalyst has been the subject of many studies [24—30], the exact reaction mechanisms 
are not clear [30], and various — partly contradictory — reaction paths are proposed. Several 
studies indicate that hydrogen peroxide may be formed as an intermediate reaction product [24— 
26]. But in general, silver is known to have a high catalytic activity towards peroxide reduction 
[31], so that a complete reduction of the oxygen and no hydrogen peroxide as a final product 
can be assumed. However, different reported reaction pathways coincide in the formation of 
the superoxide as the rate determining step [24—29], independent of pH value [26] or the crystal 
structure of the silver catalyst [25]: 


Oe 305 (2.11) 


Pinnow et al. gave the restriction that for ODCs under industrial conditions this assumption is 
not valid for cathodic overpotentials below 0.155 V [32]. 
Exemplary for a possible reaction path, the results of Ignaczak et al. [28] are discussed in the 
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following. Please note, even though this is the result of an extensive study, this mechanism is 
not validated for the specific process conditions of ODCs. In their recently released, comprehen- 
sive study combining density functional theory (DFT) with electrocatalysis theory and electron 
transfer theory, Ignaczak et al. postulats that the ORR in alkaline media on silver catalyst is 
proceeding in the steps below [28]: 


O; +H,0+e = H0; - OH (2.12) 
HO; +H,0 = 20H,4-4- OH- (2.13) 
20H,4+2e =20H- (2.14) 


According to the authors, the reaction steps (2.11) and (2.12) take place in the so-called outer- 
sphere mode [28], which means that the reactants do not adsorb on the catalyst, but are bound 
to the catalyst via a bridge of adsorbed OH and hydrogen bonds [33]. The O—O bond of the 
oxygen is split in the reaction step 2.13. For an in-depth discussion of the microkinetics, it is 
referred to the literature, e.g. Ignaezak et al. [28] or Ge et al. [30]. 

Further aspects of ODCs were inter alia studied by Botz et al., who measured an increasing ion 
concentration during ODC operation at the electrode surface [34]. Zhang et al. evaluated the 
electrochemistry of oxygen in concentrated NaOH-electrolyte [27]. 


2.3 Concluding remarks 


In this chapter, literature about chlor-alkali electrolysis with ODCs and the ODC itself was re- 
viewed. It was possible to gain insights into the electrode structure, the processes taking place as 
well as the process conditions. 

Notwithstanding their great importance, ODCs as well as the interaction of processes and limi- 
tations in ODCs are not clearly understood: Properties like size, location and distribution of the 
gas-liquid interface are not identified yet, and the time constants of the different processes remain 
unknown. Also, it is not revealed yet, if the equilibrium of the gas-liquid interface stays constant 
during operation. This hinders knowledge-based improvement of ODCs. A promising tool to 
close these knowledge gaps is the development and use of dynamic models. As shown in various 
studies before, dynamic simulation gain information on the importance and interaction of mass 
transport, reaction mechanism and kinetics or geometrical structures. This is difficult to access 
in steady state simulation. Examples for complex model-based dynamic analysis are the identifi- 
cation of reaction kinetics [35] or the discrimination between charge transfer and mass transport 
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losses [36] using electrochemical impedance spectroscopy (EIS). A notable, tangible alternative 


is given by the evaluation of potential steps (multi-step chronoamperometry), chronoamperome- 


try or potentiometry. These methods have been successfully used inter alia for identification of 


the transport and reaction processes in direct methanol fuel cells (DMFC) [37] or to investigate 


limiting factors for the ORR in metal-air batteries [12]. 


Based on the discussion in this chapter, the following points must be considered for the develop- 
ment of an ODC model: 


The model must be valid for typical industrial conditions such as high electrolyte con- 
centrations exceeding 10M, temperatures of 80 to 90°C and current densities of 4 to 
6kAm ?. 


The model structure must reproduce a realistic gas-liquid interface. 
The phase transitions have to be considered, i.e. oxygen dissolution and water evaporation. 
In the gas phase, oxygen and water are the mainly transported species. 


In the liquid phase, the characteristic mass transport path of oxygen is from the gas-liquid 
interface to the catalyst surface. The characteristic mass transport path for water and ions 
is the distance between the reactive zone of the ODC and the liquid bulk. 


As various literature sources point out, the rate determining reaction step for the alkaline 
ORR is the formation of the superoxide (cf. eq. (2.11)). It must be verified whether this is 
also the case for the specific conditions of ODCs. 


In order to develop a suitable ODC model, the above-mentioned aspects must be considered by 


choosing a feasible macroscopic modeling approach and mathematical-physical equations. That 


is the content of the following chapter. 
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3 Dynamic Modeling and Analyzing of 
Oxygen Depolarized Cathodes 


In the previous chapter the characteristic processes in and properties of ODCs were discussed. 
It has been demonstrated that mass transport in the gas and the liquid phase, the dissolution 
of oxygen and evaporation of water are crucial factors for the performance of an ODC. For a 
knowledge-based improvement of ODCs, the influences as well as the interactions of these pro- 
cesses have to be understood, and the performance driver and limiting factors have to be identi- 
fied. A suitable tool for this purpose is given by the dynamic modeling and analysis of ODCs. 
In this chapter, it is demonstrated how the mentioned processes, their interactions as well as 
their contribution to the ODC performance can be analyzed with the help of dynamic macro- 
scopic modeling. It is shown how the gas-liquid interface within the ODC can be represented 
in a macroscopic model and how ODCs have been modeled in previous studies using stationary 
methods. In the further course of the chapter, mathematical-physical approaches for modeling 
oxygen dissolution, water evaporation, reaction kinetics as well as mass transport in the liquid 
and gas phase are worked out specifically for the ODC. In addition, different dynamic elec- 
trochemical analysis techniques are presented which are suitable to evaluate the ODC and its 
processes by the aid of the macroscopic model. In the end, the chapter is concluded and an out- 
line is given, which research questions will be answered in the further course of this dissertation 
using the dynamic model. 


3.1 Mechanistic and literature models of oxygen 
depolarized cathodes 


In this section, different mechanistic approaches to reproduce the gas-liquid interface of three- 
phase systems in porous structures are presented first. Second, an overview is given, how ODC 
have been modeled in previously publications. 
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3.1.1 Mechanistic models of gas diffusion electrodes 


Different approaches of macroscopic modeling of porous GDE are reviewed by Yuh and Selman 
[38] or by Pinnow [39]. In the following the most important approaches are summarized and 
evaluated with regard to their ability to describe the gas-liquid interface sufficiently. 


Single pores 


Agglomerates 


L] Liquid electrolyte L] Electrode material 
L] Gas phase ee Gas-liquid interface 


Figure 3.1: Schematic illustration of different model approaches for three-phase systems within a porous GDE: a) single- 
pore model (flooded), based on [38]; b) meniscus-model, based on [38]; c) thin-film model, based on [38]; d) model with 
micro and macro pores, based on [40]; e), flooded agglomerates model, based on [41] f); thin-film flooded agglomerate 
(TFFA) model, based on [42]. 


The first modeling approach is the simple-pore model, where a plane gas-liquid interface is as- 
sumed to be within or at the edge of a single pore, see fig. 3.1 a). That approach is characterized 
by a low gas-liquid interface which goes ahead with mass transport limitations at low current 
densities [43]. Extensions of this approach are the meniscus (fig. 3.1 b)) and the thin-film model 
(fig. 3.1 c)). In the meniscus model, it is assumed that the size of the gas-liquid interface within 
a pore is defined by the contact angle between the electrolyte and the pore wall. Dependent on 
the surface tension the meniscus is formed, which goes ahead with a higher gas-liquid interface 
[44]. In the thin-film model, the approach of an enlarge gas-liquid interface is taken further: 
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While the pore surface is wetted with a liquid thin-film, the inner pore space is filled with gas. 
Due to the enlarged gas-liquid interface, higher and more realistic limiting current densities can 
be modeled [43]. Typical values for the thickness of the thin-films are between 25 and 300 nm 
[45]. Burshtein et al. have presented a model which considers micro and macro pores (fig. 3.1 
d)). While the larger macroscopic pores follow the thin-film approach, the smaller microscopic 
pores are filled with the liquid electrolyte. In contrast to the pure thin-film model, this approach 
is also able to satisfactorily model ionic mass transport and electrolyte conductivity. [40] 
Separately, Giner and Hunter have developed an approach for a PTFE-bonded GDEs by chang- 
ing the focus from single pores to larger structures (fig. 3.1 e)). According to their assumptions, 
the combination of a more hydrophilic catalyst and the hydrophobic PTFE creates macroscopic 
structures. These are represented by cylindrical flooded macroscopic pore networks of "flooded 
agglomerates" extending in the gas phase of the GDE. The flooded agglomerates themselves 
consist of microscopic pores filled with electrolyte. This approach has made it possible to sim- 
ulate high ionic conductivities and current densities j > 4kAm ? [41]. In order to reproduce 
mass transport more accurately, Cutlip has combined the thin-film with the flooded agglomer- 
ate approach to develop the thin-film and flooded agglomerate (TFFA) model (fig. 3.1 f)). The 
TFFA model considers flooded agglomerates within the porous GDE, which are additionally sur- 
rounded by a thin-film. The evaluation of the model has shown that the limiting current density is 
almost independent of the electrode thickness and the assumed shape of the flooded agglomerates 
[42]. 


3.1.2 Models of oxygen depolarized cathodes for chlor-alkali 
electrolysis 


In literature, a high number of electrochemical models for oxygen reduction reaction in GDEs 
is known [46], but only a few address the specific process conditions of ODCs for chlor-alkali 
electrolysis. To the best of the author's knowledge, all these ODC models are stationary. 

In 1997, Wang and Koda introduced their stationary model of an ODC as a TFFA approach with 
cylindrical flooded agglomerates [47]. In terms of mass transport, this model focuses exclusively 
on oxygen, while water and ions are not considered. The oxygen mass transport in the gas 
phase is calculated by Knudsen diffusion and in the thin-film as well as in the liquid phase by 
Fick's diffusion. The reaction kinetics of the ORR are modeled as a direct 4-e -process and are 
implemented by an Tafel approach. The results demonstrate how a variation of different model 
parameters affects the polarization of ODCs. For example, a thicker thin-film leads to a lower 
limiting current density or a lower porosity leads to higher overpotentials independent of the 
current density [47]. The work of Wang and Koda was adapted by Sudoh et al. whereby the 
essential characteristics of the model remained unchanged [48]. By the means of the model, 
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experimental data of different ODCs were compared and analyzed. It was shown, how different 
specific catalyst surface areas or electrolyte intrusions affect the performance of ODCs [48]. The 
validity of the results of Sudoh et al. and Wang and Koda remain debatable, since important 
aspects such as water and ion mass transport were not considered. 

The model of Wang and Koda was also revised and further developed by Pinnow et al. [32, 39]. 
Essential characteristics of the stationary pseudo two-dimensional TFFA model are Maxwell- 
Stefan diffusion in the gas and the liquid phase, a complete, homogeneous ORR as well as an 
oxygen solubility independent of the NaOH concentration of the liquid electrolyte. Especially 
the latter point is worth discussing, since Tromans determined a strong dependency of the oxygen 
solubility into the liquid electrolyte dependent on ion concentration [22]. In the model structure, 
a gas diffusion layer is considered at the interface to the porous electrode, but no liquid diffusion 
layer. Generally, it can be stated that the mass transport of water and hydroxide ions in the 
liquid electrolyte and their influence on the phase equilibrium seems to be underrepresented in 
the literature model by Pinnow et al. [32, 39]. The model predicts a fast oxygen mass transport 
in the gas phase and slow mass transport of the oxygen in the liquid phase which leads to a steep 
oxygen gradient within the thin-film. The limitation of the electric current is attributed to the 
depletion of oxygen in the liquid phase, which results mostly from the diffusion resistance of the 
oxygen cross the thin-film [32]. 

Although the dynamic one-dimensional three-phase model presented in this thesis was developed 
independently, the model structure, the considered processes as well as the determination of the 
material value and model parameters are based on the stationary pseudo-two-dimensional three- 
phase model by Pinnow et al. [32, 39]. 


3.2 Phase transitions 


At the gas-liquid interface within ODCs, oxygen dissolution and water evaporation take place [8, 
32]. As shown in different studies, the ion concentration in the liquid phase has a distinctive effect 
on oxygen solubility [22, 49, 50] and on water evaporation [5, 51, 52]. Botz et al. demonstrated 
that the NaOH concentration in ODCs can increase, depending on the current density, from 
10M at open circuit potential (OCP) conditions up to 20M during operation [34]. Therefore, 
the consideration of ion concentration in both processes at the gas-liquid interface appears to be 
essential. 

The solubility of oxygen in the electrolyte can generally be described by Henry's law: 


E Po, 
Ho, 


co (3.1) 
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Where co, is the concentration of the dissolved oxygen in the liquid phase, po, is the partial 
pressure of gaseous oxygen, and Ho, is the inverse Henry constant. A calculation of Ho, in 
NaOH electrolytes is proposed by Tromans [22]. Although the calculations are validated only 
up to 6 M, Zhang et al. were able to show that these calculations agree with experimental data 
up to 12M [27]. Since no data validated beyond 12 M were found despite an extensive literature 
search, the calculations have to be further extrapolated in this thesis. 
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Figure 3.2: a) NaOH concentration dependent inverse Henry constant at 80?C, calculated based on [22]; b) NaOH 
concentration dependent normalized water vapor pressure at 80 °C, calculated based on [51]. 


Fig. 3.2 a) shows the Henry constant for oxygen solubility at T = 80°C as a function of the 
NaOH concentration. In a 16 M electrolyte the oxygen solubility is 5 times lower than in a 10M 
one. Using Henry's law with a oxygen partial pressure of 1 atm, the concentration of dissolved 
oxygen drops from 26.7 x 107° to 5.04 x 10~3molm~? in the given interval. A comparison of 
further literature sources on oxygen solubility in the concentrated NaOH electrolytes is given by 
Pinnow [39]. 

The vapor pressure of water Pho in concentrated NaOH solutions is documented by Hirschberg 
[51]. The normalized water vapor pressure over the ion concentration at T — 80?C is shown in 
fig. 3.2 b). While for a 10M electrolyte the water partial pressure constitutes about 25 % of the 
gas phase pressure, the value drops to less than 10 46 for a 16 M electrolyte. 

In summary, both oxygen solubility and water vapor pressure decrease significantly with in- 
creasing NaOH concentration. Only small amounts of oxygen can be dissolved in the liquid 
electrolyte, but water partial pressure in the gas phase can reach 25 % of the total pressure. It is 
concluded that the partial pressure of water in the gas phase influences the availability of oxygen 
in the liquid phase, but the dissolved oxygen does not have a significant impact on the water 
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evaporation. However, the ion concentration dependencies of water evaporation and oxygen 
dissolution have to be considered in the model. 


3.3 Reaction kinetics in oxygen depolarized 
cathodes 


In literature studies, the reaction kinetic of ORR in a strong alkaline milieu has been modeled by 
its microscopic steps [24] or as lumped kinetic [29, 32]. 

A common way to calculate the reaction rate of ORR is the Butler- Volmer Equation (BVE) [46]. 
Detailed derivations are provided in the literature, e.g. by Bard and Faulker [53] or Hamann and 
Vielstich [4]. Herein, the BVE equation is adapted to the specific problem of the ODC, which is 
its applicability for multi-step reactions and simplifications due to large over-potentials. 

The current density jp for a general, reversible electrochemical reaction eq. (3.2), in which an 
oxidized species O reacts under the acceptance of n electrons to the reduced species R, can be 
described using Faraday’s law, where the faradaic current density jp is the sum of the oxidation 
Jox and the reduction jox current densities [53]. 


k 
O+ne ==R (3.2) 
red 
jp = Jox — Jed = OF (Pox — Fred) (3.3) 


Where Fis the Faraday constant and rox and req are the reaction rates. To describe the reaction 
rates of a multi-step reaction — as the ORR -, it has to be kept in mind, that the overall reac- 
tion rate of (electro-)chemical reactions is limited by its slowest, rate determining step (RDS). 
Thus, the kinetics of the RDS limits the total reaction rate: rox = rgps, ox and Fred = RDS, red- 
The forward and backward rates of the RDS as single reaction step can be resolved according 
to the Butler-Volmer equation, depending on the overpotential 7, the kinetic reaction rate con- 
stant Kored/ox and the activity apps red/ox,i Of the reactants involved in the RDS and v as their 
stoichiometric coefficient [53]: 


aFz 
TRDS,ox = Less; t ko,ox exp (= n) (3.4) 
1 — a)Fz 
l'RDS,red = Less; 3 Ko.red exp (- m) (3.5) 
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Current density j / A m? 
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Figure 3.3: Dependency of current density on overpotential close to the OCP, calculated by Bulter-Volmer and 
Tafel equation. Following data are used: [[a&ps.4; = [l4Rps oxi = b: n — z — l; T = 298.15 Ki Korea = koox = 
0.1 m? mols-!; q =0.5. 


Where & is the (anodic) symmetry factor, R the universal gas constant, T the temperature and 
z the transferred electrons in the RDS. Substituting the reaction rates into Faraday's law, the 
Butler-Volmer equation for a multi-step reaction is obtained as: 


. aFrz 
JF aN: El [aks yea. $ Ko,ox exp (s ) 


] — «)Fz 
—n- F] [anps.o; ` ko,red exp (- m) 


(3.6) 
With this, the current density of an electrochemical multi-step reaction can be described with 
respect to the RDS. It should be noted that this only applies if the RDS is an electrochemical 
and not a chemical reaction [53]. It is also possible that the RDS may differ depending on the 
electrode potential [32] or on the activity of the reactants. 

The current dependency due to the BVE on the overpotential is given in fig. 3.3, based on 
eqs. (3.3) - (3.6). The figure shows the total current density jr, its oxidative jo, and its reduc- 
tive component jreq depending on the overpotential. As illustrated, at potentials near the OCP 
—0.05 V < 7 x 0.05 V, both partial currents have a considerable influence on the total current. At 
high positive overpotentials y > 0.05 V, the total current density jr is governed by the oxidative 
current jox, while the reductive current jreg dominates at negative overpotentials yn < —0.05 V. 
These properties can be used to simplify the BVE in case of high negative or positive overpoten- 
tial. At high negative overpotentials n «& 0 V, in eq. (3.5), exp (- Eom) becomes a high 


value, while exp (ss ) of eq. (3.4) approximates zero. Thus, at high negative overpotential, 
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Jox can be neglected and the Butler- Volmer equation can be reduced to the so-called Tafel equa- 
tion, where j = — jred. For high positive overpotentials n << O0 V it works the other way round. 
However, this leads to the fact that the description of the reaction rate loses its validity near the 
OCP [53]. Finally, the current density of the ORR under the given conditions as an electrochem- 
ical multi-step reaction at higher negative overpotentials n < —0.05 V can be simplified to the 
Tafel equation, which takes the total number of electrons transferred and the reaction rate through 
the RDS into account: 


m) (3.7) 


j=—n-Fko- [Jaks ox. exp (- RT 


The overpotential used in the BVE and the Tafel equation is defined as the difference between 
the potential in the reaction zone E and the open circuit potential of the electrode Eo: 


n =E -— Ep (3.8) 


The open circuit potential results from the thermodynamics of the electrode and describes the cell 
potential at a net current flow of j =0Am'?. Ep can be calculated from the tabulated value of 
the standard electrode potential Eoo corrected by the activities a; ox reg of the oxidized or reduced 
species by the Nernst equation [4]: 


Tal 
zen ) (3.9) 


Ev = E È 1 ( 
0 = £00 n 
nF [vil 
IIei rea 
t 
n is the number of electrons transferred in the overall reaction. Regardless of the RDS, in a 
multi-step reaction the activities of the total reaction are used in the Nernst equation to determine 
the OCP. In contrast to the kinetics, not only the activities of the reactants but also the activities 
of the products influences the thermodynamics of the reaction [53]. 


3.4 Mass transport in oxygen depolarized cathodes 


In the gas phase of ODCs, water and oxygen are transported. The liquid phase is characterized 
by multi-component mass transport of water, hydroxide ions and the dissolved oxygen. The 
complexity is further increased by the phase transition, the electrochemical reactions and the 
high ion concentration of the liquid electrolyte. 

In an electrochemical system, a distinction between three different mechanisms of mass transport 
is made [53]: 
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1. Diffusion: Mass transport due to a gradient in chemical potential, in idealized systems 
also due to a gradient in concentration or partial pressure. 


2. Convection: Mass transport due to a gradient in density (natural convection) or as a forced 
convection due to pressure gradients. Additional convection can apply due to reactions or 
phase transitions in form of the so called Stefan-flow [54]. 


3. Migration: Mass transport of charged species due to a gradient in electrical potential. 


In a binary NaOH-electrolyte, due to the constraint of electroneutrality, the electrical potential 
or rather the gradient of electrical potential can be eliminated as driving force from the mass 
transport equation. Therefore the migration can be neglected [55] and is not considered further. 
For the diffusion process it must be distinguished between two basically different approaches: 
Fick's diffusion and Maxwell-Stefan diffusion. In contrast to Maxwell-Stefan diffusion, Fick's 
diffusion is based on the assumption of an idealized system [56]. Moreover, Fick's approach can 
only be transferred to multi-component diffusion in concentrated solutions to a limited extent 
[57, 58]. 

The bi-mixture of oxygen and water in the gas phase can be idealized in terms of mass transfer, 
thus the gradient in partial pressure can be used to calculate the diffusion flux. An additional 
convection may occur due to the oxygen dissolution and water evaporation. This superposing 
convection is considered as the Stefan-flow. The mass transport in the liquid electrolyte is char- 
acterized by high ion concentrations and thus a strictly non-ideal behavior, multi-component 
mass transport as well as phase transitions. Consequently, these properties are considered by the 
Max well-Stefan diffusion in combination with the Stefan-flow in the herein developed model. 
In the following, the Maxwell-Stefan equation for the specific system of ODCs is derived. Sub- 
sequently, the Stefan-flow is described on the example of Fick's diffusion. 


3.4.1 Maxwell-Stefan diffusion 


In this chapter, the dynamic Maxwell-Stefan relation for the specific system of the liquid phase 
within an ODC is derived beginning from the ordinary Maxwell-Stefan diffusion for the station- 
ary case. Herein, the diffusion is discussed exclusively; the species balance, including sinks and 
sources, are considered in chapters 4, 5 and 6. 

For a better understanding of the Maxwell-Stefan relation, the derivation of the Maxwell-Stefan 
diffusion for the stationary case is attached in appendix A and briefly summarized hereafter. 
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As illustrated in fig. 3.4, the basic idea of the Maxwell-Stefan diffusion is that the drag force of 
species i based on the gradient in its chemical potential is equal to the forces acting on species i 
due to its relative velocity to the surrounding species j (i Z j) [59, 60]. 


Z 
Driving force acting on species i = Force extended by species j on species i 


Figure 3.4: Schematic illustration of Maxwell-Stefan diffusion for the binary case. Partly based on [56]. 


The drag force di- RT - Cora] acting on species i due to its gradient in chemical potential is given 
by eq. (3.10) [60] (cf. A.11). 


dui 
di: RT - Ctotal = "de * Ctotal (3.10) 
z 
Where Ctota 18 the total concentration of all species, ai is the gradient in the chemical potential, 
x; is the mole fraction and dj is a coefficient for the forces acting on the molecules of species i. 
The force di- RT - Crotay Which is applied on the molecules of one species i due to their relative 
velocity to the other present species j is calculated by eq. (3.11) [59] (cf. A.10). 
"xix? (vi—-vj 
di- RT - Ctotal = y EM | j) -RT - Ctotal Q.11) 


j=l Lj 
¡Aj 


Here Dj ; is the binary Maxwell-Stefan diffusion coefficient of species i and j, and v; the average 
velocity of species i. 

The Max well-Stefan diffusion for the stationary case is obtained by equaling both forces given 
in eq. 3.10 and 3.11 and using the relation c; = xj: Ctota (cf. A.12): 


Cj du; 4 cjci-(vj—vi) 
RT d 


(3.12) 
j=l CrotalDi,j 
ii 
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The transformation of the stationary Maxwell-Stefan diffusion into the dynamic form begins with 
the calculation of the molar flow densities: The net flux N; of the species i can be calculated by 
[59]: 


N; = — Y Bij-cuu:d;j (3.13) 
j=l 
izj 


By substituting d; j by the right hand side of equation 3.11, we get: 


xx (vi — vj) 


n 
Ni = — Y Di; Crotal — (3.14) 
hn ij 
m 
which can be simplified to 
N; = p cer (Vi vj) (3.15) 
j=l Ctotal 
if 


Now the resulting flux N; of the species i is given as a function of the relative velocities (v; — 
vj) causing momentum exchange between different species and the frequency of the collisions, 
represented by cjc;. That equation gets clarified by noting that the Maxwell-Stefan relation 
only allows n — 1 independent fluxes for n interacting species [59]. The velocities v; and v; 
depend on the driving forces of the species. Since in the binary NaOH-electrolyte the migration 
can be neglected and convection only may occur as a Stefan-flow, which leads to a shift of the 
system and not to a change of the relative velocity (v; — vj) between the different species, an 
exclusive diffusion as driving forces for v; and v; can be assumed. For an exclusive diffusion in 
a concentrated electrolyte, the velocity v; can be calculated by eq. 3.16 [60] The calculation of 
vj is done analogously. 


. Bij di 
“RT dz 


(3.16) 


Please note, that eq. (3.16) and all following steps would change, if e.g. migration-based fluxes 
respective to the surrounding medium would occur. An overlaying convective flow, which moves 
as species with the same velocity and direction is still valid. From eq. (3.15) and (3.16), eq. 
(3.17) is obtained. 


: "Di; cic; [ duj du; 

Ni= en 3.17 
2 RT Ciotal | dz dz S 
Aj 
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For the diffusive mass transport, the change of concentration depending on the time equals the 
derivation of the flux N; over the location z: 
dc; dN; 


qe (3.18) 


Using eq. (3.17) and (3.18), the change of the concentration depending on time can be calculated. 


Since cj, cj, Wi, Lj and Crora are depending on the location, the total differential of (3.18) is 


needed, which is given in eq. (3.19). A constant diffusion coefficient 757 = -g =0 is assumed. 
de; Bi) du cic; u; cic; 
dt -L RT | oz  Ctpial 92? * Ctotal 
iAj 
" OM; de; Cj du; oci Cj 
Oz 02 Ctotal AZ AZ Ctotal (3.19) 


ou; dcj Ci OL; dcj Ci 


Oz Oz Cita OZ Oz Cota 
du; . CiCj Oca OL; . CiCj aoa ) 


oz 0% c, Oz 


2 
Ctotal total 


As discussed in chapter 4 and appendix B in the range of industrial relevant working conditions 
for the advanced chlor-alkali electrolysis, the total concentration Ctota) of the liquid phase can be 
handled as constant. With a constant Ctotal, eq. (3.19) can be simplified to: 


dc; = a Dij 9?p; ou; CiCj 
dt => RT (( ax gg ) 


Ctotal 
Aj 


| (7 pI Cj (3.20) 
Oz  Oz/ Oz Cota 


< a) Ci 
"Naz 927 Oz Ciotal 


Equation (3.20) describes the time-dependent change in concentration of species i, respecting 


interaction with all species j Z i and assuming a constant total concentration. The driving force 
is exclusively diffusion, based on gradients in chemical potential. 


A simplification of the equation (3.20) can be achieved, if an ideal behavior of the involved 
species is assumed. Using activity coefficients y; = 1, the concentration is obtained as the driving 
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force instead of the chemical potential. As explained in detail in appendix A, in this case equation 
(3.20) can be simplified to equation (3.21): 


dc; de e; c; e; cj 
Hi y ea TES (3.21) 
dt 2, E | oz Ctotal oz Ctotal 

¡Aj 


Eq. (3.21) is the Maxwell-Stefan diffusion which regards the interaction between the species 
i and j and the gradient in concentration as the driving force. All further before mentioned 
assumptions are still valid. 

With the aid of the equations (3.20) or (3.21) developed here, multi-component mass transport in 
the liquid phase can be described dynamically for the specific system of ODCs. 


3.4.2 Fick's diffusion and Stefan-flow 


In the gas phase of an ODC, oxygen and water are present, which can be assumed to behave as 
ideal gases. Thus the diffusion can be calculated by Fick's law. The aim of this section is to 
expand the pure Fick's diffusion by the overlaying convective Stefan-flow, which may occur in 
ODCs as a result of the water evaporation and oxygen dissolution. Herein, Fick's first law is 
taken as given. 

Dependent on the properties of a system, three different mechanisms of diffusion can occur [54]: 


1. The bilateral equimolar diffusion 
2. The bilateral non-equimolar diffusion 
3. The unilateral diffusion 


In the first case, the bilateral equimolar diffusion, the number of molecules moving in opposite 
directions are equal to each other. This mechanism can be illustrated by two closed vessels 
connected by a tube. The vessel 1 is mostly filled with the gas A, the vessel 2 is mostly filled 
with the gas B, no gradients in temperature or pressure are present (see fig. 3.5 a)) [54]. 


b) 


Gas B 


[a |n 


Liquid A 


Liquid A 


Figure 3.5: Schematic illustration of a) bilateral equimolar diffusion (based on [54]); b) bilateral non-equimolar diffu- 
sion; c) unilateral diffusion. 
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In an idealized system, it can be assumed that gas A diffuses through the connecting tube towards 
vessel 1 and gas B diffuses through the connecting tube towards vessel 2. The fluxes of both gases 
can be calculated by the mean of Fick's first law, with Na and Ng as the fluxes of the respective 
species, pa and pg as the partial pressures, z as the location and Zap as the Fick's diffusion 
coefficient: [54] 


_ Za dpa 


Ni = 3.22 
^ RT dz ( ) 
. Dap dp 
BEN SL na 2 
Np RT dz (3.23) 
Since no gradient in total pressure occurs, Pa + pp = Protal = const., this results in dea =- aps, B. 


Further taking Ya B = Zg.A into account, it becomes obvious that Na = —Ng and that the number 
of the transported molecules of both species are the same [54]. This is not the case for bilateral 
non-equimolar diffusion. 


The bilateral non-equimolar diffusion is particularly important when the diffusion process is 
associated with sinks and sources. Examples for sinks and sources can be (electro-)chemical re- 
actions or phase transitions. The mechanism is illustrated by the example of the phase transition. 
Here, a model system can be a beaker half filled with a liquid A and half with a gas B, whereby 
the liquid evaporates into the gas phase and the gas dissolves into the liquid (see fig. 3.5 b)). It 
can be expected that the flux of the evaporating liquid Ng and the flux of dissolving gas NA will 
be unequal in amount. The bilateral non-equimolar diffusion is based on the assumption, that the 
inequality of both diffusion fluxes is compensated by superimposed convective fluxes vs pa and 
vs - pp, the so-called Stefan-flow, with vs as the velocity of the convective mass transport. The 
fluxes of both species can not longer be described by a pure diffusive term as in eq. (3.22) and 
eq. (3.23), but also with a convective term: [54] 


: Da. dpa PA 
Na = > "E 3.24 
3-7 RT dz | RT pew 

ZABpB PB 


(3.25) 


Here again, at a constant pressure and temperature, the absolute values of the actual diffusion 
fluxes of both components are equal to each other: Day La = -Ip A Ya. In contrast, the 
convective Stefan-flow and thus the total fluxes are unequal, so that Na 4 —Ng. The velocity 


of the Stefan-flow v, has to be determined — according to the system — from the fluxes across 
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the interface or the stochiometry of the reaction. The change in partial pressure of the species 
dependent on the time can be determined by eq. (3.18) and is given for the general case with 
Vs — const., as [54]: 


9? pj 
92 


dpi 
ot 


=D j (3.26) 
Since the Stefan-flow is a shift of the whole system (same velocity and flow direction for all 
components), the approach of the Stefan-flow can also be combined with the Maxwell-Stefan 
diffusion. 


The last case, the unilateral diffusion, represents a special case of the bilateral non-equimolar dif- 
fusion. Examples are given by a two-gas mixture in a chamber sealed with a selective membrane, 
or by a two phase system with either evaporation or dissolution (see fig. 3.5 c)) [54]. However, 
this case is not considered further here, since this mechanism is not expected in ODCs. 


3.5 Dynamic analysis of oxygen depolarized 
cathodes 


The performance of ODCs is determined by a complex interplay of its reaction kinetics, mass 
transport in gas phase and mass transport in liquid phase which is difficult to access in steady 
state simulation. One way to separate the influences of the different processes on electrode per- 
formance from each other is the utilization of dynamic methods [35]. The principle of dynamic 
methods is based on the excitation of the system by a time-dependent change of at least one 
process parameter (e.g. potential, current, temperature). From the time-dependent system re- 
sponse, inter alia, the characteristic time constants 7 of the excited processes can be determined. 
In fig. 3.6 different processes which may take place within ODCs and also suitable dynamic 
electrochemical methods to study them are assigned to their characteristic time constants. 
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Figure 3.6: Time constants of processes in ODCs and suitable methods to detect them, Az is the characteristic diffusion 
length, partly based on [13, 35]. 


The time constants of the processes within ODCs, can range from milliseconds to several min- 
utes. In ODCs diffusion occurs in gas and in liquid phase. For both cases, the time constant is 
significantly increasing with the diffusion length Az. But for the same diffusion length, the time 
constant of gas and liquid diffusion are magnitudes apart. This can be explained by the equation 
of the diffusion time constant, which is derived from Fick's second law [13]: 

Az? 


fmt ^: > (3.27) 


Where Tmt is the time constant of the diffusion process, D is the diffusion coefficient and Az is 
the mass transport length. Here it can be seen that the time constant increases quadratically with 
the diffusion length and decreases proportional with the diffusion coefficient. Typical values 
for the gas phase diffusion coefficients are about 1 x 1076 to 1 x 107^ m*s~! and for the liquid 
diffusion about 1 x 1071? to 1 x 1078 m? s^! , which explains the strongly different time constants 
of diffusion in the two phases [61]. 

Another characteristic time constant is received due to the charge and discharge of the double 
layer. This time constant can differ depending on further electrode processes. For example has 
the time constant of double layer charge parallel to faradaic reactions to be calculated differently 
than the time constant of double layer charge of a blocked electrode surface without faradaic 
reaction. For a parallel reaction — as expected in ODCs - the following equation is given [62]: 


Tar = Car: Ret (3.28) 
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With the characteristic time constant Tq, the double layer capacitance Ca; and the charge transfer 
resistance Ret. The charge transfer resistance is a complex parameter that is determined by the 
reaction kinetics of the electrode. The double layer capacitance can be effected by many factors, 
such as catalyst loading, ion concentration in the electrolyte or applied potential [4]. 

Additional processes which may be dynamic with characteristic time constants are e.g. the heat 
conduction for non-isothermal systems [13], or degradation processes. Óhl et al. showed that 
degradation of the silver catalyst in ODCs progress in the range of weeks [63] and is therefore 
negligible for this work. 


Fig. 3.6 also gives an overview of suitable methods to analyze the processes depending on 
their time constants. Herein, dynamic methods are focused on the electrochemical impedance 
spectroscopy (EIS) and the potential steps. An overview about further dynamic electrochemical 
methods and their field of application are provided by Bard and Faulker [53] or Hamann and 
Vielstich [4]. While for faster time constants EIS is a suitable, high resolution tool, potential steps 
are feasible to detect slower time constants, especially for diffusion processes and cover a wider 
range of time constants. Electrochemical impedance spectroscopy has been used successfully 
e.g. for identification of reaction kinetics in DMFC [35], or to analyze catalyst poisoning in 
proton exchange membrane fuel cell (PEMFC) cathodes [64]. By applying potential steps on 
DMFC, the influences of kinetics and transport have been examined [37]. 

The potential step — also called (multi-)step chronoamperometry — is an instationary method, in 
which the electrode potential is changed abruptly at a defined time as input signal. The time- 
dependent current density is obtained as output signal [53]. Input, output and the triggered 
processes within the ODC are illustrated in fig. 3.7. 
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electrode 


Figure 3.7: Schematic illustration of changes in states during a potential step: a) potential dependent on time as input 
signal; b) double layer charge / discharge (illustration for ion concentrations > 1 M) as fast process; c) mass transport due 
to electrochemical surface reactions as slow process and d) current density dependent on time as output signal. Partly 
based on [4, 53]. 


The change in potential (fig. 3.7 a)) triggers the double layer charge (fig. 3.7 b)), the faradaic 
reactions and thus the mass transport (fig. 3.7 c)). Other processes may also be stimulated, but 
will not be considered further here. The current density (fig. 3.7 d)) is characterized by a peak 
immediately after the step, followed by a slower relaxation. The peak is assigned to the fast 
process of the double layer charge. The subsequent relaxation is the result of a slower diffusion 
process. The time constants of the double layer charge / discharge Tq, and the diffusion Tmt can 
be inferred directly from the time-dependent response of the current. 

As also shown in fig 3.6, each characteristic time constant is associated with a characteristic fre- 
quency fchar. The time-behavior of the individual processes can be transferred into the frequency 
domains by: 


1 
char = Int (3.29) 
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The quasi-stationary EIS makes use of this relationship. Herein, a short overview about EIS is 
given, deeper insight into the theory of EIS is presented by Orazem and Tribolet [62], the effects 
of different processes on the impedance spectra are evaluated by McDonald [65]. The general 
principle of EIS is sketched in fig. 3.8. The electrochemical system is stimulated by a sinu- 
soidal potential E with the frequency f and a small amplitude, the sinusoidal current density j 
is considered as the output signal [4]. At a certain frequency, the processes with the associated 
characteristic frequency are triggered, cf. fig. 3.8 a) - d). To excite different processes according 
to their specific frequency the procedure is repeated for several orders of magnitudes of frequen- 
cies. 

For each frequency f, the impedance magnitude |Z(f)| and the phase shift p(f) is determined. 
As illustrated in 3.8 e), the phase shift p(f) is obtained from the time shift Ar(f) between the 
sinusoidal curve of the potential and the sinusoidal curve of the current. 


Qf) =At(f)-2af (3.30) 


The impedance magnitude |Z(f)| is calculated from the amplitude of the sinusoidal potential E 
and the sinusoidal amplitude of the current J(f) [62]. 
E 
IZ(£) ^ 5 (3.31) 
JCf) 
Based on |Z(f)| and 9 (f), the impedance can be divided into its real Zp (f) and imaginary parts 
Z(f): 


Zr(f) = cos(p(f))  IZCf)] (3.32) 


Z(f) = sin(e(f)) ZN (3.33) 


One way of presenting the resulting spectra is the Nyquist plot, in which usually the negative 
imaginary part of the impedance for all frequencies —Z is plotted over the real part Zr. In 
fig. 3.8 f) a Nyquist plot for a electrochemical system with both processes, double layer charge 
/ discharge parallel to a faradaic reaction as well as mass transport, is illustrated. Please note 
that the contribution of a process to the electrode dynamics can be negligibly small, in such 
a case the process will not be visible in the Nyquist plot. In the Nyquist plot, each process 
is usually represented by a semi-circle. It should be mentioned that other shapes, such as a 
Warburg line with a 45? slope, are possible; however, since these are not expected in an ODC, 
here they will not be considered any further. Two essential characteristics of the semi-circles 
are their diameters and their characteristic frequencies. The characteristic frequency fchar 18 the 
frequency at which the imaginary value of the semi-circle reaches its maximum. The diameter of 
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Figure 3.8: Schematic illustration of electrochemical impedance spectroscopy analysis and the related dynamic state 
changes: a) sinusoidal potential as input signal and current as output signal in high frequency range; b) double layer 
charge / discharge (illustration for ion concentrations > 1M) as fast process; c) sinusoidal potential as input signal 
and current as output signal in low frequency range; d) mass transport due to electrochemical surface reactions as slow 
process; e) time shift Ar and amplitude of sinusoidal potential Ê and current J; f) Nyquist diagram with an charge transfer 
(Re) and an mass transport (Rmt) semi-circle. Partly based on [4, 53]. 
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the semi-circle reflects the resistance R of the process. However, if semi-circles of two processes 
are overlapping, the analysis is getting more complex. From the distance between the real part 
of the impedance at the highest frequency Zr (fmax) and the ordinate, the ohmic resistance of the 
system can be determined [62]. 


3.6 Concluding remarks 


In this chapter, fundamentals for developing a dynamic macroscopic ODC model were discussed 
and a selection of dynamic analysis methods suitable for ODCs was presented. In conclusion, 
following points must be considered for developing and analyzing a dynamic macroscopic model 
of an ODC: 


In order to simulate the measured high current densities of ODCs, a TFFA approach seems 
to be suitable. 


The Henry-constant for oxygen solution and the vapor pressure of water for water evapo- 
ration have be set in dependency of the ion and water concentration in the liquid phase. 


The ORR can be described using lumped kinetics and the Tafel equation. 


Due to the multi component diffusion and the high ion concentration in the electrolyte, 
liquid mass transport have to be calculated by Maxwell-Stefan diffusion. 


Migration can be neglected. 


Due to the phase transitions, a Stefan-flow has to be considered in the liquid and the gas 
phase. 


As already known from previous studies, the reaction rate is strongly dependent on local 
concentration gradients [32]. Thus, the model has to be discretized. 


The interaction of the different processes can be analyzed by their characteristic time con- 
stants. While the potential steps cover a high range of time constants, the EIS is high- 
resoluting tool especially for fast processes. 


Based on the outcome of this chapter, in chapter 4 a basic model will be developed which depicts 
all relevant processes. In consideration of (intermediate) results, the model will be extended 
systematically in chapters 5 and 6. Each development stage of the model is customized to answer 
specific research questions: 
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Chapter 4 : Basic model 
By means of the basis model, the time constants of the processes of an ODC will be deter- 
mined. It will be shown that the current density is limited by the availability of dissolved 
oxygen in the liquid electrolyte and that this is closely linked to ion and water mass trans- 
port. 


Chapter 5 : Expanded liquid transport model 
The model focuses further on mass transport in the liquid phase. The significance of water 
and ion mass transport will be evaluated in more detail and how these govern the dynamics 
and performance of ODCs. 


Chapter 6 : Inhomogeneous model 
The inhomogeneous model enables to consider the distribution of the electrolyte and thus 
the distribution of the gas-liquid interface within the ODC. It will be shown that only 
certain areas in the ODC are electrochemically active, although the three-phase boundary 
extends inhomogeneously over the entire ODC thickness. 


The characteristics, implemented processes and assumptions of each model are given in the cor- 
responding chapters. An overview about the differences between the models is given in table 3.1. 


Table 3.1: List of properties varying in the different model expansion stages. 


Property / Process Basic model Expanded liquid Inhomogeneous 
mass transport model 
model 
Activity of O5 ideal non-ideal non-ideal 
Activity of NaOH not included non-ideal non-ideal 
Electrolyte distribution homogeneous homogeneous inhomogeneous 


Reaction kinetics 


Gradient for diffusion in 


liquid phase 


l-step kinetics based 
on gross reaction 


concentration 


l-step kinetics based 
On gross reaction 


concentration 


l-step kinetics based 
on rate determining 
step 

chemical potential 
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4 Processes and Their Limitations in 
Oxygen Depolarized Cathodes’ 


In this chapter, the basic dynamic three-phase model of ODCs is introduced and parameterized 
by using literature data, in order to analyze the processes and limitations within ODCs. 


4.1 Introduction 


To explore the performance limiting processes of ODCs and to lay the foundation for a sys- 
tematic improvement, in this chapter a dynamic one-dimensional three-phase model of porous 
ODCs for the ORR is introduced. To the best of my knowledge, it is the first dynamic model of 
an ODC and the first ODC model in general, which takes the phase equilibrium dependent on 
the water and ion concentration in the liquid electrolyte consequently into consideration. These 
attributes are e.g. beneficial for estimating the location of the gas-liquid interface. In contrast to 
the stationary models, e.g. the state-at-the-art thin-film flooded agglomerate model by Pinnow et 
al. [32], time constants can be simulated for a deeper analysis of the limiting factors and their 
interaction. 

In this chapter, polarization curves from literature are used to parameterize the model. Next, 
dynamic simulations are evaluated to determine the relevant time constants and the limiting pro- 
cesses of the system. Finally, the evaluation is completed with a sensitivity study to show the 
impact of model parameters on the overall performance of the ODC. Through this procedure, I 
will show that the mass transport of water hydroxide ions in the liquid electrolyte is determined 
as the overall restricting factor. 


5 Parts of this chapter have been published in Röhe et al., ChemSusChem, 2019 [1]. 
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4.2 Basic model: Model concept and assumptions 


Under operating conditions, the ODC is impinged by oxygen from the one side and by the liquid 
electrolyte from the other side (cf. fig 4.1)). Due to the combination of the hydrophilic silver 
catalyst and the hydrophobic PTFE, only a part of the ODC gets flooded so that three different 
segments within the ODC are expected: The section facing the oxygen inlet contains a pure gas 
phase, while the section facing the electrolyte is electrolyte flooded. In between the interface 
of both phases, the so called three phase area, is located. Processes occurring directly at the 
interface are the evaporation of water and the dissolution of oxygen into the electrolyte. Due to 
the poor solubility of oxygen into the electrolyte it is estimated that the ORR (eq. (1.1)) takes 
place in the first few nanometers of the electrolyte [32, 34]. In addition to the transport of the 
oxygen into the reaction zone, the provision of water and the removal of hydroxide ions are mass 
transport processes. 


To transfer the complex behavior of ODCs to a suitable dynamic model, the following assump- 
tions are made: 


1. No pressure drop occurs within the electrode. 
2. The temperature in the electrode is constant and homogeneous. 


3. Complete ORR (eq. (1.1)) is assumed. The possible incomplete ORR with hydrogen 
peroxide as final product is not considered, due the high catalytic activity of silver for 
reducing hydrogen peroxide [31]. The reaction is modeled as lumped one-step kinetics. 
So it depends on the availability of all reactants. 


4. Due to the high overpotential, the anodic current is negligible, hence the kinetic can be 
described with the Tafel-equation. 


5. The location of the gas-liquid interface is stationary over time. 
6. Within the electrolyte electroneutrality is assumed. 


7. Since the reaction zone in the electrode is very thin and the catalyst as well as the elec- 
trolyte are highly conductive, no internal voltage distribution is considered. 


8. The high concentration of NaOH results in a low viscosity of electrolyte, thus a liquid 
diffusion layer between the ODC and the electrolyte bulk has to be considered, even if the 
electrolyte bulk itself is ideally stirred. 
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9. A gaseous diffusion layer is not relevant. As will be shown later, the diffusion length in the 
gas phase does not contribute significantly to the total mass transport resistance and would 
not affect the ODC performance. 


10. In all liquid sections, the multi component diffusion is described by Maxwell-Stefan dif- 
fusion. Because the concentration of hydroxide ions and water is each about 10° times 
higher than the oxygen concentration, the influence of the oxygen diffusion on the water 
and hydroxide diffusion is negligible. 


11. The change of the electrolyte concentration due to the ORR causes a change of electrolyte 
density (cf. fig. B.5). With the further assumption of locally stationary phases, an overlay- 
ing convective flux is included, to consider the change of electrolyte volume. To calculate 
this flux, the simplification Crotal = Ctotal(z = 0,1 = 0) is used (cf. fig. B.1). The error is less 
than 0.8 %, a detailed error estimation is given in the appendix. Further, convection caused 
by the water evaporation and the oxygen dissolution is implemented in both phases. 


12. The concentration of Na* is calculated from the equation for electroneutrality. To consider 
its influence on the mass transport, the binary diffusion coefficient B sun of NaOH 
vs. H50 within the liquid electrolyte is calculated with the concentrated solution theory, 
which takes OH” as well as Na* into account [60]. 


13. In the ODC, the gradients in through-plane direction (z-direction) are decisive. A homo- 
geneous distribution of the phases and their components in x- and y-direction in-plane is 
assumed. 


The one-dimensional three-phase model is divided into four spatially resolved sections. As 
shown in figure 4.1, three of them are located within the porous ODC: A pure gas phase, a 
thin-film and a flooded porous structure, named flooded agglomerates in the following. The thin- 
film represents the interface of liquid and gas phase within the electrode. The fourth section is a 
liquid boundary layer beyond the electrode. For the thin-film, the flooded agglomerates and the 
liquid diffusion layer, hereafter the collective term liquid phase is used. 

The gas phase contains exclusively gases, namely oxygen and water vapor. While the oxygen 
is provided by the gas bulk phase on the left, the water evaporates from the liquid phase at the 
right hand side. Additionally oxygen is dissolved in the thin-film and transported into the flooded 
agglomerates where the ORR occurs. The thin-film itself is an exclusive transport layer, where 
no reaction takes place. The upstream liquid boundary layer is included, since a significant diffu- 
sion resistance is expected [34]. Since a homogeneous distribution of the gas-liquid interface and 
of the concentration in x- and y-direction is assumed, the specific surface area of the interface 
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Figure 4.1: Schematic view of the model including reaction and mass transport. 
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Figure 4.2: Mass transport processes at the gas-liquid interface. 


Sphase,interface 18 defined, which represents the phase interface area per geometric cross sectional 
area of the electrode: 


Aphase.interface (4.1) 


Sphase, interface = pw 
This parameter enables the one-dimensional model approach to meet the characteristic of a large 
gas-liquid-interface within the ODC, as displayed in figure 4.2. Since a steep oxygen gradient 
within a few nanometers distance from the gas-liquid interface is assumed [34], the meaningful 
area for dissolution and subsequent mass transport of oxygen is defined to be the specific interface 
area (cf. fig. 4.2). The thin-film represents the liquid side of the gas-liquid interface, hence the 
interface area is also relevant for the transport of water and hydroxide ions through the thin-film. 
In contrast, the area for transport in the gas phase and for the transport of water and the hydroxide 
ions in the flooded agglomerates is given by the geometric cross section. 
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4.3 Model equations 


4.3.1 Electrochemical reaction kinetics, double layer 
capacitance and internal resistance 


The Faradaic current density jp is modeled by using Farady's law (4.2) and the Tafel-equation 
(4.3). All given current densities are referenced to the geometric cross section: 


jp =—F-4-r (4.2) 


Where the reaction rate r is calculated by 


(1 — a) Be Arc (4.3) 


= ko: de OX 
r 0 ao, do »( R.T 


Where ko is the kinetic reaction rate constant, a; are the activities of the species i, Nreaction 18 the 
overpotential, & is the anodic symmetry factor and T is the temperature. 
The concentration and temperature dependent activity ag,o is based on empirical studies of Balej 
[66], while for the oxygen activity ao, an approximately ideal behavior is assumed: 

O, 


Ci 
ao, = fo, on (4.4) 
Co, 


Hence ao, is calculated with the fugacity coefficient fo, = 1, the standard concentration da, = 
I1 kmol m^? and the concentration of the dissolved oxygen in the electrolyte CO): 

The overpotential of the ORR Nreaction 18 defined as the difference between the open cell voltage 
Eo and the potential in the reaction zone Ein: 


Nreaction = Eint — Eo (4.5) 


Ohmic potential losses can be calculated according to Ohm's law as the difference between the 
external applied and the potential in the reaction zone: 


Eext — Eint = Rspecific i J (4.6) 


Here Ein: is the potential in the reaction zone, which equals the external applied potential Eext 
under consideration of the specific Ohmic resistance Rspecific- The current density j is the sum of 
the Faradaic and the capacitive current density j = Jr + ja. 
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Assuming a constant double layer capacitance Cg, the charge balance in the reaction zone is 
given, as: 
dE 


Ca = j-F-4r (4.7) 


4.3.2 Mass transport in the gas phase and processes at the 
gas-liquid interface 


For modelling the mass transport in the gas phase, the evaporation of water and the dissolution 
of oxygen at the gas-liquid interface have to be taken into account as major effects. Since the 
absolute values of these two fluxes are typically unequal, in addition to diffusion an overlaying 
convection has to be taken into account. This process is described by the bilateral non-equimolar 
diffusion, characterized by Fick’s diffusion with an overlaying convective Stefan-flow [54]. The 
resulting species balance for i € {H,O,O,} in the gas phase is: 


dpi 9°p; gas,eff dp; 
ae © ah er ds 


With p; as the partial pressure of the species i, t as the time, € as the porosity, z as the location 
in Z-direction, CUBE as the effective binary diffusion coefficient and vgas as the velocitiy of the 
convective flow. ges is calculated from the free diffusion coefficient DES, porosity € and 
tortuosity c: 


€ 
DI = Oe (4.9) 


The tortuosity is approximated with the Bruggemann relation x = 8705 [67]. 

For calculating the Stefan-flow, constant pressure Pxotal is assumed. The only sinks and sources 
for both components are dissolution and evaporation at the gas-liquid interface z = Zgas and the 
further inclusion of the ideal gas law results in the same convective velocity over the whole gas 
phase: 


dvgas 


de =0 (4.10) 


Thus the velocity at every point in the gas phase exclusively depends on the mass transfer over 
the gas-liquid interface: 


‘convection __ xyevaporation ‘dissolution 
total = Ngo — No, (4.11) 
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Figure 4.3: Illustration of model on water evaporation and oxygen dissolution immediately at the gas-liquid interface. 
While the oxygen is solved from the gas phase into the liquid electrolyte (top), the water evaporates from the liquid 
electrolyte into the gas phase (bottom). 


convection 


The resulting convective mass flow Ñ OI 


can also be described by the velocity of the flow 
in gaseous phase: 


comets — += "is (4.12) 


The evaporation of water is described by Raoult's law, where the partial pressure of water pj; o 
directly at the gas side of the interface equals the vapor pressure Pho: which is a function of the 
concentration of sodium hydroxide cy,oy and of water cj, o immediately at the interface. The 
vapor pressure dependence Pro (CNaOH CELO» T) is taken from by Hirschberg. [51] cx,oy and 
Cjj, o are calculated by considering the mass conservation during the evaporation process with the 
use of a steady state mass balance (see fig. 4.3): 


“TF,out ^ x,evaporation _ a in 
MN MD (4.13) 


Here, aoe is the mass flow of evaporating water, which equals the sink of water in the 


159"! as well as the source of water in the gas phase N55". By solving the given 


liquid phase Na 
equations, a dns implementation of Raoult's law can be obtained. " detailed description of 
this procedure can be found in Schróder et al. [12]. For technical implementation, this approach 
is used to interpolate the concentrations and in particular, the partial pressures from the center 
of adjacent gas volume element and thin-film volume element at the gas-liquid interface to the 
common border between these volume elements (cf. fig 4.3). 


A similar method is implemented for the dissolution of oxygen in the electrolyte, where the flux 


leaving the gas phase No, Gas.out ¡s equal to the flux entering the liquid electrolyte No, T 
NG out _ Ne nee in (4.14) 
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Here Henry's law comes into effect to determine the oxygen concentration Co, directly at the gas- 
liquid interface depending on the oxygen partial pressure po, where the inverse Henry constant 
for oxygen dissolution Ho, is a function of cy, 9, and Cj, o as well [22]: 


* 
Po, 


CO; = Fy TER F 0X (4.15) 
° Ho, (Chaon: Cno) 


The concentration dependent functions of the water vapor pressure Puo (CNaOH; CH,0) (fig. B.3) 
and the Henry constant for oxygen dissolution Ho,(cNaoH; CH;0) (fig. B.2) are shown in the 
appendix B. 


4.3.3 Mass transport in the thin-film 


The mass transport in the thin-film is characterized by a diffusion process with an overlaying 
convection of all species due to the evaporation. This results from the assumption of a non- 
moving interface and a constant thin-film volume. Since there are three species i in the thin-film, 
i € (OH, H50,0,], for the diffusion process the Maxwell-Stefan diffusion with the concentra- 
tion gradient as driving force is taken into account. 


dc; = y EE Cj 9?c; Ci ) pia de » (4 16) 
a, 0 j ^ 3- ' "liq.evaporation E 
ot j=l oz Ctotal oz Ctotal ae 


i 
LT gg 
if 
Again, the oxygen concentration is very low, thus no reasonable effect from the oxygen diffusion 


to the diffusion of the water and hydroxide ions is expected, which leads to j € [OH ,H50]. As 
in the gas phase, the liquid evaporation velocity viiq.evaporation 18 constant throughout the thin-film. 


dv liq,evaporation 


Ps =0 (4.17) 


Thus, the velocity of the convective flow is directly proportional to the evaporation rate. 


‘gas, in 
Niño 


(4.18) 


Vliq,evaporation = E: 
HO 


4.3.4 Mass transport in the flooded agglomerates 


The mass transport within the flooded agglomerates is also described by the Maxwell-Stefan 
diffusion with overlaying convection vjig. The effective diffusion coefficient takes the porous 
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structure into account and is calculated by equation (4.9). Moreover the sinks and sources due to 
the ORR have to be modeled. This yields the following species balance: 


dc; z 9?cj Cj 9?c; Cj liq, eff dc; 

or d L | 92 Ctotal Az um) "RE T Oz Maps ind 
14) 

With i € (OH, H20, 07}, j € (OH,H20). 

In addition to the above described convection caused by the evaporation, there is a further effect 

due to the reaction. The ORR causes a local change of water and hydroxide concentration, which 

consequently leads to a change in liquid density (cf. fig. B.5). Although the molecules of water 

and hydroxide per volume change depending on the reaction, the total concentration Crotal stays 

almost totally constant (cf. fig. B.1). Consequently, the additional Stefan-flow based on the 

density change between the reaction area and the electrolyte bulk phase can be described by the 

removal of the additional produced ions: 


a riz 
"isse (2) = (SI = v) f° zE az (4.20) 
ota. 
where 
dviiq react 
DUAE cp 4.21 
oe X (421) 


Here, r is the reaction rate, VORR the stoichiometric factor of water and VORR stoichiometric 
factor of hydroxide ions in the ORR. 

Finally the local total velocity vii for each location z can be calculated by the velocity of the 
evaporation Vjig,evaporation and the location-dependent velocity of reaction Vjig, react: 


Vliq (z) = Vliq,evaporation “+ Vliq, react (z) (4.22) 


4.3.5 Mass transport in the liquid diffusion layer 


The mass transport in the liquid diffusion layer is described by the same equations as the mass 
transport within the flooded agglomerates. In this layer, free diffusion occurs, so the free binary 
diffusion coefficient is used: 


dc; n e; Cj 9?c; Ci li dci 
ae L Fo) abe td 2 Sy) 4.23 
3r ( ) ij + By "d (4.23) 


j=l 022 Ctotal de Ctotal 
Aj 
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With i € {OH , H20, 03}, je {OH ,H,O}. 
Since there are no sinks and sources in the liquid boundary layer, the convective velocity vıaı is 
equal to that at z = z (cf. fig 4.1) of the flooded agglomerates: 


Vidi = Viig(z = zi) (4.24) 
, Where 
vigi 
=0 4.25 
de (4.25) 


4.3.6 Boundary conditions 


Due to the high flux of oxygen in the gas bulk phase, the following boundary conditions are 
defined: 


PO; (z = 0) = Ptotal (4.26) 


pujo(z = 0) =0 (4.27) 
Since the electrolyte bulk phase is assumed to be ideally stirred, oxygen is expected to be zero. 
co, (z= %) = 0 (4.28) 


The amount of water in the electrolyte bulk phase can be calculated as a function of sodium 
hydroxide concentration [68], which is constant. 


CNaOH (Z = Zt) = Cnaon(Z = 2,1 = 0) (4.29) 


cuyo lz = zx) = f(Cuaon(z = 24,1 = 0)) (4.30) 


4.4 Parameterization 


The model was implemented in MATLAB. For parameter identification, data by Pinnow et al. 
was used. The examined ODC was manufactured by a spray method with a content of 97 wt— 96 
silver and wt—% PTFE. The experiments were performed at 80°C and with a NaOH molality 
of muaon = 11.25 mol kg! [32], which under the given conditions is a NaOH concentration of 
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€NaoH(z = zi) = 10.1 x 10? mol m ? [51, 68]. All further constant model parameters, set to their 
respective values from the experiment, are given in table 4.2, and concentration dependent pa- 
rameters are listed in table 4.2. The literature data were measured from j = 0 to j = 4kA m~? 
and simulated above j = 4kA m ? by the same authors [32, 39]. 

Identified model parameter are the kinetic reaction rate constant ko, the specific interface area 
Sphase,interface, the specific resistance Rspecific and the length of flooded agglomerates zu. The 
manually adjusted parameters which best reproduce the polarization curve documented by Pin- 
now et al. are given in table 4.3. The corresponding curve in figure 4.4 shows, that the data 
matches very well, thus the model can reproduce the steady state polarization curve. 

The identified specific gas-liquid interface area Sphase,interface 18 about 2.5 times higher than the 
interface modeled by Pinnow et al. [32]. This deviation may originate from the here applied strict 
consideration of the ion concentration dependent equilibrium at the gas-liquid interface. The ion 
concentration dependency of the phase equilibrium also enables an estimation of the location 
of the gas-liquid interface: The average value of the distributed gas-liquid interface is located 
11.5 um from the end of the electrode facing the liquid side, so that only a small part of the ODC 
is flooded with the liquid electrolyte (fig. 4.2). Since the diameter of the pores of comparable 
ODC is typically less than one micro meter [69] and the PTFE has strong hydrophobic attributes, 
this value seems to be valid. Present theories about the shape of the gas-liquid interface are e.g. 
cylindrical electrolyte flooded agglomerates which are surrounded by the gas phase [32] or mi- 
croscopic gas channels reaching into the flooded electrode part [70]. When assuming a regular 
cylindrical shape of the catalyst agglomerates in the electrode which are oriented perpendicular 
to the electrode surface [39], agglomerate diameters of dag = 0.62 um are obtained. This com- 
parably small diameter may indicate microscopic gas channels into the flooded agglomerates, an 
non plane surface following the catalyst particles or a moving interface. 

As discussed in chapter 4.5.3, all identified parameters have a characteristic influence on the po- 
larization curve, which is indicated by different sensitivities of the overpotential to the parameters 
at kinetic, ohmic and mass transport governed operation conditions. Nevertheless, the identified 
parameter set may not be unique. Reliable parameterization would require additional, especially 
dynamic measurements, which enable a strict separation of slow and fast processes. However, 
based on the comparison to literature values, I assess the parameter set to be reasonable. 
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Figure 4.4: Simulated polarization curve in comparison to literature data [32]. 


4.5 Results and discussion 


4.5.1 Steady state analysis 


To develop a deeper understanding of the processes within the ODC, the states underlying the 
simulated polarization curve (cf. fig. 4.4) are analyzed in the following section in detail. The 
electrode potential is referred to the open cell potential Eo of the simulated process conditions: 
E = Eext — Eg. As shown in figure 4.5 a), the oxygen concentration depletes strongly within the 
electrode close to the gas-liquid interface. The drop gets steeper with increasing negative current 
densities. Even at medium current densities of about 5kA m "?, dissolved oxygen can be found 
only in the first 200 nm of the flooded electrode section. Although only a small fraction of the 
catalyst is involved in the reaction, the kinetic losses are comparably low, which indicates a high 
catalytic activity of the silver for the ORR in alkaline media. At higher current densities, only the 
region directly behind the interface stays electrochemically active, which supports the findings 
of [32, 34]. Here, two major effects are identified, leading to the transport limitation of oxy- 
gen: Firstly, for increasing current, more oxygen is consumed by the ORR, which is indicated 
by the steeper gradient of oxygen concentration in the thin-film and the flooded agglomerates. 
Secondly, the higher reaction rate at higher currents causes an enrichment of ion concentration of 
OH as the reaction product of the ORR and an accompanying decline in water concentration as 
an educt of the ORR (cf. figure 4.5 b). Since the Henry constant is a function of ion concentration 
(cf. fig. B.2), the solubility of oxygen decreases at higher current densities and limits the amount 
of dissolved oxygen in the liquid electrolyte. This effect can be recognized by the decreasing 
oxygen concentration at the gas-liquid interface. The high sensitivity of oxygen solubility to the 
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Table 4.1: List of constant parameter and operation conditions. 


Name Symbol Unit Value 
Electrode geometry 
Electrode thickness [32] Zelectrode m 300 x 1076 
Thin-film thickness [32] Zee m 60 x 10-? 
Liquid diffusion layer thickness^ Zldl m 50 x 10-6 
Porosity [32] € — 0.4 
Tortuosity [67] T — 1.58 
Operation conditions 
Pressure in gas chamber? [32] P =Po, Pa 1x 10° 
NaOH concentration in electrolyte bulk? cnaou (z = 21) molm” 10.1 x 10? 
[32] 
H,O concentration in electrolyte bulk [51] cH,0(2 =) molm” 49.7 x 10? 
Temperature [32]? T K 353.15 
Electrochemical and kinetic data 
Open circuit potential [71] Eo V 0.216 
Symmetry factor“ [24] a — 0.15 
Double layer capacitance [72] Ca Em? 50 
Diffusion coefficients 
Binary gas diffusion coefficient [73] 2E o m?s=! 2.98 x 10? 
Binary liquid diffusion coefficient“ [73,74] Po 4.0 ms! 2.12 x 107? 
Binary liquid diffusion coefficient? [73,74] ^ Do Naon m?s”! 3.16 x 107° 
Binary liquid diffusion coefficient" [60] DE sn m?s”! 1.19 x 107? 


@ Own estimation. ^ Experimental operation conditions from [32]. ^ Measured with 6.5 M NaOH. 
4 Calculated with modified Wilke and Chang equation, © = 3.9 and a = 0.5 chosen, to fit the in [74] 
given experimental data best. ^ No NaOH data available, modeled for KOH with NaOH viscosity from 


[68], as seen in [32]. 


Table 4.2: Electrolyte concentration depended variables. 


Name Symbol Unit 
Inverse Henry constant^ [22] Ho, Pam? mol”! 
Activity of H,O? [75] aÑo = 
Density of electrolyte [68] p ! kgm? 
Water vapor pressure [51] PED Pa 


Dependencies are shown in fig. B.1-B.5 
^ Extrapolated from cyaon < 6M, but as shown in [27] extrapolated values fit very well to experiments 


cNaon € 12M. ? Extrapolated data from 343.15 K to 353.15 K. 


ion concentration also explains the in previous studies documented mass transport impact even 
at low current densities of about 0.1 kA m"? [39]. Water concentration only declines by 14.3 96 
to 43.5 x 10? mol m? at j = 8.AkA m ?, whereas the ion concentration increases by 37.3 % to 
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Table 4.3: Identified parameters. 


Name Symbol Unit Value 
Reaction kinetic rate constant ko mols~!m~2 165 
Specific interface Sphase,interface m? m~? 134 
Overall specific resistance Rspecific Qm? 1.9 x 10? 
Length of flooded agglomerates Zfa m 11.5 x 1076 
Length of gas phase ^ Zgas m 288 x 10-6 


^ Calculated by Zgas + Ztf + Zfa = Zelectrode 


16.1 x 10? mol m^?. Hence, significant consequences do not result from the lower water concen- 
tration being available for the reaction, but the mentioned change of the phase equilibrium. The 
change of water concentration at the interface of the flooded agglomerates to the liquid diffusion 
layer shows that not only the mass transport within the ODC is important, but also the mass 
transport across the liquid diffusion layer. 

For analyzing the driving forces of the mass flow in the liquid phase, the fraction of mass flow 
driven by diffusion Ne in relation to the total mass flow N; = Aum en of the 
species i is calculated by x; = |Nffusion| /|N;|. For all species the mass flow in the liquid phase is 
dominated by diffusion. As shown in figure 4.6, the oxygen mass transport in particular is almost 
100% diffusion driven, which is due to the steep concentration gradient and the low concentra- 
tion in the flooded agglomerated. Since oxygen depletes with increasing depth into the flooded 
agglomerates, an evaluation of the oxygen mass transfer is only directly behind the thin-film 
possible (cf. fig. 4.6). For water and hydroxide ions lower values are obtained: xy,0 ~ 70% 
and Xog- = 90% respectively. Since the concentration gradient of water and hydroxide ions are 
practically the same but in opposite direction within the flooded agglomerates, the diffusion mass 
flow is in the same range for both species as well. In contrast, the convection-based mass flow 
of each species is proportional to the concentrations, which is higher for the water than for the 
hydroxide ions. Accordingly, relative to the total mass transfer, the diffusion driven mass transfer 
for water is the lowest and for oxygen, the highest. The high ratio of convective mass flow of 
water and hydroxide at low currents close to the thin-film is caused by the evaporation of water. 
The concentration dependent equilibrium at the gas-liquid interface also affects the partial pres- 
sures and the mass transfer in the gas phase. The vapor pressure of the water decreases by 
a factor of 2.93 (cf. fig. B.3) and the inverse Henry constant increases by a factor of 5.81 
(cf. fig. 4.3) when increasing current to maximum. As a consequence of the changes equilib- 
rium for water evaporation and oxygen dissolution, the velocity of convection decreases from 
Vgas(J = 0) = 7.18 x 107? ms^! to vgas(j = jmax) = 1.78 x 10->ms—!. As a result, the partial 
pressure of water is decreasing with higher current densities. 

To get a deeper impression of the outlined findings and to separate the different effects, as well 
as getting information about their interaction, dynamic results are presented in the next section. 
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Figure 4.5: Dependence of simulated a) oxygen and b) water concentrations profiles within the liquid phase on current 
density. 


4.5.2 Dynamic analysis 


In this section, dynamic simulations with potential steps will be presented, which enable a separa- 
tion and identification of the interaction of the reaction and the different mass transfer processes. 
To cover the typical operation current densities of —4 to -6kAm~? [20], potential steps from 
E = —0.55V to E = —0.65 V and from E = —0.65 V to E = —0.55 V vs. OCP have been chosen, 
which correspond to steady state currents of —3.76 and —5.30kA m ?. The applied potential 
E and the current density j as the system response is depicted in figure 4.7. Two significant 
different time constants can be distinguished: A rapid and strong change of the current within 
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Figure 4.6: Fraction of mass flow driven by diffusion in comparison to the total mass flow in the flooded agglomerats as 
a function the current density at various locations. 


milliseconds and a slow relaxation of smaller amplitude within 5 s. The response signal to step- 
up and step-down in potential are identical, but with opposite sign. The rapid change observed 
for the positive as well as for the negative step results from the charging or discharging of the 
double layer. 
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Figure 4.7: Electrode potential and current dependent on the time for a potential step of —0.55 to —0.65 V vs. OCP at 
t — 10s and back at t — 20s. 
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The second observed, significantly larger time constant results from the mass transfer within the 
ODC following the consumption and production of species, but also from the potential sensitive 
charge transfer. In figure 4.8 a) it can be seen that after the potential step, both water and oxygen 
concentration need approximately five seconds to reach the new steady state, no overshooting is 
observed. Since the diffusion pathway of water respectively the hydroxide ions between reaction 
zone and bulk phase is much longer than the diffusion pathway of oxygen through the thin-film 
into the reaction area, at least two time constants for the mass transfer were expected. Indeed in 
figure 4.8 a) one may see that the oxygen concentration change is significantly steeper in the first 
second, but then slows down and reaches steady state at the same time as water. Comparing both 
plots in that figure, the concentration of water and oxygen at the boundary between thin-film and 
flooded agglomerates shows the same relaxation time as the current. This supports the conclu- 
sion, that the shift of the phase equilibrium couples both processes and dominates the dynamic 
response. 

To analyze the time behavior of oxygen mass transport, the oxygen partial pressure in the gas 
phase has to be taken into account (see fig. 4.8 b)). Oxygen partial pressure oddly rises simul- 
taneously with stronger current. Two aspects can be identified: Firstly, although more oxygen 
is consumed at higher current densities, the partial pressure of oxygen in the gaseous phase is 
increasing. Secondly the comparison of the dynamics in the gas phase with the dynamics in the 
liquid phase shows almost the same shape and time constant, suggesting similar rate determining 
steps. The rise in partial pressure can be attributed to the lowering of water activity due to the pro- 
duction of hydroxide ions and the consumption of water, which decreases the oxygen solubility, 
but also decreases the vapor pressure of the electrolyte. For deeper analysis, the time constants 
of the different mass transport processes can be approximated by equation (4.31), under the valid 
assumption that the mass transfer is dominated by diffusion (see. fig. 4.6). 


Az 
mt ^: a (4.31) 


Since the oxygen partial pressure in the gaseous phase is increasing with higher current densities, 
the sinking concentration of oxygen in the liquid phase at high current densities can be attributed 
to the mass transport through the liquid phase, but not by the gas phase. Inserting the thickness 
of the thin-film as diffusion length, a time constant of about Tmt,o, W 1.70 x 10~°s is calculated, 
which explains the significantly steeper fall of oxygen concentration than of water concentration 
in the first milliseconds after the potential step. For the mass transport of water and hydroxide 
ions between the reaction zone and the bulk phase, the overlapping time constants for the effec- 
tive diffusion in the porous structure of the flooded agglomerates eh uo and the free diffusion 
in the liquid diffusion layer T HO have to be determined separately. Here, [T =2.10s and 
under consideration of the porous structure (cf. eq. (4.19)), THO & 0.176s are calculated. Due 
to the more rapid decreasing of the water concentration within about 1 s after the potential jump 
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on the one hand and to the comparably long relaxation time to a new steady state of about 5 s on 


the other hand (cf. fig. 4.8 a)), it can be concluded that the observed current density response is 


dominated by both, the mass transport in the flooded agglomerates and the mass transfer in the 


liquid diffusion layer. 


Combining the given information it can be concluded, that the process governing the dynamic 


behavior is the slow mass transport of water respectively hydroxide ions through the liquid phase. 


Thus their interplay with the phase equilibrium dominates the overall performance of the ODC. 
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Figure 4.8: a) Concentration responses of oxygen and water at the thin-film-liquid interface for a potential step of —0.55 
to —0.65 V vs. OCP at t — 10s and back at t — 20s and b) response of oxygen partial pressure at various locations. 
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4.5 Results and discussion 


4.5.3 Parameter sensitivities 


To analyze the influence of the parameters for the evaluated electrode on the model output sys- 
tematically, a local sensitivity analysis has been performed. The relative sensitivity s; of the 
electrode potential Eg; to the parameter ç is determined as follows: 


Eexı—Eexı 
Eext 
Sc = 4.32 
[4 c—c? ( 3 ) 
Ge 


Where c? denotes the value at reference point and ç the one for an increased or decreased pa- 


rameter. Every parameter is varied + 596: ç = ç° - (1 +0.05). The results given in figure 4.9 
are the mean of individual results of the positive and the negative modified parameter. A high 
sensitivity of a parameter indicates a strong impact on the electrode potential, while a sensitivity 
close to zero indicates a negligible influence [76]. 

Sensitivity analysis is presented for three characteristic regimes (cf. fig. 4.4): Current densities of 
AkA m ? as the industrial operating point, 0.5 kA m7? as a kinetically governed operating point, 
and 8kA m? as a mass transport governed operating point were analyzed. For the evaluation, 
it is important to notice that the slope of the polarization curve (cf. fig. 4.4) at 0.5kAm ? and 
8kAm ? is much higher than the slope at 4AKA m ?. Consequently, the parameters may exhibit 
smaller sensitivities at 4kA m”? than at the other two points of the polarization curve. The varied 
parameters are further segmented into three different groups: Kinetic parameters, structural and 
geometric parameters and operation conditions. The results are given in figure 4.9. 

As expected, the reaction rate constant Ko affects the potential most at low currents. But although 
the performance at high current densities is dominated by mass transport, the performance shows 
significant sensitivity to high catalyst activity and a high catalyst surface. The symmetry factor 
a has a medium-high influence on the ODC over the entire operation range. 

To evaluate the sensitivity of the thickness of the different electrode sections individually, the 
boundary condition Zgas + zit + Zfa = Zelectrode Was ignored. The thickness of the liquid phase zır 
has the largest influence on the cell potential, while the length of the gas phase z, is not selec- 
tive. This supports the prior conclusion of a limitation caused by the mass transport of water and 
hydroxide ions across the liquid phase, while the mass transport in the gas phase is not limiting. 
ODC performance may thus be improved by a reduction of the flooded part, especially at higher 
currents. This could be reached by a more hydrophobic behavior of ODCs, using either smaller 
pores or a higher amount of the hydrophobic PTFE. But since a lower porosity and tortuosity due 
to smaller pores as well as a lower reaction rate due to a smaller catalyst surface would lead to 
a decrease of cell performance, an optimum of the mentioned parameters has to be found. Fur- 
ther, the length of the flooded part zr, can not be arbitrarily reduced, because a certain space for 
large specific interface is needed. The same holds for the diffusion layer in the liquid electrolyte: 


53 


4 Processes and Their Limitations in Oxygen Depolarized Cathodes 


i = BE: = 0.5 kam?) 
Ms 6 eam) 


Sj — 8kA- m?) 


S interface | 


1 
i a 
- S 

4 


Zab 


-1 -0.5 0 0.5 1 1.5 2 
Parameter sensitivity S, l- 


Figure 4.9: Relative sensitivity sç of the electrode potential to parameter variation of + 5% for three different constant 
current densities. When increasing cwaog, the current density of j = 8 kA m ? could not be reached, thus only results for 
negative variation are shown. 


A thicker liquid diffusion layer has especially at higher current densities a negative impact to 
operation conditions of the electrode, here an approach to improve the electrode performance is 
given by an forced convective flow to reduce the diffusion layer thickness zig. Moreover, the 
thin-film acts as an additional diffusion resistance for oxygen dissolved by the gaseous phase 
[38], thus the high sensitivity of the thin-film thickness at the diffusion limiting current densities 
is explainable. It has been shown that especially the thin-film thickness dominates the shape of 
the polarization curve [45]. Finally, the specific interfacial area Sphase,interface between liquid and 
gas phase is particularly important at high current densities, where oxygen solubility limits the 
performance. 

The large sensitivity to the electrolyte concentration cNaon (z = z;) emphasizes the previously dis- 
cussed solubility limitation by the ion concentration at the gas-liquid-interface. Previous studies 
have shown that the highest current densities for the ORR could be achieved with electrolyte 
concentrations of 5 M [70]. I suggest that process conditions with high concentrations over 5 
M leads to a mass transport limitation. But as the hydroxide ions are the desired product of the 
ORR in advanced chlor-alkali electrolysis, the only way for a higher efficiency of the process is a 
faster removal of ions. The potential is sensitive to oxygen partial pressure. The decrease of the 
electrode potential with increasing pressure Pgas results from lower mass transfer limitations of 
the oxygen due to the higher amount of dissolved oxygen into the electrolyte. Although a better 
performance with higher pressure is documented [70], it has to be mentioned that the pressure 
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of the oxygen can only be increased in a defined range to avoid a break through of the gas to the 
electrolyte side [9]. Also notable is the sensitivity for the pressure even at low current densities, 
which meets the results of former studies, that the mass transport influences the ODC losses [39] 
even at kinetically controlled current densities. 


4.6 Conclusions 


The first dynamic three-phase model of a porous ODC for advanced chlor-alkali electrolysis has 
been presented. Suitable model parameters were partly obtained from literature and partly iden- 
tified from published polarization curves. 

It was shown that at industrial operation conditions, the ORR only takes place close to the gas- 
liquid interface, because of depletion of the dissolved oxygen. The consumption of water and the 
production of hydroxide ions leads to a significant change of the phase equilibrium with current 
density, which causes a lower solubility of the oxygen into the liquid electrolyte. This effect 
intensifies at higher reaction rates and causes a total depletion of oxygen and thus the limiting 
current. Water is sufficiently available for the reaction and does not limit the maximum current 
density in a direct way. A deeper dynamic evaluation of this behavior pointed out that the re- 
moval of ions from the reaction area is the slowest process. It has a significant impact not only 
on stationary performance but also on the full dynamic behavior of all concentration profiles and 
thus of the dynamic response. In this context, it was also shown that only a few micrometers 
of the ODC are flooded with the liquid electrolyte under industrial working conditions. Hence, 
dynamic analysis is a highly potent method to reveal the interaction of processes in ODCs and to 
determine its state. Dynamic measurements should be conducted for reliable, unique parameter- 
ization and for further validation of the model. 

A sensitivity analysis of the kinetic, geometric and process parameters indicated that an improve- 
ment of the ODC performance can be affected most with a fast mass transport in the liquid phase, 
e.g. via a smaller level of flooding within the electrode or a smaller liquid diffusion layer. 
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Transport and Water Activity' 


In this chapter, the previously presented basic model is enlarged and used to analyze dynamic 
and steady-state measurements, in order to determine the role of water activity and water mass 
transport on ODC performance and dynamics. 


5.1 Introduction 


In the previous model-based study in chapter 4, I analyzed the limiting processes of ODC. I 
have shown that the consumption of water and the production of hydroxide ions induced by 
the electrochemical ORR lead to an enrichment of hydroxide ions at the gas-liquid interface 
within the liquid electrolyte. This causes a reduction in availability of water, a change of the 
gas-liquid equilibrium and in consequence a substantial decrease of oxygen solubility. Based 
on these results I hypothesized that a fast mass transport of water and hydroxide ions across the 
porous electrode and the electrolyte diffusion layer outside the electrode are key drivers for a 
high ODC performance. Those findings are supported by experiments of Clausmeyer et al. who 
achieved higher current densities in a 5 M compared to a 10M NaOH electrolyte utilizing single 
nanoparticle studies [70], and by Botz et al. who measured an accumulation of hydroxide ions 
in the diffusion layer in close proximity to the ODC surface at higher current densities [34]. 
Furthermore, Zhang et al. presented a decreasing oxygen diffusivity in line with an increasing 
NaOH concentration and a decreasing water activity as a further limiting effect in ODCs [27]. 

The aforementioned studies demonstrate that the electrolyte has a significant influence on the 
ODC performance. Under careful consideration of the non-ideal behavior of water as a solvent, 
the impact on the electrocatalytic properties of an ODC may be explained: In highly concentrated 
aqueous electrolytes, the activity of water, and its thermodynamics [77], cannot be described by 
dilute solution theory. Instead, the activity coefficient strongly depends on the ion concentration 
[75]. The underlying mechanism of a decrease in water activity at increasing NaOH concentra- 
tion is the formation of solvation shells, in which water molecules are bound to ions. In alkaline 


7 Parts of this chapter have been published in Róhe et al., ChemElectroChem, 2019 [2]. 
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media, for instance, the first layer of the hydration shell of a hydroxide ion consists on average 
of four water molecules, whereas the coordination number decreases to about three at higher ion 
concentrations [78]. A decreased water activity goes, apart from the lower solubility of oxygen 
[50], along with a change in mass transport [79], a lower vapor pressure of the water [51] and 
slower reaction rate [27]. To completely elucidate the influence of the water activity on the ODC 
performance, a thorough analysis is required. 


In this chapter I will confirm the hypothesis that water activity and water mass transport rep- 
resents a crucial factor in governing the ODC performance. For this purpose, the model was 
extended in order to resolve the influence of the water and hydroxide activity more precisely. 
Stationary and dynamic measurements in two different electrochemical setups were conducted 


by cooperation pats and analyzed with the extended one-dimensional dynamic three-phase 
model measurements. Both setups had comparable conditions concerning the oxygen gas supply, 
yet they differed in the respective electrolyte conditions: The first experimental setup was op- 
erated under convective electrolyte flow, whilst a stagnant electrolyte was present in the second 
one. Additionally, in the latter setup, the analysis was assisted by scanning electrochemical mi- 
croscopy (SECM), which allows the determination of hydroxide and water activities at a distance 
of about 1 um above the ODC surface. The differences in the dynamic behavior and steady state 
performance of both experimental cells are analyzed using the dynamic model. By comparing the 
two cells, the influence of liquid phase mass transport and water activity in the liquid phase was 
assessed model-based. It is shown, that both systems differ significantly in their performance and 
dynamics. The model-based analysis of dynamic measurements and the SECM measurements 
demonstrated that the performance and the dynamics of ODCs are highly dependent on the mass 
transport of water and hydroxide ions between the bulk phase and the electrode. 


5.2 Experimental 


Two different electrochemical cells, as schematically depicted in fig. 5.1, were utilized and both 
were operated using identical ODCs (see table 5.1). The measurements have been performed 
by the Institute of Chemical and Electrochemical Process Engineering, Technische Universitát 
Clausthal, Germany and Analytical Chemistry — Center for Electrochemical Sciences, Ruhr 
University Bochum, Germany. In this section, the ODC preparation, measurements and the 
experimental setups are summarized shortly. A detailed description is given in the appendix C. 


5 Measurements performed by Institute of Chemical and Electrochemical Process Engineering, Technische Universität 


Clausthal, Germany and Analytical Chemistry — Center for Electrochemical Sciences, Ruhr University Bochum, 
Germany, published in [2]. 
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In both setups polarization curves as well as potential steps were measured in a 10M NaOH 
electrolyte at 50°C. The main factor differentiating both setups was the mass transport at the 
liquid side offered in the respective system. While the electrolyte was driven by a convective flow 
in the first setup (con), it was stagnant in the second one (stag). In the first cell, the electrolyte 
temperature was controlled by regulating the temperature of the external electrolyte reservoir. 
In the second cell, an external heating circle had to be installed to enable temperature control 
without electrolyte flow. Furthermore, both setups differed in their geometrical cross section and 
in their electrode orientation respective to the ground (cf. fig. 5.1). A detailed documentation 
of the setup with the convective electrolyte flow is given in ref. [63], and for the cell with the 
stagnant electrolyte in ref. [34]. In the latter setup, SECM was used to quantify the hydroxide 
activity under stationary operating conditions. 


Table 5.1: Electrode characteristics. 


Name Symbol Unit Convective Stagnant elec- 
electrolyte cell trolyte cell 

Electrode thickness Zelectrode m 300 x 10-6 295 x 10-8 

Geometrical cross section Ageo cm? 3.14 0.264 

Catalyst loading OAg mgcm ? 130.0 129.9 

PTFE loading OPTFE mgcm ? 4.0 4.1 

OCP (measured) Eo V vs. RHE 1.095 1.10 


5.3 Modelling and parameterization 


In this section, first the expansion of the model is described. Subsequently the parameterization 
of the model using the dynamic and staionary measurements is given. 


5.3.4 Expanded liquid mass transport model 


The equations of the basic model have been introduced in the previous chapter 4. Since it is 
known from the last chapter, that ion concentration and water activity are crucial for the ODC 
performance, their influences are investigated by means of the expanded liquid mass transport 
model. Therefore following extensions have been made: 
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Figure 5.1: Schematic view of the used measuring cells with a) convective electrolyte flow and b) stagnant electrolyte. 


* The ion concentration, water activity and also the oxygen availability affect the thermody- 
namics of the ODC and thus OCP. Considering this, the Nernst-equation (see eq. (5.4)) 


has been included into the model. 


e To solve the Nernst-equation, the NaOH activity anaou is needed. It is calculated based 
on the mean activity coefficient y+ of the electrolyte, assuming electro-neutrality, cf. eq. 


(5.6) and (5.7) [66]. 


* The oxygen activity ao, is modeled to be non-ideal and have been set into dependency of 


the ion concentration of the electrolyte, cf. eq. (5.5) [50]. 


The differences to the basic model are also summarized in table 3.1. All equations implemented 


in the expanded liquid mass transport model are given in table 5.2. 
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Table 5.2: Model equations for expanded liquid mass transport model. (Model schema for reference: fig: 4.1.) 


Electrochemical reaction kinetics, double layer capacitance and internal resistance 


Faraday current Jp--F4r (5.1) 
1— a). F- i 
Reaction rate r = ko : ao, ao “exp ( ( Hn ise ) (5.2) 
Reaction overpotential Nreaction = Eim — Eo (5.3) 
R.T 40, 44,0 
Nernst equation Eo = Eoo + 4.F In ( A ) (5.4) 
NaOH 
Activity of O, 40, = X0, WO, (5.5) 
. bwaoH 
Activity of NAOH @NaOH = pg — (5.6) 
NaOH 
NNaOH 
Molality of NaOH DNaOH on (5.7) 
Ohmic drop Eext — Eint = R specific ^ J (5.8) 
dE int ; 
Charge balance Cai d —-j-F4r (5.9) 
t 
Mass transport in the gas phase and processes at the gas-liquid interface 
dpi 9? pj gas,eff dpi 
Species balance "e ee dE ‘Dee FE (5.10) 
with i € (H50,05] 
gas, ff s € 
Effective diffusion coefficient gr me gr 2 (5.11) 
Bruggemann correlation qug to (5.12) 
x Po, 
l ds 2 (5.13) 
Henry’s law > Ho, ( a um ch " 
Protal i ioni — M 
Stefan M RT Vegas = Nao — NGissolution (5.14) 
A hacian 
Specific phase interface area Sphase,interface = a. (5.15) 
geometric 
Mass transport in the thin film 
: de &/8ci cj ec c ia, OC; 
Species balance m ( tor J * ) liq. € IL ae 
i ot L Oz* Crta 022 Ctotal ijt gy  Hlevaporation 
ij 
with i € {OH ,H,0,0,}, j € {OH ,H,O} (5.16) 
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iin velocity due to evapo- iier =0 (5.17) 
Mass transport in the flooded agglomerates 
2. ; 
Species balance a e- (Fo Fy) otto 
ij 
with i € {OH , H20, 02}, je {OH ,H,0} (5.18) 
Velocity due to reaction Vliq,react (2) = (vo: = [VSRR y [ i A dz (5.19) 
Total velocity of liquid phase Vliq (z) = Vliq,evaporation + Vliq, react (Z) (5.20) 
Mass transport in the liquid diffusion layer 
Species balance un = Y EE ej 2 E ) Du (521) 
t EMO Coa 0% Cowi J az 


ifj 
with i € (OH^,H50,05), j € (OH”,H,0) 


Velocity Vid] = Vliq (z = zı) (5.22) 


Boundary conditions 


O, in gas bulk Po, (z = 0) = Pia (5.23) 
H,O in gas bulk Pu,0(z = 0) =0 (5.24) 
O; in electrolyte bulk co,(z =z) =0 (5.25) 
NaOH in electrolyte bulk CNaOH (Z = 21) = cNaon(z = 21,1 = 0) (5.26) 
H,O in electrolyte bulk“ c,0(2 = a) = f(CNaon(2 = 21,1 = 0)) (5.27) 
Electrode geometry Zgas + Ztf + Zfa = Zelectrode (5.28) 


@ Calculated as in [51] 


5.3.2 Parameterization and parameter evaluation 


Constant model parameters are adjusted to the operating conditions of the measurements and are 
listed in table 5.3, concentration dependent parameters are given in table 5.4. For the manual 
parameter identification, dynamic potential steps from 0.8 to 0.7 V vs. RHE and back to 0.8 V as 
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well as an experimental polarization curve coupled with concentration measurements were used. 
The experimental data and simulation results will be discussed in the next section. The agree- 
ment between experiment and simulation is satisfactory for the polarization curves as well as the 
potential steps. The parameter values for both cells are presented in table 5.5. The combination 
of stationary and dynamic measurements enables one to cover a wide operation range but also 
to determine the different processes in the electrode separately due to their characteristic time 
constants. 


As expected, significantly different values for both systems are obtained for the length of the 
liquid boundary layer zi in and the electrolyte bulk phase: zjj^ = 10 um for the setup with the 
convective electrolyte flow and rd = 270um for the setup with the stagnant electrolyte. For 
the reaction rate constant Ko, the specific interface Sphase,interface and the average length of flooded 
agglomerates zg, both systems exhibits values in the same order of magnitude. The larger specific 
interface in the cell with the convective electrolyte flow may occur due to the higher motion in 
the liquid electrolyte and is associated with temporary local shifts of the gas-liquid interface on 
the pore scale level. 

The large differences in double layer capacitances may be explained by its dependency on the 
current densities [72], which differs considerably in both systems at the respective potentials. 
Additional possible measurement inaccuracies are discussed below. 

In my previous work dealing with an ODC operated at 80?C, I presented a higher value of 
Sphase,interface = 134m? m ? and a slightly higher value for zr, of 11.5um [1]. On the one hand 
this deviation may originate from temperature dependent system characteristics like viscosity or 
wettability. On the other hand, the use of dynamic measurements in my present work allows for 
a better separation and a more precise identification of the process parameters. 


5.4 Steady state analysis of oxygen depolarized 
cathodes performance 


Figure 5.2 displays the measured and the simulated polarization curves of both systems. In the 
kinetically governed regime at small current densities both cells have almost the same behavior, 
which is expected, since similar electrodes were used. Accordingly, the derived reaction rate 
constants Ko (cf. table 5.5) are in the same order of magnitude. 

In contrast, at medium current densities up to the diffusion limited current density, a strong dif- 
ference between both systems occurs: While the system with the stagnant electrolyte shows a 
significant mass transport limitation at about j ~ 1.3kA m"? at Eg = 0.1 V vs. RHE, no mass 
transport limitation proceed for the system with electrolyte convection at j = 3.8kA m ? at the 
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Table 5.3: List of constant parameters and operating conditions. 


Name Symbol Unit Value 
Electrode geometry 
Porosity [32] € — 0.4 
Operating conditions 
Pressure in gas chamber“ P= Po, Pa 1.013 x 105 
NaOH concentration in electrolyte bulk“ CNaOH(Z = zr) molm— 1.00 x 104 
H50 concentration in electrolyte bulk [51] cH,0(z = 21) molm” 5.05 x 10* 
Temperature” T K 323.15 
Electrochemical and kinetic data 
Charge transfer coefficient” [24] a — 0.15 
Stoichiometric coefficient H,O VORI — —2 
Stoichiometric coefficient OH VORR = +4 
Binary diffusion coefficients 
For gaseous oxygen/water [73] Z5, H0 m?s! 2.48 x 10? 
For liquid oxygen/water® [73, 74] Er uo ms! 1.36 x 107? 
For liquid oxygen/NaOH“ [73, 74] BO aH m?s”! 2.02 x 107? 
For liquid water/NaOH? [60] ol mes" 5.33 x 10710 


^ Experimental operating condition. P Measured with 6.5 M NaOH. ^ Calculated with modified Wilke 
and Chang equation, with O = 3.9 and a = 0.5 chosen, adjusted to experimental data from [74]. ^ No 
NaOH data available, modeled for KOH with NaOH viscosity from [68], as seen in [32]. 


Table 5.4: Electrolyte concentration dependent variables. 


Name Symbol Unit 
Inverse Henry constant^ [22] Ho, Pam? mol! 
Activity coefficient of O, » 150] Yo, — 

Activity of H,O [75] 4,0 — 

Mean activity coefficient of NaOH [66] Yi kgmol^! 
Density of electrolyte [68] p kgm ? 
Water vapor pressure [51] Pu Pa 


7 Extrapolated from cnaoH < 6M, but as shown in [27] extrapolated values fit very well to experiments 
with cyaom € 12M, b Extrapolated from cyaon < 5.5 M, but as shown in original research, calculations 


for e.g. cog € 13.5M are valid [50]. 


same potential. 
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Table 5.5: Parameters identified from electrochemical experiments. 


Name Symbol Unit Convective Stagnant elec- 
electrolyte trolyte 

Reaction rate constant ko mols Im”? 15.42 x 104 9.20 x 104 
Specific interface Sphase, interface m? m~? 56.0 36.4 
Liquid diffusion layer Zl m 10 x 1076 270 x 1076 
thickness 
Overall specific resistance Rspecific Qm? 1.16 x 1074 1.82 x 10-4 
Length of flooded Zfa m 7.90 x 10-6 8.14 x 10-6 
agglomerates 
Thin-film thickness ^ auf m 80 x 107? 
Double layer capacitance Cal Fm? 80 720 


“ Set as equal for both systems. Typical values for the thin-film thicknesses: 25 - 300nm [45]. 
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Figure 5.2: Simulated polarization curves in comparison to experimental polarization curves for convective electrolyte 


flow and stagnant electrolyte cell. 


The ORR can proceed at locations within the flooded agglomerates where the reactants oxygen 


and water both are present. In the following, this section is called reaction zone. In figure 5.3, 


simulated oxygen concentration profiles over the thin-film and the first micrometer of the flooded 


agglomerates for both systems at different current densities are shown. For both systems the 


reaction zone, facing the gas-liquid interface, decreases in size with increasing current densities 


due to the depletion of oxygen. Similarly to what was presented earlier [1], no oxygen associated 


limitation takes place due to the mass transport in the gas phase, therefore all oxygen mass 


transport limitation can be attributed to the liquid phase. 
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Figure 5.3: Dependence of simulated oxygen concentrations profiles within the liquid phase with respect to the current 
density for a) convective electrolyte flow and b) stagnant electrolyte. 


In the system with the stagnant electrolyte, the reaction zone shrinks to almost zero at the mass 


transport governed current density of j ~ 1.3kAm ? 


. This shows that the limiting current density 
is determined by the availability of the oxygen. In contrast, for the system with the convective 
electrolyte flow, the reaction zone amounts to several nanometers even at a current density of 
j — 3.8kAm ? so that no oxygen mass transport limitation occurs. 

The vast disparity in oxygen availability is caused by the difference of current dependent water 
activity inside the electrodes. The ORR consumes water and produces hydroxide ions within the 
reaction zone close to the gas-liquid interface. Since three to four water molecules form a solva- 
tion shell around the hydroxide ion [78], the amount of free water molecules and thus the water 


activity decreases significantly due to the production of hydroxide ions. This change of water 
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Figure 5.4: Simulated and experimental NaOH concentration in the liquid diffusion layer 1 um above the ODC surface 
for cell with stagnant electrolyte. For comparison, NaOH concentration in the cell with the convective electrolyte flow is 
also given at the same location. 


activity leads to a change of the phase equilibrium of water and oxygen, and thus to a reduced 
oxygen solubility [50] and a higher vapor pressure of water [51]. There is also experimental evi- 
dence for the accumulation of hydroxide ions in the system with the stagnant electrolyte: Figure 
5.4 shows the hydroxide ion concentration detected via SECM 1 um above the ODC surface (cf. 
fig. 5.1), in comparison to the simulated values at the same location. It can be deduced that both, 
the measurement as well as the model, show a significant increase of hydroxide concentration 
and subsequently a decrease of the water activity. Differences between experiment and simula- 
tion amounts to maximum 10 46 and might result from additional convection caused by a density 
gradient — due to temperature or ion concentration gradients — or local inhomogeneities of ion 
concentration in-plane direction of the electrode in the experiment [34]. The mass transport re- 
sistance for the exchange of water and hydroxide ions between the ODC and the electrolyte bulk 
is much higher for the thick liquid diffusion layer in the system with the stagnant electrolyte than 
in the system with the convective electrolyte flow. It may be that the mass transport in the liquid 
diffusion layer is not only influenced by the forced convection, but also by the vertical orienta- 
tion of the ODC. A thicker diffusion layer constitutes a larger diffusion resistance and causes a 
higher accumulation of hydroxide ions in the electrode. Thus, the accumulation of the hydroxide 
has a detrimental effect on the oxygen solubility, which is higher in the system with the stagnant 
electrolyte than in the system with the convective electrolyte flow. This in turn explains why the 
system with the stagnant electrolyte is limited by oxygen availability at E — 0.1 V vs. RHE. 
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5.5 Dynamic analysis 


In this section the influence of the reaction-induced oxygen solubility limitation on the dynamics 
of the ODCs is discussed in detail based on dynamic measurements and simulations of both 
configurations. In figure 5.5, the current response to a potential step from 0.8 to 0.7 V vs. RHE 
and back is displayed for both systems. 
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Figure 5.5: Simulated and experimental time-dependent current density response for a potential step from 0.8 to 0.7 V vs. 
RHE at t = 30s and a) back at t = 60s for convective electrolyte flow cell and b) back at t = 90s for stagnant electrolyte 
cell. 
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Although the same potential is applied to both electrodes, the system with electrolyte convection 
offers a higher cathodic current density, which agrees with the polarization curves and has been 
discussed in the previous section. Furthermore, both systems exhibit a very different dynamic 
response to the potential step, which will be discussed individually. Two distinct processes can 
be separated: The charging and discharging of the double layer and the mass transport of water 
and hydroxide ions across the liquid phase. As discussed below, the dynamics of the oxygen 
follow the latter process. 

The charging and discharging of the double layer is the reason for the positive and negative cur- 
rent peaks directly after the positive and the negative potential step. For both systems, different 
time constants are obtained, which is attributed to the dependency of the double layer capacitance 
on the current density [72]. Additional a lager area outside the reactive area due to small leakages 
between the sealing ring and the ODC of the double layer may be penetrated in the smaller cell 
with the stagnant electrolyte. Since the electrode circumference / area ratio is approx. 3.5 times 
larger for the smaller cell than for the lager cell with the convective electrolyte flow, this effect 
may be more significant. 

The subsequent relaxation of the current density occurring in the range of seconds is dominated 
by the mass transport of water and hydroxide ions between the reaction zone and the electrolyte 
bulk phase [1]. The corresponding time constant can be approximated with equation (3.27), 


derived from Fick's second law: j 
Az 


(5.29) 


Tmt © 


With Tmt as the time constant of water and hydroxide ion mass transport, D as the diffusion co- 
efficient and Az as the mass transport length. For both systems the mass transport time constant 
is a function of the diffusion through the porous structure of the flooded agglomerates and the 
free diffusion through the liquid diffusion layer (cf. fig. 4.1), which can be calculated sepa- 
rately. Since the length of the flooded part of both electrodes is almost the same, similar time 


constants are obtained: e mo = 0.19s for the convective electrolyte flow and rf aes = 0.20s 


m 
for the stagnant electrolyte. In contrast, the time constants for the mass transport across the liq- 


uid diffusion layer differs by orders of magnitude between both systems: Qum s — 0.19s and 


m 
HN = 181.3s. The large time constant is in line with the observation that no steady state is 


reached within the 60s duration of the dynamic test for the system with the stagnant electrolyte, 
but within 1s for the system with the convective electrolyte flow. In conclusion, the water and 
hydroxide ion transport across the liquid diffusion layer dominates the current response to the 
potential step and causes the large difference in the current response of both electrodes. The cor- 
responding time-dependent ion concentration and water activity at the gas-liquid interface and 
1 um away from the ODC surface are presented in figure 5.6. 
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Figure 5.6: Simulated NaOH concentration and water activity within the liquid phase dependent on the time for a 
potential step from 0.8 to 0.7 V vs. RHE at t = 30s and a) back at t = 60s for convective electrolyte flow cell and b) back 
at t = 90s for stagnant electrolyte cell. 


Due to the relationship between water and hydroxide activity as discussed above, the activities of 
both species show opposite trends. Furthermore, water activity has a strong influence on oxygen 
solubility [75]. This effect is illustrated in figure 5.7 for the system with the stagnant electrolyte. 
The right-hand side of the plot displays the steady state water activity gradient within the liquid 
diffusion layer for t = 30s directly before the negative potential step and at selected time points 
afterwards. The plot in the middle shows the same process within the flooded agglomerates. 
A comparison of the curves suggests that the gradient over the whole electrode and the liquid 
diffusion layer change with a delay. This is due to the fast reaction rate compared to slow mass 
transport of water and hydroxide ions. The plot on the left-hand side of figure 5.7 depicts the 
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Figure 5.7: Henry constant at the gas-liquid interface (left) at steady state current at 0.7 V vs. RHE and its change at 
selected time points after a potential step to 0.8 V vs. RHE. Corresponding gradient of water activity within the ODC 
across the flooded agglomerates (middle) and in the liquid diffusion layer (right). Cell with no electrolyte convection. 


corresponding change of oxygen solubility for the different points in time, indicated by the rela- 
tive change of the Henry constant at the gas-liquid interface during the step. The change of the 
Henry constant directly affects the oxygen solubility. In the cell with electrolyte convection the 
liquid diffusion layer is smaller and the new steady state of the phase equilibrium and thus of 
oxygen concentration is reached much faster. 

The resulting response of oxygen concentration with and without electrolyte convection is 

depicted in figure 5.8 for different locations. It can be inferred that due to the ion dependent 
solubility the oxygen concentration is governed by the same time constant as the ion concentra- 
tion. Therefore, the oxygen supplied to the reaction zone and thus the current density exhibits 
the same time constants as the liquid phase mass transfer. 
In figure 5.5, deviations between experiments and simulations are visible. For the system with 
electrolyte convection, a slight decrease of the negative current density with time after the nega- 
tive potential step at 30s can be identified in the experiment. For the system without electrolyte 
convection, the measured current density is more cathodic than the simulated one. Different 
factors which are not included in the model might have contributed to this behavior: Although 
the experiments were conducted under temperature control, a moderate increase of the ODC 
temperature which would improve kinetics and transport processes, cannot be excluded. In that 
regard, heating effects can be caused by ohmic losses or the entropy effect, for instance [8]. The 
same effects may occur for the steady state results shown above. 
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Figure 5.8: Simulated oxygen concentration within the liquid phase dependent on the time for a potential step from 0.8 
to 0.7 V vs. RHE at t = 30s and a) back at t = 60s for convective electrolyte flow cell and b) back at t = 90s for stagnant 
electrolyte cell (fa: flooded agglomerates). 


From the dynamic analysis, I can conclude that the performance as well as the dynamics of 
ODCs are determined by the mass transport across the liquid phase. This insight is corroborated 
by two main aspects: On the one hand, in the system with the stagnant electrolyte, water activity 
decreases only by less than 7.5% as a consequence of the negative potential step, but oxygen 
solubility decreases by 18% (cf. fig. 5.7). This suggests a minor influence of the water activity 
on the kinetics in a direct way but a high indirect impact due to the oxygen solubility. On the 
other hand, due to the faster water mass transport, the oxygen solubility of the system with the 
convective electrolyte flow only decreases by 7%. The comparison of both systems indicates 
that the mass transport within the liquid electrolyte has a significant influence on the electrode 
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performance. It is expected that also electrode design and temperature will have significant 
impact, as they are known to affect transport phenomena. This leads to the conclusion that 
the performance of technical ODCs can be improved by accelerating the mass transfer in the 
liquid phase, e.g. via convective electrolyte flow or electrolyte flow field optimization. Already 
now, in industrial application quite high performances are reached due to convective electrolyte 
transport in so-called percolators [3], further performance increase might be reached by a even 
further optimized mass transport. 


5.6 Conclusions 


Experimental polarization curves and potential steps demonstrate that the performance and dy- 
namic behavior of ODCs deviate strongly when the electrodes are operated with and without 
convection of the liquid electrolyte. Without electrolyte convection, a mass transfer limitation is 
reached at a current density of j 1.3 kAm~?, whereas this does not occur at current densities 
of j > 3kAm ? with electrolyte convection. Furthermore, the dynamic response of the current 
density to a potential step exhibits a significantly slower relaxation behavior without electrolyte 
convection. 

By comparing both systems utilizing simulation and measurements, the hypothesis that water 
activity plays a key role was confirmed. The accumulation of hydroxide ions leads to a reduction 
in water activity and thus in oxygen solubility. Since three to four water molecules are bound to 
every hydroxide ion as a solvation shell, the water activity decreases rapidly as a function of the 
ionic strength, and the accompanying change of the phase equilibrium leads to a lower oxygen 
solubility. The model-based analysis suggests that electrolyte convection increases mass trans- 
fer rates in the liquid phase which in turn reduces the hydroxide ion accumulation and thereby 
increases both, water activity and oxygen solubility. Without electrolyte convection, a thicker 
liquid diffusion layer is established causing a stronger accumulation of hydroxide ions and a 
lower water activity. Not all processes can be explained by mass transport and phase equilibrium 
phenomena, so that further research on catalyst activity and local reaction rates within porous 
structures is needed. 

In conclusion, a high mass transport rate of water and hydroxide ions is essential for high ODC 
performance and needs to be ensured for technical applications of these kind of electrodes. A 
thorough understanding of the transport processes is also indispensable for understanding the 
dynamic behavior of the electrodes. 
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Depolarized Cathodes Structure 
and Electrolyte Distribution’ 


In this chapter, a modeling approach of ODCs to evaluate the influence of inhomogeneities in 
electrolyte distribution on electrode performance and dynamics is presented. It is investigated 
how the location and size of the gas-liquid interface affect the local reaction rates within the 
ORR. 


6.1 Introduction 


Recently, Franzen et al. [23] analyzed in their experimental study ten ODCs with different com- 
positions and structures, and found significant performance differences. In detail, they evaluated 
ODCs with a silver content of 90 to 99 wt—% (10 to 1 wt—% PTFE, respectively) and identi- 
fied a silver content of 98 wt—% as best. The amount of silver had a very low influence on the 
polarization curve at low, e.g. kinetically governed current densities, but a significant influence 
on the performance at high current densities, i.e. in the mass transport dominated range of the 
polarization curve. It was hypothesized that the hydrophilic ODC with a high silver content of 
99 wt—% was flooded by the liquid electrolyte, whereas ODCs with a high PTFE content above 
5 wt—% were too hydrophobic and electrolyte intrusion was insufficient. Both effects may lead 
to a decrease in the extent of the gas-liquid interface and thus an insufficient oxygen availability 
in the liquid phase [23]. 

Paulisch et al. [80] investigated the electrolyte distribution within ODCs by operando X-ray ra- 
diography and chronoamperometric measurements. They showed that the pore structure as well 
as the distribution of the electrolyte during operations is inhomogeneous: some pores were barely 
flooded by the electrolyte, whereas a local breakthrough of the electrolyte to the oxygen side was 
observed in other pores [80]. These findings suggest that the gas-liquid interface may also be 
inhomogeneously distributed over the entire thickness of the electrode. The identified locally 


9 Parts of this chapter have been published in Róhe et al., Electrochimica acta, 2021 [15]. 
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high electrolyte saturation seems surprising, as previous studies have shown that even very small 
amounts of PTFE lead to a non-wettability of a pore system [81]. 

So far in this dissertation, it was demonstrated that the ODC performance is limited by the avail- 
ability of dissolved oxygen in the liquid phase, whereby the oxygen solubility rate is significantly 
affected by the length of the water and hydroxide mass transport path between the gas-liquid in- 
terface and the liquid bulk. 

The findings of Franzen et al., Paulisch et al. and the results of this work raise new questions: 
How does the phase distribution affect the ODC dynamics and performance? Does the distributed 
mass transport length of water and ions in the liquid phase lead to locally different oxygen sol- 
ubilities and thus to different local reaction rates? Is this effect related to the silver / PTFE ratio 
dependent ODC performance? To answer these open-ended questions and to quantify the in- 
fluence of the electrolyte distribution, a model is needed that considers the inhomogeneities in 
electrolyte distribution within ODCs. The importance of understanding the distribution of the 
liquid phase in porous electrodes was already emphasized by various experimental and stimula- 
tive studies inter alia for proton-exchange membrane fuel cell (PEMFC) [82-84]. 

Model-based studies to investigate inhomogeneities were inter alia performed by Levi et al. who 
connected two different diffusion paths in parallel to study diffusion time constants for particle 
size distributions in non-uniform composite graphite electrodes for lithium-ion batteries [85] or 
by Jüttner and Lorenz who investigated the influence of surface inhomogeneities on corrosion 
processes for different metal electrodes in sodium hydroxide electrolyte [86]. Chevalier et al. 
interconnected single cell models to draw conclusions about the local distribution of the state-of- 
health of a PEMFC stack from the stack impedance spectra [36]. Krewer et al. have shown that 
local inhomogeneities in concentration in direct methanol fuel cells (DMFC) lead to characteris- 
tic changes of the dynamic behavior [87]. 

In this chapter, I present the first model for describing the influence of inhomogeneities in elec- 
trolyte distribution on ODC performance. To this end, the in chapter 5 introduced dynamic 
expanded liquid mass transport model is extended to incorporate multiple electrode domains, 


each with different representative characteristics. The inhomogeneous ODC model is used to 
10 
analyze experimental electrochemical impedance spectra (EIS) as well as polarization curves 


of ODCs with systematically varied silver catalyst content. It will be shown how to interpret 
the observed flattened semicircles in the impedance spectra regarding to local inhomogeneities. I 
demonstrate how the structure and silver content of ODCS affect the electrolyte distribution, elec- 
trode performance and dynamics. The aim of this study is to infer the macroscopic distribution 
of the electrolyte from the performance and impedance spectra of the ODCs. The methodology 


10 Measurements performed by Institute of Chemical and Electrochemical Process Engineering, Technische Universität 
Clausthal, Germany, published in [15]. 
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used does not allow conclusions to be drawn about microstructures such as pore size distribu- 
tion or local electrolyte distributions. It will be discussed that not only the size of the gas-liquid 
interface but also the location is of importance for the performance. These are key findings to 
understand the influence of the electrode structure on the electrode performance and dynamics. 
EIS is shown to be a highly sensitive tool to diagnose the electrolyte distribution within ODCs, 
while this information is hardly accessible with potential steps or polarization curves. 


6.2 Experimental 


The measurements have been performed by Institute of Chemical and Electrochemical Process 
Engineering, Technische Universität Clausthal, Germany. Detailed descriptions of the exper- 
iments and also the electrode preparation are given in the appendix D and by Franzen et. al 
[23]. In short: The electrochemical experiments have been performed in a Gaskatel half-cell 
(FlexCell HZ PP01, Gaskatel GmbH) with a geometrical area of 3.14 cm?. After a start-up pro- 
cedure, pseudo-galvanostatic impedance measurements have been performed at current densities 
of 2kAm ? with a 10 mV amplitude. Subsequently polarization curves have been measured by 
linear sweep voltammetry, starting at open cell potential (OCP) with a scan rate of 0.5 mV s7!. 
Chemical characterizations of the ODCs are given in [23]. 


6.3 Analysis of measurement results 


In this section, the measurement results which will be used for validating the simulations are 
explained and interpreted. As discussed in detail by Franzen et al., the Ag ratio of ODCs has 
a major effect on the stationary performance of ODCs [23]. Thus, the polarization curves of 
four ODCS with distinct silver-PTFE ratio have been chosen for further analysis from previously 
published data [23]: 


* 92 wt—% Ag: This composition exhibits a medium performance and hydrophobic proper- 
ties (Electrode thickness: 292 um) [23]. 


e 97 wt—% Ag: The structure of an ODC with 97 wt—% Ag has been analysed by Neumann 
et al. [69], so that identified structure parameter can be used in this study (Electrode 
thickness: 297 um [23]). 


* 98 wt—% Ag: This composition exhibits the best performance [23] (Electrode thickness: 
335 um). 
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e 99 wt—% Ag: This composition exhibits a low performance and hydrophilic properties 
[23] (Electrode thickness: 356 um). 


In the following, the ODCs are denoted by their silver content; the PTFE content is equal to rest 
of the mass fraction. The polarization curves of these ODCs are given in fig. 6.1. The main 
differences between the electrodes occur at current densities above 2000 A m”?. 
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Figure 6.1: Experimental [23] and simulated polarization curves of ODCs with different Ag-content. 


To explore the reasons of the different performances, electrochemical impedance measurements 
are conducted in this study, the results of which are given in fig. 6.2 a). All impedance spectra 
show distorted semi-circles. For the ODC with 99 wt—%, a second semicircle emerges at high 
frequencies. The ODCs mainly differ in terms of ohmic and charge transfer resistance. The 
ohmic and the charge transfer resistance Ret of the different impedance spectra is determined as 
[62]: 

Ra = ZR(fmin) — ZR (fmax) (6.1) 


With Zp as the real part of the impedance |Z| and f as the frequency, where Zg (fmax) is the 
ohmic resistance. There is no systematic correlation between ohmic resistance and silver content 
or performance of ODCS; the differences may result from slightly different distances between the 
ODCs and the reference electrode, since sealing elements are tightened by hand and therefore 
their compressions differ. The resulting values for Ret are shown in fig. 6.2 b). The smallest 
charge transfer resistance is obtained for the ODC with 98 wt— 96, the largest for 99 wt— 96, which 
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Figure 6.2: a) Simulated and experimental EIS at j = 2000 Am”? of ODCs with different Ag-content, b) calculated 
charge transfer resistance and c) calculated time constants from experimental and simulated EIS. 


agrees with the performance ranking in stationary measurements. To identify the associated 
processes, the time constants were determined from the impedance data by eq. (6.2): 


1 


T= Z K 
27 f(—Zimax) 


(6.2) 
f(—Zimax) is the frequency where the negative imaginary part of the impedance reaches its 
maximum. The results are given in fig. 6.2 c). Two trends are visible: First, the time constant is 
increasing with an increasing Ag-content, but does not correlate with ODC performance. Second, 
for all four ODCs the characteristic time constant is in the range of 0.0015 s to 0.03 s, which 
agrees with previous results from potential step measurements [2]. It should be noted that the 
potential step measurements cover a larger range of time constants than EIS [35], and that no 
time constants slower than those found here were observed [2]. The charge transfer resistance 
reaches a minimum at 98 wt—% Ag, but the time constant decreases further until 92 wt—% Ag. 

In principle, the measured semi-circle and the corresponding time constants can be caused by 
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transport processes, by the interplay of reaction and double layer charging, or a combination of 
both. However, mass transport can be excluded in this case. The time constant of a diffusion 
process Tmt for species i in j can be approximated from Fick’s second law [13]. Equation (6.3) 
gives the diffusion time constant within a porous media, which is proportional to the square of 
the characteristic diffusion length, i.e. the electrode thickness Az, the porosity € and inversely 
proportional to the binary diffusion coefficient DE of species i in j. 


x —z € (6.3) 


The diffusion processes in ODCs were analyzed in a previous chapter. As demonstrated, a 
long mass transport way in the liquid phase leads to a distinct time constant in chronoamper- 
ometric measurements due to the water and hydroxide ion mass transport, whereas no distinct 
time constants for oxygen transport — in the gas as well as in the liquid phase — were found 
[2]. For j as water and i as NaOH, a flooding degree of the ODC from 5% to 50% with 
Di = 3.20 x 10719 m?s"! (cf. table 6.5) and € = 0.4 would lead to time constants of 0.28 
to 28.2s, which is orders of magnitudes higher than the measurement results. Also, the results in 
fig. 6.2 suggest that the observed time constant is not determined by oxygen transport in the gas 
phase, since the time constants become faster with increasing hydrophilicity and thus decreasing 
mass transport length in the gas phase. Although the current density of ODCs is limited by the 
availability of dissolved oxygen in the liquid electrolyte [1, 2, 32], oxygen transport in the liquid 
phase can also be excluded as the process causing the measured semi-circles for two reasons: 
First, as shown before, the diffusion length of the oxygen in the liquid phase is only several tens 
nanometers [32, 88] up to a few hundred nanometers [89], and the associated time constant is 
smaller than 10 x 1079 s [1]. This is orders of magnitudes faster than the time constants measured 
here. Second, the time constants of all ODCs differ while oxygen diffusivity in the electrolyte is 
not expected to be changed with the electrolyte distribution or electrode composition. 

In shutdown measurements with comparable ODCs, double layer capacitances of about Cg, = 
S0Fm”? were identified [72]. The time constant of the double layer charge / discharge taking 
place in parallel to a faradaic reaction is given by eq. (6.4) [62]: 


Ta = Cal -A geo Ret (6.4) 


With the characteristic time constant Ty, the double layer capacitance per geometrical area 
Ca, and the geometrical electrode area Ageo. Using Ca = 50Fm~, Ra = 0.2Q and Ageo = 
6.28 x 10? m? (cf. fig. 6.2 b), Ta = 0.003s is obtained, which is in the range of the measured 
time constants (cf. fig. 6.2 c). This hypothesis is also supported by the fact that the measured 
time constants correlate with the silver content — which increases the available surface area and 
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thus the double layer capacitance and decreases the time constant — and not with the ODC per- 
formance. 

In conclusion, the semi-circles of the EIS measurement represent only the double layer charge 
and discharge and no mass transfer. In this context, the distorted shape of the semi-circles is sur- 
prising because the interplay of double layer capacitance and a single reaction in a homogeneous 
electrode would lead to a non-distorted semicircle. Although two parallel reactions pathways for 
the alkaline ORR may take place, both have the same rate determining step [26, 30]. Thus, no 
influences on the dynamics of complex kinetic processes at the electrode surface are expected 
and the distortion of the semi-circles is most likely not caused by two overlaying reaction pro- 
cesses. The distortion of the semi-circles could, however, be caused by inhomogeneous electrode 
properties. As studied by Lukács, a distribution of process parameters within an electrochem- 
ical system leads to detached semi-circles [90]. Paulisch et al. measured an inhomogeneous 
electrolyte distribution in the ODC [80]. Based on these studies, I assume that the distorted 
semi-circles of the ODC impedances reflect locally different reaction rates and double layer ca- 
pacitances due to inhomogeneities of the electrolyte distribution. Gradients in the local reaction 
rates have inter alia already been reported for DMFC anodes and cathodes [91], or for acid ORR 
in PEMFC cathodes [92]. Based on this analysis of the experimental data, I employ a dynamic 
three-phase model which encompasses the effects of local inhomogeneities on the electrode to 
further analyze the processes in ODCs. 


6.4 Modelling 


Itis known, that the electrodes are inhomogeneous in particle size, pore size distribution as well 
as in the local catalyst and binder loading [69]. These local inhomogeneities in turn lead to a dis- 
tribution of hydrophilic and hydrophobic properties of the different pores [81] and may explain 
the previously observed electrolyte distribution [80]. The processes in the inhomogeneous three- 
dimensional GDE are very complex and include the effects of microkinetics [29], electrowetting 
[80] or closed pores [23]. Thus it is necessary to simplify the structure and processes in macro- 
scopic model approach to describe electrode performance. The presented model focuses on the 
effect of an inhomogeneous electrolyte distribution and three-phase area on electrode perfor- 
mance. Determining an electrolyte distribution from the electrode microstructure would require 
a detailed understanding of factors such as pore size distributions, local surface properties inside 
the pores or local tortuosity, but also two-phase flow in porous media, and is out of scope of this 
work. The electrolyte distribution is particularly important because it governs the mass transfer 
paths for the water and the hydroxide ions in the liquid phase, the oxygen dissolution via the 
size and location of the gas-liquid interface and the area for reaction and double layer charge / 
discharge via the size of the wetted catalyst surface per volume. 
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The model of an inhomogeneous ODC is developed by dividing the ODC in multiple domains 
perpendicular to the through-plane coordinate z. The domains are connected in parallel and dif- 
fer in the penetration depth of the electrolyte and the size of the gas-liquid interface, as shown in 
fig. 6.3 a). Each domain sub-model is discretised in through-plane direction (z), see fig. 6.3 b). 
This model is based on the assumption that the preferred mass transport is along the z coordinate 
in the through-plane direction, between the gas and the liquid bulk. Since the mass transport 
ways in in-plane direction are comparatively long (electrode thickness: 0.3 mm thick; electrode 
diameter: 20 mm) and run parallel to the two bulk phases, mass transport in in-plane direction is 
assumed to be negligible. The aim of the model is to demonstrate how the inhomogeneity influ- 
ences the performance and the dynamics of the ODC. The used approach does not allow one to 
replicate electrode structures in all aspects, but to derive important conclusions about the impact 
of the inhomogeneous electrolyte distribution. 

In the following, first the domain sub-model is summarized shortly. Subsequently, the intercon- 
nection of the domain sub-models to an inhomogeneous ODC is outlined. Model equations are 
given in table 6.1. 


6.4.4 The domain sub-model 


The domain sub-model is based on the expanded liquid mass transport model as introduced in 
chapter 5. With the help of the expanded liquid mass transport model, it was shown that — 
although sufficient water is available for the ORR - the dynamics and performance of the ODC 
are significantly affected by the mass transport of water and ions in the liquid phase. In order to 
consider these findings, following extensions have been made: 


* As demonstrated, even under unfavorable conditions enough water is available for the ORR 
and does not limit the reaction kinetics in any scenario. This indicates that the common 
assumption is valid for ODCs, which suggest that the formation of superoxide is the rate 
determining step for ORR in alkaline media at a silver catalyst (cf. eq. (2.11)) [24— 
29]. Under consideration of the discussion in section 3.3, the reaction rate of the ORR is 
implemented as 1-step kinetic depending on the rate determining step, see eq. (6.6). 


It has been shown that strongly non-ideal conditions exist in the liquid electrolyte due to 
high ion concentration. This may affect the mass transport in the liquid phase. To consider 
this, Maxwell-Stefan diffusion is modeled in dependence on the gradient in chemical po- 
tential. The chemical potential of NaOH is calculated respecting the law of electroneutral- 
ity using the mean activity coefficient y+, (see eq. (6.21)) [32]. Fick's diffusion coefficient 


of oxygen is adapted by an approach given by Bird et al. [93], (see eq. (6.23)). 
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Figure 6.3: a) Schematic view of the model of an inhomogeneous ODC comprised of multiple homogeneous sub- 
domains; b) schematic view of the domain sub-model including reaction and mass transport; c) zoom at the gas-liquid 
interface. 


* The calculation of the convective velocity within the liquid diffusion layer has been 
changed with regard to the model of an inhomogeneous ODC. It is discussed in detail 
in the next section. 


The changes compared to the models in chapter 4 and 5 are summarized in table 3.1. All imple- 
mented model equation of the domain-sub model are given in table 6.1. 
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Table 6.1: Model equations for domain-sub and inhomogeneous model. (Model schema for reference: fig. 6.3.) 


Electrode charge and reaction 


1 
Faraday current jg = —F-4- / r(z)dz (6.5) 
Ja-—zn 
: 1— a).F- i 
Reaction rate [29] r(z) = ko,cat * Scat * 40, * eXp ( ( x em (6.6) 
Reaction overpotential Nreaction = Eint — Eo (6.7) 
Lae 
Nernst potential [29] Eyed Ec In | 9% “#20 (6.8) 
4.F AN 
aOH 
Activity of O, ao, = XO, * YO, (6.9) 
b 
Activity of NaOH anaon = Ya (6.10) 
DNaOH 
n 
Molality of NaOH Buon eR. (6.11) 
Melectrolyte 
Ohmic drop Eext — Eint = Rspecific 2 (6.12) 
a dEint z a 
Charge balance ri Cdl.cat * Scat dz dr we F.4. T r(z)dz (6.13) 
Z| —Zfa ZI fa 
Mass transport in the gas phase and processes at the gas-liquid interface 
Opi 9? pj gas,eff dp; 
Species balance e e JZ "n + Jz Ves (6.14) 
with i € (H50,05] 
a © E 
Effective diffusion coefficient gu = gr “= (6.15) 
: d ux 
Bruggemann correlation :-g 0 (6.16) 
PO, 
Henry's law Co, = : (6.17) 
> Ho (chaom C850) 
Stefan flow xz ` Vgas = ae = ie (6.18) 
A ha S 1 fa 
Specific phase interface area Sphase,interface = Pase, ienee (6.19) 
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Mass transport in the thin film 


Species balance 


dci zy Di = u Cit} 
dt E RT dx Ox? / Ciotal 
ij 
(os e Cj Hn eed Ci 
Ox  Ox/ OX Cota Ox dx OX Ctotal 
| oc; 


F = ' Vliq,evaporation 
Oz g,evapı 


with i € {OH ,H,0,0,}, j € {OH ,H,O} (6.20) 
b 2 
[9] NaOH 
aOH = + + Hom = +RTIn 2 ( ) 
Chemical potential of NaOH POE ne ol S EINIGE 5 OF oH 
[32] (6.21) 
Chemical potential of H,O and W= ie +RT + In (a;) (6.22) 
O5 with i € {O5,H,O} 
ol 
Maxwell-Stefan Diffusity [93] ^ Do,,; = Zo, Ru (6.23) 
2 2 dIn(ao,) 
with je {OH ,H,O} 
| | nes 
Velocity due to evaporation Vliq,evaporation = e = (6.24) 
H,O 
Mass transport in the flooded agglomerates . 
n pilaeff 2 2 
dci D; j 9^uj O HjN cicj 
-e=y ( ) 
dt = RL dx Ox? / Cota 
¡Aj 
. res e) Cj ree a Ci 
Species balance Ox 0x7 OX Ctotal Ox ox OX Ctotal 
dc; 
er “Vig — T 
with i € {OH , H30, 03}, j € {OH ,H,O} (6.25) 
Z1 7 
Velocity due to reaction Viig,react (z) = (|vQRR = [vos ) f ud dz (6.26) 
H Z1—Zta total 
Total velocity of liquid phase Vliq (z) = Vliq,evaporation + Vliq,react (z) (6.27) 
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Mass transport in the liquid diffusion layer 


Ctotal 


dj y pie 9u Ou cic; 
i Ł (Fz 209 


ee ee cj ps a Ci 


Species balance Ox Ox OX Cotal Ox ax OX Ctotal 
^ dc; 
az Vidi 
with i € {OH , H50,0;], j € {OH ,H,0} (6.28) 
Interconnection of the ODC domains 
a a _ Adomain,k 
Specific domain area Sdomaink = — —— (6.29) 
Ageo 
: : =5 
Molar flow density out of / into ja" fa 
N* = ) Nik" Saomai (6.30) 
liquid diffusion layer L SER 
n 
. : f Diicic; (du; du; 
Molar flow density out of / into N Y BT - J ( = 2 J ) FCi* Vligk (6.31) 
flooded agglomerates a Ctotal ^ CX * c 
iz 2=2 
n=5 
Velocity viat = Y, Sdomain,k * Vig (Z =21) (6.32) 
k=1 
n=5 
Total current density Jopc = y Jdomain,k ‘ Sdomain,k (6.33) 
k=1 
Boundary conditions 
O; in gas bulk Po,(Z = 0) = Pia (6.34) 
H,O in gas bulk Pujo(z = 0) =0 (6.35) 
O; in electrolyte bulk co, (z= a) =0 (6.36) 
NaOH in electrolyte bulk CNaOH (Z = Zt) = Craon (Z = 2t,1 = 0) (6.37) 
H50 in electrolyte bulk“ colz = zt) = f(Cnaon(Z = zi. t = 0)) (6.38) 
Electrode geometry Zgas + Zif + Zfa = 21 (6.39) 


“ Calculated as in [51] 


6.4.2 The inhomogeneous ODC model 


This model connects the single sub-domain models to an inhomogeneously flooded ODC model. 
Each domain k takes up a specified proportion Saomain,k of the total geometrical electrode area 
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Aseo (eq. (6.29)). The domains are separated in through-plane direction, but connected by a 
common liquid diffusion layer (cf. fig. 6.3 a)). The concentrations c; at the boundary between 
the flooded agglomerates and the liquid diffusion layer are obtained by solving eq. (6.30). Ac- 
cordingly, the velocity of the convective flow in liquid diffusion layer is calculated by eq. (6.32). 
Since it has been shown that oxygen transport in the gas phase does not significantly affect the 
dynamics and performance of the ODC [1], it was decided not to interconnect the mass trans- 
port in the gas phase to simulate the inhomogeneities. Electrically, the domains are connected 
in parallel circuit, so that the total current density of the cathode equals the sum of all cathode 
domains, weighted by their share of the total area, eq. (6.33). 


6.4.3 Properties and parameters of oxygen depolarized 
cathodes and representative electrode domains 


In this section, the properties of the different ODCs and their representative domains are dis- 
cussed. The model parameters are classified on three levels: 


1. Parameters that are assumed to be identical for all ODCs and to have no distribution in 
ODCs: The catalyst surface area specific reaction rate constant Ko car and the catalyst sur- 
face area specific double layer capacitance Cal cat- 


2. Parameters that are assumed to be constant for each ODC, but differ between the individual 
ODCs: The ohmic resistance Rspecific, the porosity € and the specific surface area of the 
catalyst Scat. 


3. Parameters that are assumed to be distributed within the ODC: The intrusion length zfa of 
the electrolyte and the local size of the gas-liquid interface Sphase interface AS Well as relative 
size of the domains Saomaink- 


As it is clear that properties are distributed and that the goal of the study is to identify a distribu- 
tion that can reproduce the experiments via numerical parameterization, the number of simulated 
domains is limited to five. This number allows to have flooded, empty and several differently 
flooded pores; areas with similar flooding degrees are lumped together, as the domains do only 
interact via the joint bulk phases. Since some parameters are independent while other parameters 
are connected through physical interrelationships to others, their (in)dependencies are discussed 
in the following. 

First, an electrolyte distribution that covers a three-phase zone which extends across the entire 
ODC is fixed. For this purpose, the flooding length of the different domains zr, as well as the 
length of the three-phase areas AzapA — calculated by their left z3pA ; and right boundaries z3pa ; — 
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Table 6.2: Geometrical dimensions of the three-phase area of 
the different domains. 


Domain Z3PAJ^ Z3PA,r Zfa AZ3pA 

— um um um um 

1 0.01 zj 0-:z 0.005  z 0.01 - z 
2 0.032- zı O001l-q 002.3 0.022- z 
3 0.1:z 0.032 - zı 0.066 - zı 0.068 - z; 
4 0.316-z1 0.1-2 0.208 - zı 0.216: z 
5 1.0:z 0.316. z 0.658- zı 0.684- z 
a 


Logarithmic scale from 0.01 to 1: 


are considered: As illustrated in (cf. fig. 6.3 c)) the maximum intrusion depth of the electrolyte 
in each representative electrode domain marks the left side boundary of the three-phase area 
Z3pA.. It is assumed that some pores are blocked by hydrophobic PTFE at the beginning of the 
electrode, while other pores are completely flooded. Thus the maximium intrusion depth of the 
electrolyte is assumed to be in the rang between 1 % up to 100 % of the electrode thickness. A 
logarithmic scale was chosen since it is expected that due to the high hydrophobicity of PTFE 
many pores are impermeable to the liquid electrolyte, leaving a large part of the electrode area 
with almost no electrolyte intrusion. In table 6.2, the values of z3pa for the five representative 
electrode domains are shown. The flooding degree increases monotonically from domain 1 (min- 
imum flooding) to domain 5 (maximum flooding). A second assumption is that the three-phase 
area of the full electrode extends over the entire electrode thickness, but is split between the do- 
mains: For the more flooded domains, it is located closer to the gas bulk, and for the least flooded 
parts it is located close to the liquid diffusion layer. The right side boundary of the three-phase 
area z3pA ; Of one domain equals the left side boundary of the next less flooded domain so that 
there are no overlaps or gaps in the three-phase area. The thickness of the three-phase area is 
calculated by Az3pq = ZapA1 — Zapa,r. The mean length of the flooded agglomerates of a domain 
Zfa 18 the mean value of the left side and the right side boundary of the three-phase area of the 
domain (cf. fig. 6.3 a), more detailed in fig. 4.1 and fig. 4.2). The oxygen gradient within the 
liquid electrolyte is very steep, [32, 34], thus almost no ORR outside of the three-phase area 
takes place, especially at higher current densities [1]. The resulting values for zapA 1, z3pA .; and 
Zfa are given in table 6.2. 

Next the dependency of the local reaction rate (cf. eq. (6.6)) as well as the double layer 
capacitance (cf. eq. (6.13)) on the catalytically active Ag surface Seat area are discussed. It 
can be assumed that the surface area specific reaction rate constant ko,cat as well as the surface 
area specific double layer capacitance Calcat are constant material properties. Thus differences 
in total catalytic activity and double layer capacitance between the different ODCs only result 
from differences in electrode design, and thus in the specific catalyst surface area Scar and in 
the wetted share of the catalyst. Neumann et al. obtained an average specific surface area of 
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Scat = 6x 10° m? m~? Ag catalyst per electrode volume for the ODC with 97 wt—% Ag [69]. 
Taking these values into account, eq. (6.6) and eq. (6.13) can be used to identify ko cat and Cal,cat 
from the measurements of the ODC with 97 wt—%. The values for ko cat and Cal car are used for 
all investigated ODCs, and only the specific catalyst surface area is adjusted for each electrode. 
As shown before, the double layer capacitance depends on the current density [72], therefore, 
the value of Calcat determined here is only valid for current densities of j = 2000 A m2. 

Within ODCs the porosity € may be distributed. Future work might elucidate how the lo- 
cal porosity is related to the silver/PTFE distribution and the electrolyte intrusion. Since this 
complex relationship is not clear yet, I use the simplifying assumption of equally distributed 
porosities, so that average porosity € is set to be constant for all domains of a given electrode. 
The porosity, however, is dependent on the electrode design and electrode structure. It has been 
calculated to be 0.355 for the ODC with 97 wt—% Ag [69] and is adjusted for all further investi- 
gated electrodes and domains to reproduce the measurements. 

The area-specific internal resistance of the different electrodes Rspecific Was calculated from the 
impedance at the intersection with the real axis at high frequencies. The electric conductivity of 
the ODC is assumed to be high due to the high amount of Ag [32]; thus, variations in electric 
conductivity can be neglected and Rspecific 18 set to be constant for all domains. The proportion 
of electrode area attributed to a domain Sdomain and the specific gas-liquid interface Sphase interface 
of each domain have been adjusted individually to reproduce the measurements. 

The model was implemented in Matlab. Measured polarization curves (cf. fig. 6.1) as well as 
impedance spectra (cf. fig. 6.2) have been used for manual model parameter identification. The 
parameters were adjusted one by one in multiple iteration loops to minimize discrepancies be- 
tween experiments and simulation results. The identified parameters values for non-distributed 
electrode properties are given in table 6.3. The identified sizes and the specific gas-liquid in- 
terfaces for all ODC domains are given in table 6.4. All constant model parameter are given in 
table 6.5, all concentration-dependent parameters are given in table 6.6. The parameters and the 
model output are discussed in the next section. 
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6.5 Comparative model-based analysis of inhomogeneous oxygen depolarized cathodes 


Table 6.5: List of constant parameters and operating conditions. 


Name Symbol Unit Value 
Operating conditions 
Pressure in gas chamber“ P= po, Pa 1.01 x 10? 
NaOH concentration in electrolyte bulk“ CNaOH(Z = zi) molm ? 9.68 x 10? 
H50 concentration in electrolyte bulk [51] cg,o(z = zi) mol m? 50.0 x 10? 
Temperature^ D K 353.15 
Electrochemical and kinetic data 
Charge transfer coefficient? [24] a — 0.15 
Stoichiometric coefficient H,O VERE — —2 
Stoichiometric coefficient OH” vORR — +4 
Open circuit potential [23] E? V vs. 1.13 
RHE 
Binary diffusion coefficients 
For gaseous oxygen/water [73] HD ms! 2.95 x 10? 
For liquid oxygen/water“ [73, 74] Te ie ms! 2.18 x 107? 
For liquid oxygen/NaOH" [73, 74] do mH ms! 3.25 x 107° 
For liquid water/NaOH? [60] PA O NaoH ms! 1.26 x 107? 
Geometrical parameter o 
Thin-film thickness [32] Zte m 60 x 107? 
Thickness of liquid diffusion layer“ [2] Zldl m 10 x 1076 


^ Experimental operating condition. > Measured with 6.5 M NaOH. © Calculated with modified Wilke 
and Chang equation, with € — 3.9 and a — 0.5 chosen, adjusted to experimental data from [74]. 4 No 
NaOH data available, modeled for KOH with NaOH viscosity from [68], described in [32]. * Estimated 
to be independent of current density. Taken from my previous study where it was identified to be 10 um 
for measurements with vertically orientated ODC [2], the measurement results in the set-up used here 
correspond to those in the reference cell [23]. 


6.5 Comparative model-based analysis of 
inhomogeneous oxygen depolarized cathodes 


In fig. 6.2 a) the experimental and simulated impedance spectra are presented. The simulations 
reproduce the measurements well regarding the distorted shape of the semi-circles, the trend in 
charge transfer resistance with silver content (fig. 6.2 b)), and the time constant (fig. 6.2 c)). Ad- 
ditionally, in fig. 6.1, the experimental and simulated polarization curves are given, which also 
agree with slight deviations in the range of very high and very low current densities. Deviations 
in the low current range may result from the change of the Tafel-slop at a overpotential of about 
0.16 V vs. OCP (IR-corrected) [23, 32]. As discussed in [32], the reason for this behavior is not 
clear, possible explanations may be different adsorption isotherms of intermediates in low the 
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Table 6.6: Electrolyte concentration-dependent variables. 


Name Symbol Unit 
Inverse Henry constant“ [22] Ho, Pam? mol! 
Activity coefficient of 0, [50] Yo, — 

Activity of H,O“ [75] 44,0 — 

Mean activity coefficient of NaOH [66] Yr kgmol-! 
Density of electrolyte [68] p kgm ? 
Water vapor pressure [51] PLO Pa 


7 Extrapolated from cyaon < 6M, but as shown in [27] extrapolated values fit very well 
to experiments with cnaoHu < 12M, b Extrapolated from cyaoy < 5.5 M, but as shown in 
original research, calculations for e.g. ckog < 13.5 M are valid [50], ^ Extrapolated data 
from 343.15 to 353.15 K. 


potential region [94], or influences by electronic or ionic conductivities which also could lead 
to a change in Tafel-slope [95]. Even though the experimental setup is temperature-controlled 
[23], the electrode may heat up at high currents. This could explain the deviations in the higher 
current range. Moreover, it is possible that the inhomogeneities are not limited to an inhomoge- 
neous electrolyte distribution. This may lead to effects not considered in the model. The model 
is additional validated against EIS spectra of the high performing ODC with 98 wt—% Ag at cur- 
rent densities of 1kKAm”? and 4kA m7. As shown in the appendix D (fig. D.1) the simulation 
data agree well with the experimental impedance spectra in term of shape, time constant and 
charge transfer resistance. All in all, the presented model can describe the electrode behaviour 
reasonably well. 

In the following, first the influence of inhomogeneities on electrode performance will be dis- 
cussed by analysing the contributions of the representative electrode domains to the polarization 
curves. Subsequently, the influence of the inhomogeneities on the dynamic behaviour will be 
discussed by analysing the EIS spectra. The parameters of the individual ODCs in table 6.3 and 
6.4 are evaluated and compared in more detail during the analysis. 

In fig. 6.4, the individual contributions of the different domains towards the overall performance 
are shown for the ODCs with 92 and 98 wt—% Ag. 

The local current density of each domain is plotted versus the potential. In both ODCs the 
representative domains strongly differ in their performance. For the ODC with 92 wt—% Ag, 
the current density in domain 1, which has the lowest electrolyte intrusion, approaches a plateau 
at around 0.5 V, indicating a depletion of oxygen in the liquid phase. In contrast, the medium 
flooded domains 2 and 3, amounting to ca. 50 46 of total electrode area, are most active and have 
not yet reached their limiting current density at 0.5 V vs. RHE. The highly flooded domains 4 
and 5, amounting to ca. 40 % of electrode area, contribute only little to the total current density. 
Even though the ODC with 98 wt—% has comparable proportions of active regions, it reaches 
a significantly higher performance. To understand this discrepancy in performance between the 
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Figure 6.4: Simulated current density dependent on the potential in each domain k and for total ODC, for ODCs with a) 
92 wt—% Ag and b) 98 wt—% Ag. 


93 


6 In-depth Analysis of Electrolyte Distribution 


two ODCSs and their different domains, the differences of electrolyte intrusion as well as the gas- 
liquid interface have to be taken into account. In fig. 6.5 a), the proportion of the pore volume 
that is filled with electrolyte, py, is displayed over the dimensionless electrode thickness for all 
ODCs. The proportion is calculated according to fig. 6.3 a) by accumulating the relative sizes 
of the domains k, which is equal to the relative volume fraction of the respective domains, that 
are flooded at a respective location. It is assumed that pores are flooded to different degrees in 
the three-phase zone. For the calculation an average degree of 50 % flooding was used, since 
no further information on the degree of flooding in the three-phase area is available. Thus, only 
half of the relative size of the domain is taken into account if the three-phase area is present at 
the respective location. Eq. (6.40) describes the calculation of the electrolyte filling over the 
electrode depth, zg is the distance from the end of the electrode facing the liquid diffusion layer: 


1 
Li Sdomain,k —0.5- Sdomain,1 for za < 0.01 -z 


nee E? y Sdomaink — 0-5: Sdomain2 for 0.01 - zı < za < 0.032 zi 28 


iiem =05 Sims for 0.316- z1 < Za < z1 


Among the electrodes a trend of decreasing electrolyte saturation with higher PTFE content 
is clearly visible. This agrees with the results of a radiographic study of Paulisch et al. on 
similar electrodes, who found that a higher Ag content leads to a higher electrolyte saturation 
[80]. The effect can be attributed to the hydrophobic properties of the PTFE. According to the 
simulation results, also the porosity of the ODCs decreases with increasing PTFE content — this 
trend agrees qualitatively with FIB/SEM analyses of ODCs with varying Ag amount [23]. The 
porosity change may have an additional effect on the wettability of the electrodes as smaller pore 
diameters cause higher capillary forces. Although the wetted catalyst surface is the largest in the 
99 wt—% Ag, this ODC has the lowest performance. This is interpreted as an indication that the 
actual electrochemical reaction rate is not the performance limiting step, even at medium current 
densities. The poor solubility of oxygen in liquid electrolytes is a well-known phenomenon [22]. 
Therefore the size of the gas-liquid interface can limit the performance. With the simulation 
data, an estimation of the size of the interface dependent on the location within the ODCs has 
been made, which is given in fig. 6.5 b). The specific area of the gas-liquid interface within a 
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Figure 6.5: Identified properties of ODCs with different Ag-content: a) Electrolyte distribution dependent on the dimen- 
sionless location within the gas diffusion electrode (calculated by eq. (6.40)); b) average area of gas-liquid interface per 
electrode volume dependent on the dimensionless location within the gas diffusion electrode (calculated by eq. (6.41)); 
c) Schematic illustration of the electrolyte intrusion and gas-liquid interface distribution for reference. 
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given volume at each location zg is calculated by scaling the specific gas-liquid interfaces of the 
respective domains according to their share of the total electrode area (eq. 6.41): 
Sphase interface, 1 
AZ3PA,1 


Sphase interface,2 
Sphase interface? 6. an2 for 0.01 -z1 < za < 0.032- z 
Sy(za) = m ze (6.41) 


* Sdomain,1 for za < 0.01 -z1 


Sphase interface,5 


re -Saomains for 0.316- zı < za < z 


According to fig. 6.5 b), the ODC with 98wt—% Ag exhibits the largest gas-liquid interface 
and — as the reaction mostly takes place close to the phase interface — this electrode has thus the 
largest active three-phase region and the lowest charge transfer resistance. This causes the high 
current densities presented in fig. 6.4 b). It seems that the electrode composition of 98 wt—% Ag 
and 2 wt—% PTFE leads to an optimal balance between hydrophilic and hydrophobic properties 
and thus the highest performance. In contrast, the ODC with 99 wt—% Ag is highly flooded so 
that the low performance can be attributed to missing gas phase. On the other side, the ODC with 
92 wt— 96 Ag is barely flooded. Thus the gas-liquid interface is too small due to insufficient elec- 
trolyte intrusion. Furthermore, the model suggests what would be intuitively expected: Stronger 
hydrophobic properties of the electrodes lead to a shift of the three-phase boundaries towards 
the gas side. Almost the full gas-liquid interface, and thus active area, is located within the first 
10% of the electrode from the liquid side whereas much of the Ag electrode at the gas side is 
not electrochemically active. This also suggests that there is mostly a small distance only from 
the gas-liquid interface to the liquid bulk, and thus little mass transport resistance, which agrees 
with previous studies [2]. 

The claim of negligible mass transport limitation is also supported by the time constant analysis 
of the EIS measurements in section 2. As presented in a previous chapter, long mass transport 
ways in the liquid phase between the electrochemical active zone and the electrolyte reservoir 
lead to characteristic time constants [2]. Due to the fact that no mass transport time constants 
in any ODC were measured and due to the fact that the most flooded electrode with 99 wt— 96 
Ag performs worst, it is concluded that the flooded domains 4 and 5 of the ODCs contribute 
very little to the total electrochemical activity. The conclusion that the highly flooded domains 
are inactive is counter intuitive, since in this domains both the gas phase and the liquid phase are 
present, thus the presence of a physical gas-liquid interface is expected (cf. fig. 6.3). The electro- 
chemical measurements as well as the macroscopic model strongly suggest that the three-phase 
areas in the middle of the ODCs and close to the gas side are electrochemically almost "dead". 
To clarify the mechanism, further simulations on a pore scale or detailed experimental studies 
that link the electrode microstructure with the local performance would be necessary. Possible 
reasons for the low activity of largely flooded pores might be long mass transport ways across 
the ODCs or bottlenecks in the pore system leading to a high ionic resistance. Perhaps also local 
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salt accumulations due to the ORR [34] lead to a quasi insolubility of oxygen. The availability of 
oxygen affects the OCP and the reaction rate of the ORR. In fact, the here presented gas-liquid 
interface should be interpreted as the active interface of ionically and electrically well-connected 
pores with exposed silver surface, and not the physical gas-liquid interface. 

Now the impact of inhomogeneities on EIS is discussed. Using the presented model, the different 
local electrolyte saturation and the differences in local performance are shown to be responsible 
for the distorted shape of the semi-circles in the impedance spectra. The sinusoidal current of 
the total ODC is composed of the sinusoidal currents of all single domains. It represents the 
time constants and thus the phase shifts of all domains. The domains differ in double layer 
capacitance, in local performance, and thus in charge transfer resistance, therefore they all are 
characterized by their own time constant (cf. eq. (6.4)). The impedance spectra can be examined 
in more detail by including the Bode plots, which are given for the example of the ODC with 
98 wt—96 Ag in fig. 6.6 a) and b). Besides the phase shift @ and impedance |Z| of the entire 
ODC, the impedance values for the individual domains are also plotted separately. 

The minima in phase shift of the domains extend over four orders of magnitude in frequency 
as shown in fig. 6.6 a). The phase shift y at higher frequencies is governed by the less flooded 
domains (domains 1 and 2) and the more flooded domains (domain 4 and 5) govern the lower 
frequency range. The resulting overall phase shift curve reproduces the measurements well. 

The model also is able to reproduce the trend of the EIS time constants over the Ag content 
which is evident in the experimental results. In fig. 6.2 b), it can be seen that the characteristic 
time constant is slightly increasing with the silver content in experiment and in simulation. The 
characteristic time constant equals the product of the double layer capacitance and the charge 
transfer resistance (cf. eq. (6.4)). While Ra from 92 to 99wt—% Ag only increases by about 
a factor of 1.5, the time constant increases by a factor of 20 (cf. fig. 6.2). A major part of the 
increase in the time constant can be attributed to an increase in double layer capacitance: Due to 
the more hydrophilic properties of the ODCs with a higher Ag ratio, a higher part of the ODCs is 
flooded, which leads to a larger share of wetted catalyst surface, where the double layer charge 
and discharge takes place. Additionally, the higher amount of silver catalyst within the ODCs 
also results in a higher available total Ag surface area Scat. This also explains, why the time con- 
stant of the ODC with 92 wt—% Ag is faster than the time constant of the ODC with 98 wt— 96 
Ag although the charge transfer resistance increases. The identified double layer capacitances 
are within the range of previous studies [72]. 

As mentioned before, a slightly different shape of the EIS spectrum was observed for ODC with 
99 wt—% Ag: The impedance spectrum shows a second, less pronounced semi-circle. Accord- 
ing to the model results, the higher frequency semicircle can be assigned to the active pores, 
while the lower frequency semicircle is assigned to the flooded pores. In the remaining ODCs, 
the proportion of flooded pores is lower, so the impedance spectra of the flooded parts are less 
significant in the overall impedance spectra. 
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Figure 6.6: Bode plots for EIS at j = 2000A m? for ODC with 98 wt—% Ag for full electrode and single electrode 
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6.6 Model-based scenario analysis of impedance 
spectra 


As shown so far, the combination of dynamic and steady state electrochemical measurements and 
simulations enables to detect and quantify inhomogeneities and the related parameters in ODCs; 
the influence of electrolyte distribution within gas diffusion electrodes on their electrochemical 
performance and dynamic behaviour is clearly visible. A methodology which can identify and 
quantify inhomogeneities is a highly promising tool that enables knowledge-driven, targeted 
design of electrodes and state diagnosis. To get deeper insights into effects of the parameter 
distributions and to evaluate the impedance spectra and polarization curves for their informa- 
tive value and diagnostic capability regarding the electrolyte distribution, a scenario analysis is 
presented in the following: Two hypothetical alternative model scenarios are analyzed for the 
high-performing ODC with 98 wt—% Ag, the results are compared with the above presented 
experimentally parameterized model for an ODC with 98 wt—% Ag. In the first scenario, a 
completely homogeneous ODC with only one domain is modeled: Here no influences due to an 
inhomogeneous electrolyte distribution are expected. The second scenario is an ODC with five 
a 
in table 6.2; however, here equidistant domain areas and a constant specific gas-liquid interface 
= 6.4 um? um and Sh ain =20%. The gas- 
liquid interface has the same surface area per electrode volume at each position in z-direction 


domains with different degrees of flooding, where the flooding length z¥ is identical to that given 


phase interface 


S, 
over the whole electrode are assumed: AE 


within the ODC, but the electrode maintains its different degrees of flooding. The model param- 
eters for the two hypothetical scenarios are adjusted manually to best reproduce the polarization 
curve as well as the charge transfer resistance and time constant of the EIS. Model parameters 
and model characteristics are given in the appendix D. 

The outcome of the scenario analysis is summarized in fig. 6.7 where the simulated EIS spec- 
tra and polarization curves for all three scenarios are shown. All three model variations can 
reproduce the polarization curves, only minor differences between the model outputs are visi- 
ble. However, in case of impedances, the simulated inhomogeneous ODC model matches the 
EIS spectra most accurately. In contrast, the homogeneous ODC shows a single non-distorted 
semicircle, while the ODC with semi-homogeneous properties exhibits a significant smaller 
impedance over the frequency range up to 1 kHz, i.e. it has a lower charge transfer resistance. In 
addition, the model of a semi-homogeneous ODC shows signs of mass transport impacts, visible 
in the additional arc visible at low frequencies below 1 Hz. 

The underlying effects can be explained on the basis of the Bode diagrams. In fig. 6.6 c) and 
d) the Bode plot for the homogeneous ODC is illustrated. It can be seen that the minimum in 
phase shift is spread over approximately two orders of magnitudes in frequency. The phase shift 
curve is clearly not as wide as the measured phase shift, since a distribution of time constants is 
missing. The semi-homogeneous ODC features a phase shift over a wide frequency range (fig. 
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Figure 6.7: Experimental and different simulated scenarios: a) EIS at j = 2000 A m? for ODC with 98 wt—% Ag; b) 
polarization curve. 


6.6 e)). However, its phase shift as well as the impedance (fig. 6.6 f)) of the semi-homogeneous 
ODC diverge from the measurements at frequencies below 1 Hz, i.e. where mass transport ef- 
fects would be expected. It can be seen from the lower magnitude of the impedance |Z| that 
the highly flooded domains both 4 and 5 exhibits a small |Z| and thus contribute significantly to 
the total current. In these domains, the mass transport length of the liquid electrolyte within the 
ODC is so long that a time constant of the water and hydroxide ion mass transport at frequencies 
below 1 Hz can be observed. This results in overlapping mass transport semi-circles of the highly 
flooded domains 4 and 5 that form the peculiar almost horizontal line in the Nyquist diagram. 
Please note that due to the defined diffusion length no Warburg-type behavior is expected. As 
previously discussed, as a mass-transport related behavior is not visible in the experimental data, 
it can be concluded that such a semi-homogeneous condition is not present in the experimentally 
analyzed electrodes. 


100 


6.6 Model-based scenario analysis of impedance spectra 


a) 
—— Sim. inhomogeneous 
-500 F == Sim. semi-inhomogeneous 4 
----Sim. homogeneous 
-1000 
i 
E 
< -1500 
-2000 
-2500 } 1 
-3000 
5 10 15 20 25 30 
t/s 
50 ' —— Sim. inhomogeneous 
Sim. semi-inhomogeneous 
-1000r ----Sim. homogeneous J 
«v -1050F 
E 
<x 
-= -1100 
-1150 F 
10 10.5 11 11.5 12 
t/s 


Figure 6.8: Different simulated scenarios: a) Time-dependent current density response for a potential step from 0.8 to 
0.7 V vs. RHE at t = 10s and back at t = 20s; b) and c): Zooms of a). 


Methodologically, the question arises whether the information about the inhomogeneities in the 
electrolyte distribution is also accessible by means of potential steps as used in the chapters 4 and 
5. To evaluate this, potential steps from 0.8 to 0.7 V vs. RHE and back were simulated with the 
parameter sets of the three scenarios. The resulting time-dependent current densities are shown 
in figure 6.8 a), zooms are given in fig. 6.8 b) and 6.8 c). 

As shown in figure 6.8 a), all three scenarios have a positive current peak immediately after the 
positive potential step at 10 s and a negative peak immediately after the negative potential step at 
20 s, which are indicators for a fast process; for the semi-inhomogeneous a pronounced and for 
the inhomogeneous ODC a less pronounced slower relaxation follows the peaks, which are in- 
dicators for a slow process. The different current densities of the three scenarios after relaxation 
result from the different polarization behavior of the scenarios, cf. fig. 6.7 b). Similarities and 
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Figure 6.9: Time-dependent capacitiv current density response for a potential step from 0.8 to 0.7 V vs. RHE at t = 10s 
for full electrode and single electrode domains of the inhomogeneous ODC. 


differences in the dynamic behavior of the three scenarios are discussed in the following. 

The peaks are caused by the double layer charge and discharge as the fast process. Only with a 
strong zoom on the positive current peaks (fig. 6.8 b) and c)), differences between the double 
layer charge currents of the scenarios are recognizable. The peaks of the inhomogeneous and the 
semi-homogeneous ODC are widened compared to the peak of the homogeneous ODC, which is 
a result of the parameter distribution. To investigate this more deeply, the capacitive double layer 
charge current densities jq of the inhomogeneous ODC and its domains are given in fig. 6.9. 
As can be seen, the width of the double layer peaks and thus the time constant of the domains 
increases with their flooding degree. The double layer charge current of the full electrode is the 
sum of the charge currents of its domains, weighted by their share of the total area, cf. eq. (6.33). 
In the interval 10s < t < 10.001 s, a pronounced decrease in double layer charge current of the 
ODC can be noticed, which is attributable to the fast time constants of the less flooded domains 
1 to 3. But there is still charging at 10.001s < t < 10.01s, which is caused by the comparable 
long time constants of the strongly flooded domains 4 and 5. Concluding, the contributions of 
the strongly flooded domains 4 and 5 are mainly responsible for the widen shape of the double 
layer peak of semi-homogeneous and the inhomogeneous ODC. Due to the complex shape of the 
peaks, a determination of the double layer charge time constant based on the current response of 
the applied potential steps is not possible. 

After the double layer peaks, a pronounced, slower relaxation can be noticed for the semi- 
homogeneous ODC. The same effect is also weakly emerging for the inhomogeneous ODC, but 
not for the homogeneous ODC. The underlying dynamic process is the mass transport of water 
and hydroxide ions in the liquid phase between the reaction zone and the liquid bulk (see chap- 
ter 5). Two characteristics must be considered when comparing the relaxation behavior of the 
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current densities: On the one hand it is the time constant of the process, on the other hand it is 
the magnitude of the change in current density during the relaxation. The magnitudes of the re- 
laxation amounts about 180 Am”? for the semi-inhomogeneous ODC, but only about 10 Am”? 
for the inhomogeneous ODC. Accordingly, the mass transport of water and ions in the liquid 
phase has a large influence on the dynamics of the semi-inhomogeneous ODC, but only a very 
limited influence on the dynamics of the inhomogeneous ODC. As discussed in the impedance 
analysis, the more significant mass transport in the semi-inhomogeneous is caused by the larger 
contribution of domains 4 and 5 to the total current density. In the homogeneous ODC, the elec- 
trode dynamics is only marginally affected by the mass transport, therefore this process is not 
pronounced at all. Analogous to the double layer peak, it is not possible to determine the mass 
transport time constant for ODCs with distributed parameters by using potential steps. 

In summary, the simulations of the potential steps enable to reveal the double layer charging and 
the mass transport influences. Although the information about the distribution in double layer 
capacitance and charge transfer resistance is reflected in the current responses, they are difficult 
to identify, interpret and especially to quantify. As discussed previously, this information is ac- 
cessible with impedance analysis. Thus, EIS seems to be a more suitable tool for the analysis of 
these processes. For the inhomogeneous ODC, a minimal relaxation due to the mass transport 
of water and ions in the liquid phase was detected, which remained hidden in the EIS spectra. It 
seems that the mass transport is slightly better resolved by potential steps than by EIS, although 
the time constant of this process cannot be determined by the potential step analysis. However, 
since the current changes by only about 1 96 during relaxation, the potential steps indicate that 
the mass transport in the liquid phase is of minor importance for the dynamics of the inhomoge- 
neous ODC. 


Concluding, although all scenarios match the steady state data, only the inhomogeneous model 
can reproduce the impedance spectra. The information about the inhomogeneities is difficult 
to access by means of potential steps. In the impedance spectra, the homogeneous electrolyte 
distribution does not meet the distorted semi-circles; the semi-homogeneous distribution of the 
electrolyte exhibits an additional time constant due to the active regions close the end of the 
ODC facing the gas side. The ability of EIS for detecting time constant and their distributions 
implies that it is a suitable tool for investigating and quantifying the gas-liquid distribution in gas 
diffusion electrodes and its impact on the electrode performance. 


6.7 Conclusions 


In this chapter I investigated the influence of electrolyte distribution for ODCs with different 
binder-silver ratios on electrode performance and electrode dynamics using experimental as well 
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as modeled polarization curves and impedance spectra. 

Experimental impedance spectra of oxygen depolarized cathodes with different Ag- and PTFE- 
content have been presented which were mainly characterized by fast time constants of 7 < 0.03 s 
and distorted semi-circles in the Nyquist plot. For a deeper understanding of the impedance 
spectra, the first model of an ODC with distributed properties was presented. The fast time 
constants can be assigned to the double layer charge and discharge, whereas the distorted semi- 
circles are a result of locally varying electrolyte intrusion depth. The model results suggest that 
pore systems with a low electrolyte intrusion make a significantly higher contribution to the total 
performance than the highly flooded pores. A comparative model-based analysis of ODCs with 
different silver to PTFE ratios has shown that the electrodes differ significantly in performance 
and dynamics. Performance is not attributed to a large wetted catalyst surface but to a large active 
gas-liquid interface that is located close to the electrolyte bulk. Furthermore, the model results 
were compared with two hypothetical scenarios of a homogeneous ODC and an ODC which is 
active over the entire electrode thickness. It was shown that the homogeneous ODC exhibits 
undistorted semi-circles, whereas the electrode with the semi-homogeneous gas-liquid interface 
shows mass transport effects, which are both not visible in the experimental impedance spectra. 
The presented insight will aid experimentalists and electrode designers to better diagnose the 
state of their electrode using EIS and to judge how heterogeneous or distributed the electrolyte 
phase is. EIS will also help to monitor electrodes: As electrodes age, their hydrophobicity might 
change and thus electrolyte distribution, EIS and performance. Future studies are needed to 
investigate the local performance and the structure on a pore-scale level and to establish a direct 
link between the microstructure and ODC performance. 
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In this chapter, the dissertation is summarized and key findings are highlighted. Additionally, 
perspectives and open research questions regarding ODCs are discussed. 


7.1 Summary and conclusions 


The goal of this dissertation is to gain in-depth understanding of the processes in ODCs for 
advanced chlor-alkali electrolysis in order to identify limitations and to provide the fundamen- 
tals for a knowledge-based electrode optimization. To achieve this goal, dynamic three-phase 
models were developed, parameterized and systematically improved based on simulation results 
and current literature findings. The novelties of the models are in particular (1) the dynamic 
approach, (ii) the appropriate consideration of the water activity and water mass transport, (iii) 
the inclusion of the phase equilibrium dependent on the ion and water concentration as well as 
(iv) the representation of structural inhomogeneities of the porous ODCs. 


Initially, a basic model of an ODC was developed to investigate the processes and limitations 
in ODCs. With the help of the basic model it was demonstrated that under industrial operating 
conditions ORR only takes place close to the gas-liquid interface due to a lack of available oxy- 
gen in the liquid phase. A key factor for the low availability of oxygen in the liquid electrolyte 
is its poor solubility as a result of high ion concentration and low water activity. Since during 
ORR water is consumed and hydroxide ions continue to accumulate near the gas-liquid interface, 
the oxygen solubility decreases further with increasing current density. This effect contributes 
significantly to the total depletion of oxygen within the liquid electrolyte at high reaction rates 
and thus to a limitation in current density. A dynamic analysis was performed, where the mass 
transport of the water and hydroxide ions in the liquid phase was identified as the slowest process 
in the ODC. Since water and ion concentrations affect the oxygen solubility as well as the water 
evaporation, the dynamics of concentration or partial pressure profiles of all species follow this 
process. 

For an in-depth analysis of water and ion mass transport in the liquid phase and its influence on 
the ODC performance, the expanded liquid mass transport model was developed. With the help 
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of this model, electrochemical measurements of two Bee > differing regarding the water and 
ion mass transport in the liquid phase, were compared. The model-based analysis of dynamic 
and steady-state measurements has shown that a forced convective mass transport in the liquid 
phase accelerates the dynamics and increases the performance of an ODC compared to an ODC 
without a forced convection within the liquid electrolyte. The model suggested that the forced 
convective flow reduces the thickness of the liquid diffusion layer by a magnitude, which leads to 
faster the mass transport in the liquid phase. As a result the accumulation of hydroxide ions and 
thus the decreasing water activity and oxygen solubility is less pronounced. As demonstrated, 
water activity does hardly influence the kinetics of the ORR directly, but rather indirectly via the 
oxygen solubility. 

Next, it was investigated how the structure of an ODC and an inhomogeneous electrolyte distri- 
bution within the ODC affect the water and ion mass transport, the local oxygen availability and 
finally the electrode performance. For this purpose, the inhomogeneous model was developed. 


Using the inhomogeneous model, polarization curves and EIS ia of ODCs with 
different binder-catalyst ratios were investigated. It was demonstrated that the EIS measurements 
of all ODCs are characterized by fast time constants of T < 0.03 s and distorted semi-circles in 
the Nyquist plot. While the fast time constants can be assigned to the double layer charging 
and discharging, the distorted semi-circles seem to be a result of inhomogeneous electrolyte 
distributions within the ODCs. The model results suggested, that electrode areas with a high 
electrolyte saturation are almost electrochemically inactive and barely contribute to the electrode 
performance. The dynamics and performance of the ODCs with different binder-catalyst ratio 
were compared with aid of the inhomogeneous model. It was shown that not the area of the 
wetted catalyst, but the size and location of the gas-liquid interface within the ODCs are the key 
driver to a high electrode performance. 


Concluding this work, the understanding of processes and limitations within ODCs as well as 
their dependency on the environmental conditions and electrode structure was expanded fun- 
damentally. These insights are relevant for both academia and industry. It was shown that an 
enrichment of ions at the gas-liquid interface and the corresponding decrease in water activity 
limits oxygen availability in the liquid phase and thus ODC performance. As concluded, the 
performance of technical ODCs can be optimized by a higher mass transfer rate in the liquid 
phase, for example by a convective electrolyte flow or an optimization of the electrolyte flow 
field. Also, new requirements for high-performance ODCs can be formulated on the basis of 


Measurements performed by Institute of Chemical and Electrochemical Process Engineering, Technische Universitit 
Clausthal, Germany and Analytical Chemistry — Center for Electrochemical Sciences, Ruhr University Bochum, 
Germany, published in [2]. 

Measurements performed by Institute of Chemical and Electrochemical Process Engineering, Technische Universität 
Clausthal, Germany, published in [15]. 
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the investigations. The structure of the ODC must not only support a large size of the gas-liquid 
interface, but also a location of the gas-liquid interface close to the liquid side of the ODC. 


From a methodological point of view, the use of dynamic methods is indispensable, as this en- 
ables the processes to be separated from each other by their time constants and to determine their 
specific influence on electrode performance. Stationary methods are limited in this matter. For 
an investigation of the interaction between different processes and their effect on the electrode 
performance, a one-dimensional approach has been shown to be sufficient, since the dynamics 
of the different processes as well as their interaction can already be represented in this manner. 
However, this approach does not allow to represent inhomogeneities or parameter distributions; 
this is enabled through parallel interconnection of electrodes with different properties. 

As was shown, for ODCs potential steps are particularly suitable for identification and analysis 
of processes with slower time constants as mass transport in the liquid phase. Impedance spectra 
are appropriate for a more detailed status analysis of the complex three-phase system due to their 
high information content with regard to electrolyte distribution or local reaction rates. 


7.2 Open questions and perspective 


In this dissertation, ODCs were analyzed by using macroscopic models. It was possible to iden- 
tify performance drivers as well as limiting factors of ODCs. However, the results have also 
raised new questions, which can be answered in future studies by the utilization and further ex- 
tension of the developed models: 

For the ODC in the measurement setup with the forced convective electrolyte flow, a compara- 
tively small thicknesses of the liquid diffusion layer was identified. A reason for this phenomenon 
may lay in additional mass transport effects, which are not yet covered by the model. There may 
be, for example, an additional natural convective flows perpendicular to the electrode due to 
density gradients as a result of concentration or temperature gradients. For the clarification of 
this aspect, e.g. CFD simulations would be a suitable diagnostic tool, whereby the macroscopic 
model can provide necessary data as ion concentration profiles for different current densities and 
states of the ODC. 

Since the herein analyzed measurements were performed in temperature controlled setups, the 
power-related heating of the system was not considered. However, even at low current densities 
of j = 3000A m? at an overpotential of y = 0.8V (cf. fig. 5.2), the system is heated with a 
thermal power of P = 2400 W m ?. Therefore, it is conceivable that the temperature is not ideally 
constant in terms of time and space. In the potential step measurements under certain conditions, 
a time-dependent behavior of the current was observed, which is opposite to the effect of water 
mass transport. Also divergences between the measurements and simulations at high current 
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densities were obtained. Both may be explained by temperature influences. To study the impact 
of this effect on electrode performance, the macroscopic model can be expanded to cover the 
temperature influences. For model validation combining electrochemical measurements with 
local temperature measurements would be suitable. 

The macroscopic findings indicate that there are strong gradients in local performance within 
an ODC. To investigate the underlying mechanism and its relation to the ODC structure more 
deeply, simulations on the pore scale would be useful, taking effects as bottle necks and dead 
spaces within the pore network into account. For this purpose, the macroscopic ODC model 
could be transformed into a pore network model. 
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A Maxwell-Stefan Diffusion 


The stationary Maxwell-Stefan diffusion 

The following exposition is mainly based on Multicomponent mass transfer by Taylor and 
Krischna [59] and Mass transfer in concentrated binary electrolytes by Newman et al. [60]. 
The execution is constrained on the one-dimensional case. 

To understand the Maxwell-Stefan relations and to adapt this approach to the specific system of 
the ODCs, first the mechanism of molecular collisions have to be understood. The principle of 
molecular interaction is initially presented simplified for a binary system and generalized after- 
wards. The molecules of the species 1 (e.g. water) and 2 (e.g. hydroxide ions) move with an 
average velocity of vı and v2. Introducing further the mass of a single molecule as m, and mo», 
the momentum of the average molecule of both species can be written as vı -mı and vz - m». In 
case of a bilateral collision of molecules of both species, the velocity of both molecules changes, 
where vi and v are the velocities after the collision. Thereby the total momentum v; - m +v2: m» 
is conserved [59]: 


mi: (v4 —V4) +m: (v2 —v5) =0 (A.1) 


From eq. A.l, the momentum transferred form the molecule of species 1 to the molecule of 
species 2 (and the other way round) can be calculated by [59]: 


my - (vı — v1) = ma: (v2 — v4) (A.2) 


For determining the velocity of the molecules after collision, an ineleastic collision is assumed, 
where the molecule 1 and 2 stick together and move with the velocity v} = v^ = v’, where v’ is 
calculated by [59]: 


E mM] -V + M2 :V2 


(A.3) 


mı +m 
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From eq. (A.3) the calculation of the exchanged momentum of the collision of two single 
molecules eq. (A.2) can be enlarged to eq. (A.4). Now it is shown in eq. (A.4), that the ex- 
changed momentum is proportional to the difference of the velocity of both species (vi — v2) 
[59]. 


m, ma 
mi (vı — v4) = m: (v3) = 
mı +m 


(vi v2) (A.4) 


Three important points can be extracted from this explanation so far [59]: 
* When two molecules of different species collide, they exchange momentum. 


* Although the momentum is conserved, the total kinetic energy is lowered and mainly con- 
verted to heat [96]: 


1 1 
jm vt +m. v5) > jm +m): v 2 (A.5) 


It should be noted, however, that the kinetic energy of one of both species could increase. 
This implies that both positive and negative accelerations of molecules can be achieved 
through interaction. 


* When two molecules of the same species meet, there is no loss of momentum for this 
species [59]. Accordingly, the momentum exchange between the molecules of a species 
does not have to be taken into account. 


To transfer these considerations from single molecules to a macroscopic level, the frequency of 
collisions has to be respected: The number of collisions of molecules of the species 1 and 2 per 
unit and time is proportional to the product of their mole fractions x; -x2, with x1 = c1 /Ctota and 
X2 = 2/Ctotal, Where c is the concentration and Crotal = c1 + c2 is to total concentration [59]. 

The two approaches can be combined: On a macroscopic level, the exchanged momentum 
Apm,ı,2 between species 1 and species 2 is proportional to the exchanged momentum at the 
collision of two single molecules and the frequency of collision [59]: 


APm,1,2 = (Gi -x2), (v1 — v2)) (A.6) 


Introducing now the Maxwell-Stefan diffusion coefficient D 2 as the inverse drag coefficient, the 
quantity dı 2 is obtained [59]. Where dı > - RT : ctota is the force per volume which species 2 is 
applying on species 1 resulting from the relative motion of species 1 and species 2 [60]: 

X (v1 — v2) 


dj 5 = ————À————— A.7 
12 Bua (A.7) 


130 


A Maxwell-Stefan Diffusion 


A short excursion to Newton's second law shows that the force species 2 is applying on species 
l, is the same force as species 1 is applying on species 2 [59]. 


di? = —d»41 (A.8) 
where 
D¡2=D21 (A.9) 


To generalize equation (A.7) from a binary to a multi-component mixture with n species, the 
forces affecting species i di- RT - Ctota¡ can be calculated as the sum of all binary interactions 
[59]: 


n n Ye cms 

d-Yyd,--y (vi vi) (A.10) 
j=l j=l Dij 
¡Aj Aj 


So far, the force which is applied on the molecules of one species i due to their relative velocity 
to the molecules of other species j has been described. Additionally, due to their gradient in the 
chemical potential, a drag force acts on the molecules of species i. This drag force is defined for 
non-ideal systems as d; RT - cii [60], with [59]: 


| Xi du; 
! RT dz 


(A.11) 


It can be recognized that the drag force of the species i is proportional to its mole fraction x; and 


its chemical potential gradient Si [60]. 
Finally the right side of eq. (A.11) can be substituted in the left side of eq. (A.10): 


cj du; E = cjci (vj — vi) 
RT dz | 


(A.12) 


j=l Ctotal Di, 


itj 
The eq. (A.12) is better known as the (stationary) Maxwell-Stefan diffusion. As illustrated in 
fig. 3.4, in the Maxwell-Stefan diffusion the drag force of species i is equal to the sum of the 
interaction forces between the species i and the species j, where i % j [59, 60]. 


Idealization of the Maxwell-Stefan Diffusion 
A simplification of the equation (3.20) can be achieved if an ideal behavior of the involved species 
is assumed. In that case, the concentration can be used as a driving force instead of the chemical 
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potential. The calculation of the chemical potential is given in eq. (A.13), with u? as chemical 
potential under standard conditions and a; as the chemical activity of species i [97]. 


Hi = We +RT -In(a;) (A.13) 


Whereas the chemical activity of species i is calculated by its concentration cj, its concentration 
under standard conditions ra and its activity coefficient y, [97]: 


itu (A.14) 


With eq. (A.14), eq. (A.13) can be rewritten as: 
Ui = UÊ +RT -In(c;) - RT -In(c®) + RT -In(y) (A.15) 


Please remember, the goal of simplification is to replace the gradient in chemical potential by the 
gradient in concentration as the driving force. Accordingly, the quantity sought is the derivative 
of the chemical potential over the location. This can be gained by the total differential of eq. 
(A.15) at a constant temperature: 


9 i 3 i 
OM; IUP gnc) pp Ane?) | pp 2In(y) (A.16) 


Oz Oz Oz Oz Oz 


Since ue and e = | molL-! are constant and % is defined to be 1 in an ideal system, eq. (A.16) 
can be reduced to: 


OH ol 
Oz E o 


= (A.17) 


n Ci) = RT dc; 
zZ Ci Oz 


By substituting eq. (A.17) into eq. (3.17), the flux N; is now expressed in dependence on the 
gradients in concentration of the species i and j: 


n sus O 
s-- Eo tS (A.18) 


Oz c total dz c total 


With a transformation analog to eq. (3.18) to eq. (3.20) finally the simplified Max well-Stefan 
equation is obtained, which regards the interaction between the species i and j and the gradient 
in concentration as the driving force. All further before mentioned assumptions are still valid: 


de; _ ; eS Fj Fa 7 (A.19) 


de Ctotal Je Ctotal 
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For all concentration dependent calculations the assumption of a constant total concentration 
Ctotal = Ctotal(Z = 0,1 = 0) is used. As shown in figure B.1 the respective error is less than 0.8 %, 
1.e.: 


Ctotal (CNaOH) = Ctotal (CNaOH (z HZ, t= 0)) 
Ctotal (CNaoH (Z = Zt, t= 0)) 


< 0.008, 


for cNaoH = [1.05 x 10*...1.75 x 10 ^ molm?| (B.1) 
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Figure B.1: Total, water, and sodium hydroxide concentration of liquid electrolyte, calculated based on Hirschberg 
[51] and implemented in the model, depended on the modeled water concentration in the range of the simulation at 
T =353.15K. 


13 Parts of this chapter have been published in RGhe et al., ChemSusChem, 2019 [1]. 
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Figure B.2: Inverse Henry constant over water concentration at T = 353.15K, hydroxide concentration according to 
figure B.1. 


o 


vap 
p. o / Pa 
2 
Da ME 
N A 


E 

o N 
r r 
f , 


16} 4 


Water vapor pressuer p 


0.6 1 1 ï 1 1 1 
43 44 45 46 47 48 49 5 
Water concentration c,, o / mol.m? x10* 
2 


Figure B.3: Vapor pressure of water as a function of water concentration at T = 353.15 K, hydroxide concentration 
according to figure B.1. 


As seen in figure B.2 the inverse Henry constant is highly sensitive to the water respectively the 
hydroxide concentration and increases over the simulation range by a factor of 5.81 [22]. 
Dependent on the water and hydroxide concentration the vapor pressure of water changes de- 
creases over the simulation range by a factor of 2.91 (see fig. B.3) [51]. 

As shown in figure B.4, based on the enrichment of hydroxide ions the water activity decreases 
by a factor of 2.22 [75]. 
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Figure B.4: Activity of water as a function of water concentration at T = 353.15 K, hydroxide concentration according 
to figure B.1. 
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Figure B.5: Density of electrolyte as a function of water concentration at T — 353.15 K, hydroxide concentration ac- 
cording to figure B.1. 


The concentration change of the liquid electrolyte leads to a volume change, and thus to a con- 
vective flow (cf. B.5) [68]. 
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C.1 ODC preparation 


The preparation and characterization have been done by Barbara Ellendorff and David Franzen, 
Clausthal University of Technology, Germany. Section partly written by David Franzen. 

The investigated ODCs were fabricated by means of a wet spraying technique [9]. Silver particles 
(SF9ED Ag Flake, Ames Advanced Materials Corporation, x59 = 1.61 um), a methyl cellulose 
solution (1 wt-%, WALOCEL™ MKX 70000 PP 01) and additional demineralized water were 
mixed twice with an Ultra Turrax at 13500 rpm for five minutes until a homogeneous dispersion 
was obtained. Under a lower rotational speed of 1000 rpm, to prevent agglomeration, a PTFE 
dispersion (TF 5060GZ, 3M™ Dyneon™) was added. Up to 90 layers of the suspension were 
applied on a nickel mesh as a conductive supporting material (106 um x 118 um mesh size, 63 um 
thickness, Haver & Boecker) by using a spraying piston (Evolution, 0.6 mm pin hole, Harder 
& Steenbeck). The mechanical stability was improved in the following production steps by hot 
pressing (LaboPress P200S, Vogt, 15 MPa, 130?C, 5 min) and heat treatment (330 °C, 15 min, air 
atmosphere), burning out the methylcellulose and sintering the PTFE. This results in a mechani- 
cally stable electrode with a pore system that meets the requirements of a gas diffusion electrode. 
The electrode thickness and the catalyst load were determined using a thickness dial gauge (FD 
50, Käfer) and weighing respectively. 


C.2 Experimental set-up with electrolyte convection 


The measurements have been done by David Franzen, Clausthal University of Technology, Ger- 
many. Section partly written by David Franzen. 

Electrochemical experiments were carried out in an in-house built half-cell with a flow-through 
electrolyte and gas chamber separated by the ODC. Measurements were recorded using a Gamry 
Reference 3000 potentiostat (Gamry Instruments, United Kingdom) with a three-electrode con- 
figuration. The ODC was operated as the working electrode, while a platinum wire was placed 


14 Parts of this chapter have been published in Róhe et al., ChemElectroChem, 2019 [2]. 
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downstream in the electrolyte chamber as the counter electrode. Potentials were measured with 
respect to a reversible hydrogen electrode (RHE) (Hydroflex, Gaskatel, Germany) through a Lug- 
gin capillary ending in front of the working electrode. The electrolyte, 30 wt-% NaOH prepared 
from caustic flakes (> 99 wt-%, Carl Roth, Germany) and demineralized water, was pumped 
through a heating tube and the half-cell from a heated reservoir providing a constant temperature 
of 50°C. On the gas side pure oxygen was supplied with a flow rate of 50 mLy min”!. The 
pressure on the electrolyte and gas side was set to 1.00 x 10? Pa. Sensors on the electrolyte and 
gas side monitored pressure and temperature. After a startup procedure (—2kAm”? hold for 
10 min, —A4kA m ? hold for 50min, —2kA m"? hold for 5 min), potential steps of 100 mV in the 
range from OCP (1.05 V vs. RHE) to 0.10 V vs. RHE were applied and the current response was 
recorded. 

After 30 s, current density values were extracted to record the polarization curve. For dynamic 
measurements the potential of 0.8 V vs. RHE was held for 30s and stepped to 0.7 V vs. RHE. 
After a further 30s, the step was reversed in the opposite direction. 


C.3 Experimental set-up with SECM and stagnant 
electrolyte 


The measurements have been done by Alexander Botz, Ruhr Universität Bochum, Germany. Sec- 
tion partly written by Alexander Botz. 

For SECM the ODC was measured in an in-house built half-cell with a stagnant electrolyte above 
the electrode and the oxygen flow from below. The whole setup was built on an active damped 
table (Newport RS 2000, United States) to reduce noise. An in-house built faraday cage was 
equipped with isolating vacuumed polystyrene panels (Vaku-Isotherm, Germany) to ensure a 
constant temperature within the time scale of the experiments. The heat-conductive base plate 
was connected to a thermostat (Huber CC2, Germany) for the direct heating of the electrolyte to 
the operating temperature 50 °C. The oxygen back pressure (> 99.5 %, Air Liquide, France) was 
regulated via a pressure controller (MKS Instruments, United States). All electronical compo- 
nents and the operator (using a foot connector; Conrad, Germany) were connetcted to the same 
ground, preventing ground loops and damages at the tip's apex [98]. 

All electrochemical measurements were conducted utilizing an analogue bi-potentiostat (IPS PG 
100, Germany) and the potentials were measured against a commercial RHE reference electrode 
(Mini Hydroflex, Gaskatel, Germany). A Pt mesh counter electrode was placed behind a thin 
D5 ceramic frit to prevent platinum species to diffuse to the working electrodes. All NaOH (> 
98.9 wt-%, J.T. Baker, The Netherlands) solutions were prepared with ultra-pure water (MilliQ, 
SG, Germany). In the measurements all solutions were air-saturated and heated up to 50°C in a 
water bath to the required temperature before being filled into the pre-heated measurement cell. 
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To ensure stable conditions at the sample electrode a conditioning procedure (cyclic voltammo- 
gramms, potential range between 0 to —1 V vs. RHE, 10 cycles, 100 mVs ^!) was performed. 
After ODC conditioning the solution inside the cell was replaced to avoid an activity change of 
the electrolyte prior to the measurement. During the measurements the backside of the ODC was 
in contact with oxygen at 1.013 x 10? Pa, allowing passive oxygen transport. First the OCP of 
the ODC was determined as 1.1 V vs. RHE. Then the ODC was polarized step-wise between 
1.1 V and 0.2 V vs. RHE. After the current at the initial potential reached the steady state, the 
potential was stepped 0.1 V in cathodic direction and the current was recorded for 100s. From 
this measurement sequence the polarization and the current dependent NaOH activity were ob- 
tained. For dynamic analysis, an additional potential step from 0.8 V to 0.7 V and back to 0.8 V 
vs. RHE was recorded. Prior to the first step a steady state current was reached. The individual 
potential steps were recorded for 60s. 


C.4 SECM electrode fabrication and measuring 


As tip electrodes for the activity measurements platinum ultra micro electrodes (UME) were em- 
ployed, which were prepared according to a previously published procedure [99]. Briefly, a laser 
puller (Sutter Instruments P-2000, United States) was utilized to pull Pt wires (9: 25 um, Good- 
fellow, Germany) inside of quartz glass capillaries (Sour: 0.9 mm, Øin: 0.3 mm, E: 100 mm, QSIL, 
Germany) by sealing them in 10 to 12 cycles (20 s laser on and 405 laser off) with the following 
parameters: Heat: 800, Filament: 5, Velocity: 128, Delay: 130, Pull: 0. Directly afterwards they 
were pulled: Heat: 820, Filament: 5, Velocity: 128, Delay: 130, Pull: 160. After soldering the Pt 
wire with soldering tin (Conrad, Germany) to a copper wire, the electrode tip was polished with 
an in-house built machine that rotates polishing paper and the electrode individually as described 
previously [99]. The electrochemical behavior as well as the active electrode size were subse- 
quently determined in a solution of 5 mM [Ru(NH3)g]Cls (> 98 wt-%, Sigma Aldrich, Germany) 
in 100 mM KCI (100 wt-%, J.T. Baker, The Netherlands). 

Pt UMEs were approached to the ODC by means of shear force distance-controlled scanning 
electrochemical microscopy (SECM). This technique is independent of the current signal and 
detects hydrodynamic forces in close proximity to the substrate surface [100, 101]. Therefore, 
two piezo elements (Piezomechanik Pickelmann, Germany) were mounted to the electrode body 
at an angle of approx. 45°, ca. 1 to 1.5cm distance from each other and individually connected 
to a lock-in amplifier (Ametek 7280, Germany). While the upper piezo oscillates the tip apex the 
lower one detects the resonance movement of the tip. Characteristic resonance frequencies were 
determined by comparison of frequency spectra (range: 200 to 500 kHz) in the bulk solution and 
within the shear force interaction distance (~ 200 nm above the surface). Pre-positioning of the 
tip at different lateral positions above the ODC surface was controlled using a video microscope 
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(The ImagingSource USB-camera, Germany). A simplification of the previously reported 4D 
shear force-based constant-distance mode SECM [102] was used for precise positioning of the 
SECM tip, utilizing an analogue controlled piezo-based positioning system (PI, Germany). An 
automatic shear force-based approach curve was performed until a stop criterion of 5 96 of the 
lock-in magnitude value was reached. For all experiments, the tip (9: ~ 1 um) was subsequently 
retracted to 1 um to reach the final working distance for the cyclic voltammetry operation. The 
cyclic voltammetry was recorded at the platinum tip electrode (potential range was chosen ac- 
cording to the expected range between hydrogen and oxygen evolution potentials, v = 0.05 Vs!) 
while the ODC was step-wise biased to a potential in between 0 and —1 V against its previously 
determined OCP. The activity of hydroxide and water was determined as described elsewhere 
[34]. In a short summary: By detection of cyclic voltammograms at a Pt nanoelectrode in the 
diffusion field above the ODC surface the Platinum oxide reduction peak was monitored. Con- 
sidering that the potential of this reaction is dependent on the activity of hydroxide and water, the 
activity quotient of both species can be extracted applying the Nernst equation. Furthermore, a 
comparison to a previously recorded calibration curve in different concentrated NaOH solutions 
enables the direct evaluation of the NaOH concentration. 
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D.1 Experimental 


D.1.1 Electrode preparation 


The preparation and characterization have been done by Barbara Ellendorff and David Franzen, 
Clausthal University of Technology, Germany. Section partly written by David Franzen. 

The preparation method and physical characterization for the examined electrodes is given in 
detail in a previous publication [23]. In short: For each electrode, a suspension containing sil- 
ver particles (SFOED, Ames Advanced Materials Corp.), a methyl cellulose solution (1 wt-% 
WALOCEL™ MKX 70000 PP Ol) as pore building agent and thickener, demineralized water 
and a PTFE dispersion (TF 5060GZ, 3M™ Dyneon™) in the desired ratio was used. The sus- 
pension was supplied with a spraying piston (Evolution, 0.6 mm pin hole, Harder & Steenbeck) 
on a conductive support (nickel mesh, 106 um x 118 um mesh size, 63 um thickness, Haver & 
Boecker OHG) with a total of 80 layers. In the further production steps the electrodes were hot 
pressed (LaboPress P200S, Vogt, 15 MPa, 130?C, 5 min) and sintered in an air oven (330?C, 
15 min) to improve the mechanical stability and to form the pore system by burning out the 
methyl cellulose completely. Detailed characteristics of the ODC are also given in a previous 


publication [23]. 


D.1.2 Electrochemical characterization 


The measurements have been done by David Franzen, Clausthal University of Technology, Ger- 
many. Section partly written by David Franzen. 

Electrochemical experiments were performed using a half-cell set up (FlexCell HZ PPO1, Gaska- 
tel GmbH) with a geometrical cell area of 3.14 cm? in a three electrode configuration. A platinum 


15 Parts of this chapter have been published in Róhe et al., Electrochimica acta, 2021 [15]. 
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wire acts as a counter electrode while the potentials were measured with a reversible hydrogen 
electrode (RHE) through a Luggin capillary in front of the GDE. 

The half-cell consists of an electrolyte and gas compartment divided by the GDE. As elec- 
trolyte approximately 30 mL of a 30 wt-% NaOH solution was prepared from caustic flakes 
(< 99 wt—%, Carl Roth) and demineralized water. In the gas compartment pure oxygen with a 
flowrate of 50 mLN min”? and a small back pressure using a 1 mm water column was applied. 
The whole half-cell was heated up to 80 ?C to correspond to industrial conditions [6]. The ex- 
periments were carried out with a Zennium Pro Potentiostat (Zahner GmbH) following the same 
routine for each electrode. In a startup procedure current densities of 1, 2, 3 and 4 kA m ? were 
applied for conditioning the electrodes. Afterwards pseudo-galvanostatic impedence measure- 
ments were performed at current densities of 2 kA m"?. A 10 mV amplitude was applied using a 
single sine with a starting frequency of 1 Hz going up to the upper limit of 50 KHz and sweeping 
down to the lower limit of 100 mHz. Sampling the frequencies between in the high frequency 
range twice is a standard proceed to stabilize the measurement. Afterwards linear sweep voltam- 
metry starting at open cell potential (OCP) to 200 mV vs. reversible hydrogen electrode (RHE) 
with a scan rate of 0.5 mV s^! was performed to obtain the polarization curves. 


D.2 EIS at different current densities 


I validated the model additional against EIS spectra of the ODC with 98 wt—% Ag (ODC with 
highest performance) at current densities of 1kAm~* and 4kA m ?. As presented earlier, the 
double layer capacitance in ODCs is dependent on the cell potential or the current density re- 
spectively [72]. Since this effect is important for the comparison of the impedance spectra at 
different current densities, the double-layer capacitances for the simulations with 1 and 4AKA m? 
have been adjusted according to the in [72] measured relations, all further simulation parameters 
have not been changed: 


e Caca(j = IkAm 2) © 1.3 -Cacal j = 2kAm ?) 
* Cacal j =4kAm?) z:0.6- Carca ( j = 2kA m~?) 


Measurement and simulation results for different current densities are shown in fig. D.1. The 
simulation data agree well with the experimental impedance spectra in term of shape, time con- 
stant and charge transfer resistance. Both, the measurements as well as the simulations, indicate 
a acceleration of the time constant with increasing current density. The simulation data show that 
at all current densities the time constant is dominated by double layer charge and discharge. 
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Figure D.1: Simulated and experimental EIS at different current densities of ODC with 98 wt—% Ag. 


D.3 Model parameters of the scenario analysis 


Please note that the discussed scenarios are hypothetical and based on divergent assumptions. 
The reaction rate (eq. (6.6)) of the ORR and the charge balance (eq. (6.13)) are functions of 
the wetted catalyst area. Since the amount of flooded pores and thus the wetted catalyst area 
are varied in the scenarios, the double layer capacitance and the reaction rate constant have been 
readjusted to reproduce the measured current densities and time constants. 

Semi-homogeneous 


* Number of domains: 5 


* Specific domain area: Constant for all domains, BE usi = const. = 20% 


Sphase interface 
AZ3PA 


* Specific gas-liquid interface: = const. = 6.4 um? um ? 


e Az3pa and Led according to table 6.2. 
Homogeneous 
* Number of domains: 1 


Sphase interface 


* Specific gas-liquid interface: Anp, — const. = 1.06 um? um”?, the three-phase area 
is located directly at the liquid side of the ODC with AzapA = 34 um. 


All further parameters are given in table D.1. 
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Table D.1: Parameters identified from electrochemical experiments and ODC model for ODC with 98 wt—% Ag with 
hypothetical scenarios. 


Name Symbol Unit Semi- Homogeneous 
homogeneous 

Reaction rate constant ko.cat molm-?s-! 0.28 0.51 

Double layer capacitance Cál.cat Fm? 3.04 5:2] 

Porosity“ £ — 0.4 

Specific catalyst interface“ Scat um? um”? 0.79 

Specific resistance“? R specific Qm? 2.83 x 10? 


@ Same parameter value as identified for the inhomogeneous modeled ODC. ? High frequency resistance 
from EIS measurements [23]. 
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